Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,660,636 C→A 27.7% intergenic (‑91/+61) epd ← / ← tktA D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,660,6360CA27.7% 505.7 / 238.6 300intergenic (‑91/+61)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base C (107/108);  new base A (47/36);  total (156/144)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 306 total aligned reads displayed.

ACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACT  >  NC_002947/5660487‑5660785
                                                                                                                                                          |                                                                                                                                                      
acacAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGc                                                                                                                                                        <  1:6385807/150‑1 (MQ=255)
acacAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGc                                                                                                                                                        <  1:2025380/150‑1 (MQ=255)
   caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGaaa                                                                                                                                                      <  1:2544857/150‑1 (MQ=255)
    aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCCCAGGGGTTGATGCAAcc                                                                                                                                                    >  1:36648/1‑150 (MQ=255)
     gCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACct                                                                                                                                                   <  2:2965491/150‑1 (MQ=255)
       ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccct                                                                                                                                                  <  2:6926557/150‑3 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt                                                                                                                                                <  2:466804/150‑5 (MQ=38)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt                                                                                                                                                <  2:5600782/150‑5 (MQ=38)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  2:8202950/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  2:7386491/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  1:1053311/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  1:8408598/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  2:1119135/150‑1 (MQ=255)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                              <  1:6966463/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                              <  1:4744189/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                              <  2:7102498/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                              <  2:4998777/150‑1 (MQ=255)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                              <  2:7314824/150‑1 (MQ=255)
           gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCATGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                             <  2:5951466/150‑1 (MQ=255)
             tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                           <  2:2925202/150‑1 (MQ=38)
             tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                           <  2:6515643/150‑1 (MQ=255)
                ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCg                                                                                                                                        >  2:5805292/1‑150 (MQ=255)
                 gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                                       <  1:8376544/150‑1 (MQ=255)
                 gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGACAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCgg                                                                                                                                       <  2:4150522/150‑1 (MQ=37)
                  ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCggg                                                                                                                                      <  1:5200708/150‑1 (MQ=255)
                   cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGc                                                                                                                                         <  2:3400038/150‑1 (MQ=21)
                      aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTg                                                                                                                                  >  2:7983408/1‑150 (MQ=255)
                      aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCAGTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTGGGGGGGAGCGGGCTTg                                                                                                                                  >  2:685143/1‑150 (MQ=255)
                               gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGact                                                                                                                         >  2:8371035/1‑148 (MQ=255)
                                tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCTCCCACACGGTTTAATGCAACCACAGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                        >  1:1685636/1‑150 (MQ=39)
                                  cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGaa                                                                                                                     >  1:4970369/1‑150 (MQ=255)
                                   aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGGGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAg                                                                                                                    >  2:1960883/1‑150 (MQ=255)
                                      tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCGCCCACACGGTTTCATGCAACCACGGTGGGAGCGGGCTTGCCCGCGCATGTAGGCg                                                                                                                 >  1:5359842/1‑150 (MQ=39)
                                             ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAATCCCGCACCCCAACGGTTTGATTAAAAATATGTGGGAGCGGGCTTGCCCCCGGAAGAAAGCGACAcccc                                                                                                          >  1:297180/1‑148 (MQ=255)
                                               aCGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgccc                                                                                                        >  2:2132040/1‑150 (MQ=255)
                                                   tGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                                    >  2:2845226/1‑150 (MQ=255)
                                                       gaatgcGGGGCTCACGAGACAGTCCAATTGTCGACAGCCTGCGGAATTGGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTGGCCCGCGAATGAAGGCGACACCGCCCGCAGATCt                                                                                                <  1:7594957/145‑1 (MQ=38)
                                                       gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCGCCCCCAGGGTTTAATGCAACCGCTGTGGGAGCGGGCTTGCCCGCGCATGAAGGCGACACCGCCCGCAGATCt                                                                                                >  1:1214670/1‑150 (MQ=39)
                                                                  tCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAACCCCCCGCCCCCCCGGGTTCAAGCAACCGCGGGGGGAGCGGGCTTGCCCCCGCAAGAATGCCACCCCCCCCGCCGATCGTGAAAGCACTc                                                                                     >  1:6091155/1‑150 (MQ=255)
                                                                    aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCAACCACAGGGTTTCATGCAACCGCTGTGGGATCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTGTAAATCAGTCtt                                                                                   >  2:809185/1‑150 (MQ=255)
                                                                      gagaCAGTTCAATTGTCGTCAGCCTGCGTAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTcc                                                                                 <  1:1100273/150‑1 (MQ=255)
                                                                        gaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCCTCcag                                                                               <  1:3499510/150‑1 (MQ=255)
                                                                          ctgTTCAATTGTCGTCCGCCTGCGGAATTTGGGGTGTGG‑‑‑‑CTTTCGCGGGCAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGGGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagca                                                                             <  1:791467/148‑1 (MQ=38)
                                                                           aGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTAATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagcag                                                                            >  2:1874272/1‑150 (MQ=255)
                                                                             ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGGTTCATGCACCCTCGGTGGGCGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGct                                                                          >  2:5493507/1‑150 (MQ=255)
                                                                                aaTTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGGTTCATGCAACCACGGTGGGAGCGGGCGTGCCCCCGAATGAAGGCGACCCCGCCCGCACATCATTAAATCAGTCTTCCAGCACGTCtt                                                                       >  1:579203/1‑150 (MQ=37)
                                                                                 aTTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTc                                                                      <  1:1959906/150‑1 (MQ=255)
                                                                                 aTTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTc                                                                      <  1:2615586/150‑1 (MQ=255)
                                                                                            ccTGCGTATTTTGCGGTGTGG‑‑‑‑CTTTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                           <  1:5164901/150‑1 (MQ=255)
                                                                                            ccTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                           >  2:4906436/1‑150 (MQ=255)
                                                                                            ccTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                           <  2:5135264/150‑1 (MQ=255)
                                                                                             cTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                          <  2:7333227/150‑1 (MQ=255)
                                                                                             cTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                          <  2:1537375/150‑1 (MQ=255)
                                                                                             cTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATCAAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                          <  1:4457680/150‑1 (MQ=255)
                                                                                              tGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAgg                                                         >  2:595991/1‑150 (MQ=255)
                                                                                              tGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAgg                                                         <  2:2244781/150‑1 (MQ=255)
                                                                                                   aaTTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCAGGCAACCTCGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGtt                                                    <  2:2506843/150‑1 (MQ=255)
                                                                                                        gCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTGATGCAACCTGTGGGGGAGCGGGGTTGGCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTGTTGCAGCAGCTCTTAGGGAGTACCGAGGATGTTATCCa                                               >  2:3321247/1‑150 (MQ=38)
                                                                                                         cGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAg                                              <  2:7701427/150‑1 (MQ=255)
                                                                                                         cGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAg                                              <  2:3195608/150‑1 (MQ=255)
                                                                                                             gggg‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTg                                          <  1:522281/148‑1 (MQ=255)
                                                                                                                g‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCCGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGCTGTTCTCCAGGGTGAAg                                       >  2:4744524/1‑150 (MQ=255)
                                                                                                                   ggaaTTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATCAAACCTCTGTGGGAGCGTGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                                        <  1:6691644/147‑1 (MQ=255)
                                                                                                                     cATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGc                                      <  2:4509616/150‑1 (MQ=255)
                                                                                                                               cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                            >  2:384005/1‑150 (MQ=255)
                                                                                                                               cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                            <  1:7025535/150‑1 (MQ=255)
                                                                                                                               cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                            <  1:4279615/150‑1 (MQ=255)
                                                                                                                               cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                            <  1:2877244/150‑1 (MQ=255)
                                                                                                                                  gCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                         >  2:5672621/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                       >  1:2828218/1‑150 (MQ=37)
                                                                                                                                    ccGCTCCCACAGAGGATGCATGAAACACTGTGGCAGCGGGCGTGCCCGCGAAAGAAGGCGACACCGAACCAAGATCATAAAATCAGTCTTCCAGAACCTCTTAGGCAGTACAGAGGATGTTCTCCAGGGTGAAGCCAAACACTTCGAACa                       >  2:6552270/1‑150 (MQ=255)
                                                                                                                                     cGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                       >  1:5015867/1‑150 (MQ=17)
                                                                                                                                      gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt                     <  2:8335951/150‑1 (MQ=255)
                                                                                                                                      gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt                     >  2:5741291/1‑150 (MQ=255)
                                                                                                                                         cccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc                  >  1:2360605/1‑150 (MQ=38)
                                                                                                                                          ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                 >  1:8240071/1‑150 (MQ=255)
                                                                                                                                          ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                 >  1:406606/1‑150 (MQ=255)
                                                                                                                                           cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa                >  1:5978662/7‑150 (MQ=38)
                                                                                                                                           cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAAT‑‑‑GCCACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGcc            >  1:1453219/7‑150 (MQ=25)
                                                                                                                                           cacaGGGTATCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa                >  2:4419735/1‑150 (MQ=255)
                                                                                                                                            acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  1:1284614/6‑150 (MQ=38)
                                                                                                                                            acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  1:1827566/6‑150 (MQ=38)
                                                                                                                                            acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  2:6978739/6‑150 (MQ=38)
                                                                                                                                            acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  1:7742965/6‑150 (MQ=38)
                                                                                                                                            acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  2:6287290/6‑150 (MQ=38)
                                                                                                                                            acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  2:8390821/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  1:5621253/1‑150 (MQ=255)
                                                                                                                                            acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa               >  1:4234305/1‑150 (MQ=255)
                                                                                                                                             caGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAg              >  1:5568528/1‑150 (MQ=255)
                                                                                                                                              aGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc             >  1:6274516/1‑150 (MQ=255)
                                                                                                                                                 ggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGTAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCg           <  1:361445/150‑1 (MQ=255)
                                                                                                                                                  gTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc         >  2:811637/1‑150 (MQ=255)
                                                                                                                                                   tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc         >  1:854614/1‑150 (MQ=255)
                                                                                                                                                       aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg     >  2:1406661/1‑150 (MQ=255)
                                                                                                                                                       aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg     >  1:6519331/1‑150 (MQ=255)
                                                                                                                                                       aTGCAACCTCTGTGGGAGCGGGCTAGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAATCCGGCGCCg     >  1:810359/1‑150 (MQ=255)
                                                                                                                                                        tGCACCCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAGGCCGAAGCCGGCGCCGa    <  2:6091155/150‑1 (MQ=255)
                                                                                                                                                        tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGa    <  1:7389692/150‑1 (MQ=255)
                                                                                                                                                         gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc   <  2:3626389/150‑1 (MQ=255)
                                                                                                                                                          cAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACt  >  2:8811845/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
ACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACT  >  NC_002947/5660487‑5660785

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