Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,770,234 T→G 5.2% S136S (TCT→TCG)  PP_1574 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,770,2340TG 5.2% 1188.6 / 36.2 346S136S (TCT→TCG) PP_1574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (188/139);  new base G (10/8);  total (198/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 367 total aligned reads displayed.

CCGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGC  >  NC_002947/1770085‑1770383
                                                                                                                                                     |                                                                                                                                                     
ccGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGctct                                                                                                                                                       >  1:7291370/1‑150 (MQ=255)
    ttGACCCGGTTGAACTGGAGAAGGACCCGGAGAGTACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCa                                                                                                                                                   >  2:7112450/1‑150 (MQ=255)
         ccGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGc                                                                                                                                              >  1:4769494/1‑150 (MQ=255)
         ccGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGc                                                                                                                                              >  1:4459606/1‑150 (MQ=255)
          cGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCt                                                                                                                                             >  2:6371647/1‑150 (MQ=255)
           ggTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTc                                                                                                                                            >  2:597794/1‑150 (MQ=255)
              tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTc                                                                                                                                         <  2:8478862/150‑1 (MQ=255)
              tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTc                                                                                                                                         <  1:6189951/150‑1 (MQ=255)
              tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTc                                                                                                                                         <  2:2453327/150‑1 (MQ=255)
                aaCTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTcga                                                                                                                                       <  1:7662101/150‑3 (MQ=255)
                  cTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTcgagc                                                                                                                                     <  1:5850478/150‑5 (MQ=255)
                   tGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCAGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAATACGAAGAGATCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCtctc                                                                                                                                    >  1:2888896/1‑150 (MQ=255)
                            aCCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATcgg                                                                                                                           <  2:4292420/150‑1 (MQ=255)
                            aCCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATcgg                                                                                                                           <  1:7169922/150‑1 (MQ=255)
                             cccGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATcggc                                                                                                                          >  2:3001182/1‑150 (MQ=255)
                                ggAGAGCACCACGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCggcggc                                                                                                                       <  2:8339927/150‑1 (MQ=255)
                                ggAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCggcggc                                                                                                                       <  1:5505183/150‑1 (MQ=255)
                                      cACCAAGCGCCTGGGATCTCTTGATGGAGACCTGTCCGTCTCGATGCTGGTAGCTATTCGGGAGGTCAACGACATGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGATCTCGTACTCTCCGGCATCGGCGGCTcgacg                                                                                                                 >  2:3731657/1‑150 (MQ=255)
                                       aCCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgacga                                                                                                                >  1:2789666/1‑150 (MQ=255)
                                       aCCAAGCGCCTGGGATCTCTTGATGGAGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACTAATAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCAGGCATCGGCGGCTcgacga                                                                                                                >  1:8793457/1‑150 (MQ=255)
                                        ccAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGAt                                                                                                               <  2:7357293/150‑1 (MQ=255)
                                        ccAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGAt                                                                                                               <  2:3096235/150‑1 (MQ=255)
                                            gcgcCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcgcg                                                                                                           >  1:5381289/1‑150 (MQ=255)
                                                cTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAgg                                                                                                       >  2:1489989/1‑150 (MQ=255)
                                                cTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAgg                                                                                                       <  2:8137916/150‑1 (MQ=255)
                                                   ggATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCa                                                                                                    >  2:2755969/1‑150 (MQ=255)
                                                       ctctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGCAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGa                                                                                                <  2:2818606/150‑1 (MQ=255)
                                                        tctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaa                                                                                               >  1:5174699/1‑150 (MQ=255)
                                                        tctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaa                                                                                               >  1:3041839/1‑150 (MQ=255)
                                                          ttGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaaa                                                                                             >  2:3311794/1‑150 (MQ=255)
                                                          tagaaggCGACCAGTCCGTCGCGCTGCTGGTACCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaaa                                                                                             >  2:5816654/6‑150 (MQ=255)
                                                           tGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAc                                                                                            >  2:3924657/1‑150 (MQ=255)
                                                               ggCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                        <  2:7533364/150‑1 (MQ=255)
                                                               ggCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                        <  2:1345413/150‑1 (MQ=255)
                                                               ggCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                        <  1:4613154/150‑1 (MQ=255)
                                                                     cTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGc                                                                                  >  2:7844786/1‑150 (MQ=255)
                                                                       gTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGt                                                                                >  2:1170350/1‑150 (MQ=255)
                                                                        tCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTg                                                                               >  1:3395974/1‑150 (MQ=255)
                                                                            tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGtt                                                                           >  1:1950584/1‑150 (MQ=255)
                                                                              gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                         >  1:2364489/1‑150 (MQ=255)
                                                                              gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                         >  2:183950/1‑150 (MQ=255)
                                                                              gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                         <  2:4753523/150‑1 (MQ=255)
                                                                               cgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                        <  1:6281225/150‑1 (MQ=255)
                                                                                gctgctGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAgg                                                                       <  1:2453237/150‑1 (MQ=255)
                                                                                 ctgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGt                                                                      >  1:6507237/1‑150 (MQ=255)
                                                                                  tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTc                                                                     >  2:8405599/1‑150 (MQ=255)
                                                                                     tgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCCTCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTg                                                                  <  2:8016801/150‑1 (MQ=255)
                                                                                     tgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTg                                                                  >  2:4670762/1‑150 (MQ=255)
                                                                                      ggTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTgg                                                                 <  1:6757450/150‑1 (MQ=255)
                                                                                          gCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgaga                                                             <  2:1987660/150‑1 (MQ=255)
                                                                                           cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                            <  1:8165684/150‑1 (MQ=255)
                                                                                           cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                            <  1:1848782/150‑1 (MQ=255)
                                                                                           cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                            <  1:2210143/150‑1 (MQ=255)
                                                                                              ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                         >  1:103777/1‑150 (MQ=255)
                                                                                               tCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTAc                                                        >  2:2961015/1‑150 (MQ=255)
                                                                                                cGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACa                                                       <  2:5567942/150‑1 (MQ=255)
                                                                                                 gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                      <  2:4075899/150‑1 (MQ=255)
                                                                                                 gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                      >  2:6941096/1‑150 (MQ=255)
                                                                                                  ggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAgg                                                     <  1:6915258/150‑1 (MQ=255)
                                                                                                  ggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAgg                                                     <  1:268723/150‑1 (MQ=255)
                                                                                                    aGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCa                                                   >  1:745559/1‑150 (MQ=255)
                                                                                                         aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                        >  2:3656195/1‑124 (MQ=255)
                                                                                                         aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                              >  1:4335821/1‑150 (MQ=255)
                                                                                                         aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                              >  1:5313564/1‑150 (MQ=255)
                                                                                                         aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                              >  1:7047216/1‑150 (MQ=255)
                                                                                                         aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                              >  1:7917684/1‑150 (MQ=255)
                                                                                                           cGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGt                                            <  1:3628473/150‑1 (MQ=255)
                                                                                                               cTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCt                                        >  2:8287467/1‑150 (MQ=255)
                                                                                                               cTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCt                                        >  1:1672261/1‑150 (MQ=255)
                                                                                                                 gggAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCg                                      >  2:5438097/1‑150 (MQ=255)
                                                                                                                    aaaGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTg                                   <  1:6096627/150‑1 (MQ=255)
                                                                                                                      aGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCa                                 >  2:3839260/1‑150 (MQ=255)
                                                                                                                       gACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCAt                                <  1:1332449/150‑1 (MQ=255)
                                                                                                                        aCGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATg                               >  2:5120824/1‑150 (MQ=255)
                                                                                                                        aCGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATg                               >  2:6593515/1‑150 (MQ=255)
                                                                                                                          gAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                >  2:5726382/1‑115 (MQ=255)
                                                                                                                          gAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                <  1:5726382/115‑1 (MQ=255)
                                                                                                                           aaGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCg                            <  1:3313927/150‑1 (MQ=255)
                                                                                                                             gagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCCGGCCAAAGAAAACCTGAGCCTGCGTGAGTCCAGGTCTGGGGAGAAGTACAGGCAGCACCGAGGCTCGCTGCATGTCggg                          >  2:6187128/1‑150 (MQ=255)
                                                                                                                             gagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCggg                          <  1:3744609/150‑1 (MQ=255)
                                                                                                                              agCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGa                         <  2:2120192/150‑1 (MQ=255)
                                                                                                                                 tCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgc                      >  1:3929195/1‑150 (MQ=255)
                                                                                                                                   aCCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgcgc                    >  2:2208574/1‑150 (MQ=255)
                                                                                                                                   aCCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgcgc                    >  1:1783782/1‑150 (MQ=255)
                                                                                                                                    ccGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGcgcg                   >  2:8069171/1‑150 (MQ=255)
                                                                                                                                      gAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCCTCGGCGGCTCGACGATCGCGCCGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGAGGCGCGTg                 <  2:864747/150‑1 (MQ=255)
                                                                                                                                      gAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTg                 >  2:5289212/1‑150 (MQ=255)
                                                                                                                                      gAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTg                 <  2:6467114/150‑1 (MQ=255)
                                                                                                                                        gCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGa               <  1:8285859/150‑1 (MQ=255)
                                                                                                                                         cTCGACGAGCTCTTGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCAAAATAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCTTGGCTCGCTGCATGTCGGGATGCGCGTGGat              <  2:6745078/150‑2 (MQ=255)
                                                                                                                                         cTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAg              <  1:2291778/150‑1 (MQ=255)
                                                                                                                                          tCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgc             >  1:2589389/1‑150 (MQ=255)
                                                                                                                                          tCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgc             >  1:8613059/1‑150 (MQ=255)
                                                                                                                                            gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc           >  2:2638484/1‑150 (MQ=255)
                                                                                                                                            gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCAGGTCGGGATGCGCGTGGAgcgc           >  2:8739882/1‑150 (MQ=255)
                                                                                                                                             acgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcg          <  1:5540199/150‑1 (MQ=255)
                                                                                                                                             acgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcg          >  1:6736423/1‑150 (MQ=255)
                                                                                                                                              cgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcgc         <  1:505085/150‑1 (MQ=255)
                                                                                                                                              cgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcgc         >  2:2356228/1‑150 (MQ=255)
                                                                                                                                                     tGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGc  <  1:8373981/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGC  >  NC_002947/1770085‑1770383

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