Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,362,263 C→G 6.6% V49L (GTG→CTG)  tpiA ← triose phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,362,2630CG 6.6% 276.7 / 10.0 91V49L (GTG→CTG) tpiAtriose phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base G (3/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

CGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATG  >  NC_002947/5362120‑5362408
                                                                                                                                               |                                                                                                                                                 
cGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTCATCGAGAGTCCCTTACTCGAACCCGGCATGGCGCTGACGGTGGACAGGAGGATTTTGTTTGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgatt                                                                                                                                             <  2:1463690/150‑1 (MQ=255)
 gAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattg                                                                                                                                            >  1:2422347/1‑150 (MQ=255)
  aaTGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGGTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattga                                                                                                                                           <  1:1209609/150‑1 (MQ=255)
        cAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGTCATCAATCAGGTGATtgatggcct                                                                                                                                     <  2:1634058/150‑5 (MQ=255)
         aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                         <  1:1010705/145‑1 (MQ=255)
         aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                         >  2:1010705/1‑145 (MQ=255)
         aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACaa                                                                                                                                    >  1:215351/1‑150 (MQ=255)
         aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACaa                                                                                                                                    >  1:2228523/1‑150 (MQ=255)
         aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACaa                                                                                                                                    >  1:2228268/1‑150 (MQ=255)
           tCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTAGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAgg                                                                                                                                  >  1:441481/1‑150 (MQ=255)
           tCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGAAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAgg                                                                                                                                  >  2:801796/1‑150 (MQ=255)
            cAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGc                                                                                                                                 <  2:1649925/150‑1 (MQ=255)
                  cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                           <  2:215351/150‑1 (MQ=255)
                  cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg                                                                                                                           >  1:951239/1‑150 (MQ=255)
                   aaCTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGa                                                                                                                          >  2:2258864/1‑150 (MQ=255)
                     cTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGa                                                                                                                          <  1:579981/148‑1 (MQ=255)
                     cTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGa                                                                                                                          >  2:579981/1‑148 (MQ=255)
                      ttGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAAt                                                                                                                       >  1:814560/1‑150 (MQ=255)
                        gCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATAc                                                                                                                     <  2:2123266/150‑1 (MQ=255)
                          aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatggc                                                                                                                                       <  1:236358/130‑3 (MQ=255)
                          aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatggc                                                                                                                                       >  2:236358/1‑128 (MQ=255)
                          aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACgg                                                                                                                   >  1:390031/1‑150 (MQ=255)
                          aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACgg                                                                                                                   <  2:960172/150‑1 (MQ=255)
                           aCCTGCTTCAACCAGCTGCTCGGGAGCCACTTACACGGTCTGTCCTCCCTTTTCGGGTTGTACAGCAGAATCCTGCGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  <  1:2377108/150‑1 (MQ=255)
                           aCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  >  1:2415205/1‑150 (MQ=255)
                           aCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  <  1:1427785/150‑1 (MQ=255)
                           aCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  <  1:801796/150‑1 (MQ=255)
                              tGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGCTTGATGAACAAGGCCGGTGGAAATACGGCGAc                                                                                                               <  2:121316/150‑1 (MQ=255)
                              tGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGGACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCACTCACCTGATTGATGAACAAGGCCGGGGGAAATACGGCGAc                                                                                                               <  2:945738/150‑1 (MQ=255)
                                cTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  <  2:1508331/145‑1 (MQ=255)
                                cTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                  >  1:1508331/1‑145 (MQ=255)
                                   cAACCCCCTGCCCCGGAGCAAATTCCCCGGTCAGCGCACCCTGTTCGTGTTGCACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAAGCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCa                                                                                                          <  2:2023809/150‑1 (MQ=255)
                                      ccAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTc                                                                                                       <  1:1807974/150‑1 (MQ=255)
                                      ccAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTc                                                                                                       >  2:1998877/1‑150 (MQ=255)
                                      ccAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTc                                                                                                       >  2:1998809/1‑150 (MQ=255)
                                         gctgctAGGGAGCAACTTCCCCGGTCAGCGCACCCTGTCCGGGCTGTACAGAAGACTTCTGTGAGCAGACAGTAATTTCTTTACATGCCAGGCCATCAATCACCTGATTCATGAACAAGGCCGGGGGAAATACGCCGACTCCCACTCCGc                                                                                                    >  2:988113/1‑150 (MQ=255)
                                         gctgctAGGGAGCAACTTCCCCGGCCAGAGAACCCCGTTCGGGATGTACAGAAGAATTATATGAACCGAAAGTAATTTTTTTTCCTGCAATTCCATCAAACACCTGATTAATAAACAAGGCCGGCGGAAATACGGCGACTTCCACTCCGc                                                                                                    >  2:988535/1‑150 (MQ=255)
                                                  gAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAgg                                                                                           <  2:2415205/150‑1 (MQ=255)
                                                   aGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                          <  1:1998809/150‑1 (MQ=255)
                                                   aGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                          <  1:1998877/150‑1 (MQ=255)
                                                      aaCTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCgag                                                                                       >  1:2183941/1‑150 (MQ=255)
                                                      aaCTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAAGACAGGGGGTAAACGGCATTTTCCACTCCGCTCGGCAAGGCgag                                                                                       >  2:1840728/1‑150 (MQ=255)
                                                          tCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTa                                                                                   <  1:988113/150‑1 (MQ=255)
                                                          tCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTa                                                                                   <  1:988535/150‑1 (MQ=255)
                                                             ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGAttccgc                                                                                >  1:1140233/1‑146 (MQ=255)
                                                             ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                <  2:2422347/150‑1 (MQ=255)
                                                             ccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                >  2:1340520/1‑150 (MQ=255)
                                                              cGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                          >  2:799008/1‑139 (MQ=255)
                                                              cGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                          <  1:799008/139‑1 (MQ=255)
                                                              cGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGAGAGATTACTcc                                                                               >  2:2086611/1‑149 (MQ=255)
                                                                gTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGGAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAg                                                                             <  2:1079373/150‑1 (MQ=255)
                                                                    gcgcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            >  1:1709385/1‑115 (MQ=255)
                                                                    gcgcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTc                                                                                                            <  2:1709385/115‑1 (MQ=255)
                                                                         cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCCAGCCTTTGGtc                                                                    <  2:33299/150‑1 (MQ=255)
                                                                                 ggtgggTACAGCAGAAGTCTGTGCGCCGACAGTAATTTCTTGACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagctcagc                                                            <  2:814560/145‑1 (MQ=255)
                                                                                    ttGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGAc                                                         >  2:2368427/1‑150 (MQ=255)
                                                                                       tACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCt                                                      >  2:790863/1‑150 (MQ=255)
                                                                                       tACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCt                                                      >  2:2056030/1‑150 (MQ=255)
                                                                                             agAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGcgcg                                                >  2:1713663/1‑150 (MQ=255)
                                                                                              gAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg                                               <  1:1340520/150‑1 (MQ=255)
                                                                                               aaTTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCggg                                              >  1:1792391/1‑150 (MQ=255)
                                                                                                 ttCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTa                                            >  2:1928583/1‑150 (MQ=255)
                                                                                                          ccGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAtt                                   <  2:1469803/150‑1 (MQ=255)
                                                                                                          ccGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAtt                                   >  1:1469803/1‑150 (MQ=255)
                                                                                                             aCAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTc                                <  1:396419/150‑1 (MQ=255)
                                                                                                               aGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                              <  1:790863/150‑1 (MQ=255)
                                                                                                               aGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCAACTCCGCTCGGCAAGGCGGGATTACTCAAGCCTTTGGTCAGCCCAGCGACGCTAGCGCGGGTACCGTGCATGTTCCa                              <  2:441481/150‑1 (MQ=255)
                                                                                                                 tAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGt                            >  1:620468/1‑150 (MQ=255)
                                                                                                                  aaTTTCTTTACCTGCCAGGCCATCAATCGCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGtt                           >  1:1457495/1‑150 (MQ=255)
                                                                                                                  aaTTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCTACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGtt                           <  2:2228523/150‑1 (MQ=255)
                                                                                                                  aaTTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGtt                           <  2:2228268/150‑1 (MQ=255)
                                                                                                                    tttCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAgg                                                                                                                                  >  1:586681/1‑45 (MQ=255)
                                                                                                                    tttCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAgg                                                                                                                                  <  2:586681/45‑1 (MQ=255)
                                                                                                                    tttCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTAc                         >  2:2006738/1‑150 (MQ=255)
                                                                                                                       cTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAgcgc                                                 >  1:908281/1‑123 (MQ=255)
                                                                                                                       cTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAgcgc                                                 <  2:908281/123‑1 (MQ=255)
                                                                                                                        tttACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGc                     >  1:900466/1‑150 (MQ=255)
                                                                                                                            ccTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTAcc                 >  1:661461/1‑150 (MQ=255)
                                                                                                                             cTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCa                >  2:16017/1‑150 (MQ=255)
                                                                                                                                gcgggCCATCAATCGCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAg             <  2:951239/147‑1 (MQ=255)
                                                                                                                                 cAGGCCATGCAGCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCTAGCCTCTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAgg            <  2:1140233/150‑1 (MQ=255)
                                                                                                                                 cAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAgg            <  2:900466/150‑1 (MQ=255)
                                                                                                                                  aGGCCATCAATCACCTGATTGATGAATAAGGCCGGGGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAggg           >  1:564950/1‑150 (MQ=255)
                                                                                                                                  aGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAggg           <  1:77413/150‑1 (MQ=255)
                                                                                                                                    gCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCg         >  2:844903/1‑150 (MQ=255)
                                                                                                                                      cATCAATCACCTCACTGATGCACAAGGCCGGTGGAAATACGGCGACTTCCGCTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCATGGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc       <  2:2169447/150‑1 (MQ=255)
                                                                                                                                       atcaatcaGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACg      >  2:715773/1‑150 (MQ=255)
                                                                                                                                       atcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTAc                  >  1:220056/1‑138 (MQ=255)
                                                                                                                                       atcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTAc                  <  2:220056/138‑1 (MQ=255)
                                                                                                                                        tcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGc     >  2:760640/1‑150 (MQ=255)
                                                                                                                                           atcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATg  >  1:707660/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATG  >  NC_002947/5362120‑5362408

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