Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 994,761 (C)11→9 intergenic (+224/‑170) unknown → / → unknown Magnesium and cobalt efflux protein CorC/Apolipoprotein N‑acyltransferase / Copper homeostasis protein CutE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1994,7520C.100.0% 172.0 / NA 42intergenic (+215/‑180)unknown/unknownMagnesium and cobalt efflux protein CorC/Apolipoprotein N‑acyltransferase / Copper homeostasis protein CutE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (24/18);  total (24/18)
*1994,7530C.97.6% 165.1 / ‑1.7 42intergenic (+216/‑179)unknown/unknownMagnesium and cobalt efflux protein CorC/Apolipoprotein N‑acyltransferase / Copper homeostasis protein CutE
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (23/18);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

GAGGAAGAGGTTGTGGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTCCCCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAATATTTCACTTATCTTATCGCACTCATTTTCGGCGCTGCCGGCGT  >  1/994517‑994984
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                       
gaggaagaggTTGTGGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcaggg                                                                                                                                                                                                                            >  1:825179/1‑244 (MQ=255)
 aggaagaggTTGTGGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:529525/250‑8 (MQ=255)
 aggaagaggTTGTGGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:724638/250‑8 (MQ=255)
            gtgGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:560273/239‑8 (MQ=255)
             tgGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTc                                                                                                                                                                                                             <  1:480671/250‑1 (MQ=255)
                cGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTccccc                                                                                                                                                                                                         <  1:882293/251‑1 (MQ=255)
                         ttcttTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGa                                                                                                                                                                                                <  1:590203/251‑1 (MQ=255)
                          tcttTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:693899/225‑8 (MQ=255)
                               ccccTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAg                                                                                                                                                                                          >  1:746834/1‑251 (MQ=255)
                                    ctTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:15764/215‑8 (MQ=255)
                                    ctTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAggg                                                                                                                                                                                             <  1:516605/243‑1 (MQ=255)
                                    ctTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAg                                                                                                                                                                                     >  1:288663/1‑251 (MQ=255)
                                    ctTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAAGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGGGCCAAg                                                                                                                                                                                     >  1:754041/1‑251 (MQ=255)
                                         gAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGttt                                                                                                                                                                                  >  1:730824/1‑249 (MQ=255)
                                                            gagTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAagt                                                                                                                                                               <  1:63233/249‑1 (MQ=255)
                                                            gagTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAag                                                                                                                                                                >  1:657752/1‑248 (MQ=255)
                                                                 gAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTAc                                                                                                                                                         <  1:52326/250‑1 (MQ=255)
                                                                  aaCGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACaa                                                                                                                                                       <  1:553490/251‑1 (MQ=255)
                                                                    cGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATg                                                                                                                                                     <  1:41699/251‑1 (MQ=255)
                                                                    cGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATg                                                                                                                                                     <  1:81856/251‑1 (MQ=255)
                                                                    cGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTCCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAAc                                                                                                                                                                                                  <  1:425517/206‑1 (MQ=255)
                                                                            gTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:293278/175‑8 (MQ=255)
                                                                                      tctcCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGtt                                                                                                                                   >  1:314580/1‑251 (MQ=255)
                                                                                      tctcCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGtt                                                                                                                                   <  1:876185/251‑1 (MQ=255)
                                                                                         cccGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTc                                                                                                                                 >  1:351018/1‑250 (MQ=255)
                                                                                            gACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           <  1:642876/159‑8 (MQ=255)
                                                                                            gACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccctcagggc                                                                                                                                                                                                                           >  1:577318/1‑152 (MQ=255)
                                                                                                 caGCTGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTAc                                                                                                                                                         <  1:729028/218‑1 (MQ=255)
                                                                                                           aaaCATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAg                                                                                                                                           >  1:570037/1‑222 (MQ=255)
                                                                                                              cATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTGCGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCTCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAAc                                                                                                                                             <  1:40340/217‑1 (MQ=255)
                                                                                                              cATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAagt                                                                                                                                                               <  1:803778/199‑1 (MQ=255)
                                                                                                              cATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACGGCTGTGGtcgggaga                                                                                                                             >  1:894417/1‑226 (MQ=255)
                                                                                                               aTTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACGGCTGTGGtcgggag                                                                                                                              <  1:149800/231‑7 (MQ=255)
                                                                                                                ttATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGcc                                                                                                                                                                                        >  1:802416/1‑172 (MQ=255)
                                                                                                                ttATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAagta                                                                                                                                                              >  1:309034/1‑198 (MQ=255)
                                                                                                                  ataATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGtt                                                                                                                                   >  1:624899/1‑223 (MQ=255)
                                                                                                                  ataATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGtt                                                                                                                                   >  1:497779/1‑223 (MQ=255)
                                                                                                                   taATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTccccccccc                                                                                                                                                                                                                                  >  1:672173/1‑129 (MQ=255)
                                                                                                                    aaTTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGtt                                                                                                                                   >  1:140671/1‑221 (MQ=255)
                                                                                                                               atTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACGGCTGTGGtcgggag                                                                                                                              >  1:805056/1‑209 (MQ=255)
                                                                                                                                 ttttAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTc                                                                                                                                                                            >  1:387442/1‑167 (MQ=255)
                                                                                                                                              gCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACaa                                                                                                                                                       >  1:170041/1‑175 (MQ=255)
                                                                                                                                              gCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTgaga                                                                           >  1:3273/1‑251 (MQ=255)
                                                                                                                                                                     gTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCGCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATg                                                                                                                                                     <  1:233502/154‑1 (MQ=255)
                                                                                                                                                                     gTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAAt                                                                                    <  1:500474/219‑1 (MQ=255)
                                                                                                                                                                                          ggATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAATATTTCActtatct                                <  1:242776/250‑1 (MQ=255)
                                                                                                                                                                                           gATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑CCCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAata                                            >  1:207560/1‑238 (MQ=255)
                                                                                                                                                                                                aTCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAAt                                                                                    >  1:459356/1‑192 (MQ=255)
                                                                                                                                                                                                aTCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACgcgcgc                                                                 >  1:20837/1‑211 (MQ=255)
                                                                                                                                                                                                                       tGCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAATATTTCACTTATCTTATCGCACTCATTTTCGGCGCTGCCGGCgt  >  1:626301/1‑251 (MQ=255)
                                                                                                                                                                                                                        gCTTAAAAATTGAGATTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGc                                                                                                                      >  1:791368/1‑134 (MQ=255)
                                                                                                                                                                                                                                   gagaTTTT‑‑CCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAATATTTCActtatctt                               >  1:219079/1‑210 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                       
GAGGAAGAGGTTGTGGCGTAGTTTCTTCTTTCCCCTCTTTAGAAAAGAGGGGCGGGGGGAGAGTTGAACGATAACGGTGGCGCGCGTCTCCCGACCACAGCCGTTCGAAACATTATAATTCATTGATATTTTTAAATATCTTGCCTTCCTCCGTTTACGGGCAAGGTGCCCGAAGGGCGGAAGGGGGGATTTATCTTGATTTTTTCACAAGATCATGCTTAAAAATTGAGATTTTCCCCCCCCCCCTCAGGGCTGCGCCCAGCTCCCCCTGTGAACGAGGGGAGCCAAGTTTTTTTTCTTTTTATATAAGTAGTTACAATGTTTCGAACAGCTGTGGTTTCCCGACGCGTGCCGTTTATTTCATTAAATAACGAAAATTTAGGAATTGCGTGAGAAAACGCGCGCTATCCCGCTTATGAACAAATATTTCACTTATCTTATCGCACTCATTTTCGGCGCTGCCGGCGT  >  1/994517‑994984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑