Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,157,189 T→C 6.6% intergenic (‑21/+185) PP_5398 ← / ← PP_5399 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,157,1890TC 6.6% 285.4 / 6.4 94intergenic (‑21/+185)PP_5398/PP_5399transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/50);  new base C (3/3);  total (40/54)
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.64e-01

GCCGTGGCTGATGGCGACACTGGACACGGTTGCGCCAGGTTGCAGGCATTCCTGAACGACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTAC  >  NC_002947/6157067‑6157338
                                                                                                                          |                                                                                                                                                     
gagctcaccacgatgttggccagaagcaaaaagcttcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttccgaactccaGCATTAAGGGCGATGGTGTCCGCTTAAAAATACGCGGa                                                                                                                            <  1:3416603‑M2/38‑1 (MQ=255)
    tcaccacgatgttggccagaagcaaaaagcttcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttccgaactccaGCATTAAGGGCGATGGCGTCCGCTTAAAAATACGCGGACAcc                                                                                                                        >  2:144146‑M2/109‑150 (MQ=255)
       ccacgatgttggccagaagcaaaaagcttcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttccgaactccaGCATTAAGGGCGATGGCGTCCGCTTAAAAATACGCGGACACCATc                                                                                                                     <  1:144146‑M2/45‑1 (MQ=255)
         acgatgttggccagaagcaaaaagcttcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttccgaactccaGCATTAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGc                                                                                                                   <  1:2460503‑M2/47‑1 (MQ=255)
                    cagaagcaaaaagcttcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttcagaactccaGCATTAAGGGCGATGGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTg                                                                                                        <  2:1462826‑M2/58‑1 (MQ=255)
                                   tcctcgaataggcctgatacatagatgcaaccgggttcctatggtaagcgtccggcgaactcaccttccgaactccaGCATTAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGa                                                                                         <  2:1296437‑M2/73‑1 (MQ=255)
                                                      aaCGACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTc                                                                      <  1:3864039/150‑1 (MQ=31)
                                                         gACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCt                                                                   <  1:1219270/150‑1 (MQ=33)
                                                         gACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCt                                                                   <  1:3028393/150‑1 (MQ=33)
                                                         gACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGGTCCCCGGACGCTTACTCGCt                                                                   <  1:1560222/150‑1 (MQ=18)
                                                            cTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTa                                                                >  2:2238349/1‑150 (MQ=35)
                                                            cTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTa                                                                <  2:1832274/150‑1 (MQ=35)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  2:3556412/150‑1 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  2:433378/150‑1 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               >  1:4420635/1‑150 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  2:4773451/150‑1 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  2:814366/150‑1 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  1:1402460/150‑1 (MQ=37)
                                                             tGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  1:1965105/150‑1 (MQ=37)
                                                              gggCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATg                                                              >  1:2444968/1‑150 (MQ=38)
                                                               ggCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGc                                                             <  1:1702712/150‑1 (MQ=38)
                                                                 cTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGt                                                           >  1:2604744/1‑150 (MQ=40)
                                                                 cTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGt                                                           <  1:1706343/150‑1 (MQ=40)
                                                                  tttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGtt                                                          <  2:2017820/150‑1 (MQ=40)
                                                                  tttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGtt                                                          <  1:4901763/150‑1 (MQ=40)
                                                                  tttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGtt                                                          >  2:2884885/1‑150 (MQ=40)
                                                                  tttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGtt                                                          >  2:3153720/1‑150 (MQ=40)
                                                                   ttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTc                                                         <  2:3340616/150‑1 (MQ=255)
                                                                   ttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGGTCCCCGGACGCTTACTCGCTGTATGCGTTc                                                         <  1:939201/150‑1 (MQ=255)
                                                                   ttGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAAGCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTc                                                         <  2:3340541/150‑1 (MQ=255)
                                                                    tGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCg                                                        >  1:143901/1‑150 (MQ=255)
                                                                    tGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCg                                                        <  1:4936749/150‑1 (MQ=255)
                                                                     gatgggTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGGAATAAGGGTGATCGCGTCCGCTTAAAATTTCGCGGACACCATCGCCCTTATTGCTGGAGTTCGGATGGTGAGTTCGCCGGTCGCTTACTCGCTGTATGCGTTCgt                                                       >  2:4036847/5‑149 (MQ=14)
                                                                     gAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCga                                                       <  1:2666885/150‑1 (MQ=255)
                                                                     gAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCga                                                       >  2:4537793/1‑150 (MQ=255)
                                                                     cgatctattCGGATACGAGCTTCGTGGGCGCCTGGAAATCCTGGCAATAATGGAGACCGCGTCCGCTTCCAAATACGCGGACACCACCGCCCTAAAGGCTGGAGTTCGGAAGGTGAGTGCCCCGGACGCTTACTCGCTGTATGCGTTCga                                                       <  1:17831/143‑1 (MQ=14)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGt                                                           <  2:4054182/145‑1 (MQ=40)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGt                                                           >  1:4054182/1‑145 (MQ=40)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCgag                                                      >  1:2351871/1‑150 (MQ=255)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCgag                                                      >  2:1866370/1‑150 (MQ=255)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCgag                                                      >  1:262944/1‑150 (MQ=255)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCgag                                                      >  2:4836939/1‑150 (MQ=255)
                                                                      aaCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGAGTTCgag                                                      >  2:4836778/1‑150 (MQ=255)
                                                                       aCGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCgaga                                                     >  2:4130041/1‑150 (MQ=255)
                                                                           tttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACggg                                                 >  1:3265153/1‑150 (MQ=255)
                                                                           tttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACggg                                                 >  1:3831888/1‑150 (MQ=255)
                                                                            ttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGa                                                <  1:692639/150‑1 (MQ=255)
                                                                            ttCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGa                                                <  2:1666345/150‑1 (MQ=255)
                                                                                   aaGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGt                                         <  1:4678131/150‑1 (MQ=255)
                                                                                   aaGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGt                                         >  1:5420281/1‑150 (MQ=255)
                                                                                    agagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                                        <  1:694100/150‑1 (MQ=255)
                                                                                    agagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                                        <  1:1101474/150‑1 (MQ=255)
                                                                                    agagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGtt                                        <  2:4625572/150‑1 (MQ=255)
                                                                                     gagCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTg                                       >  2:752603/1‑150 (MQ=255)
                                                                                      agCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGt                                      <  1:4147270/150‑1 (MQ=255)
                                                                                      agCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGt                                      <  2:4157944/150‑1 (MQ=255)
                                                                                         ttCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACt                                   >  2:296163/1‑150 (MQ=255)
                                                                                             ttGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCt                               >  2:276397/1‑150 (MQ=255)
                                                                                               ggCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTgg                             <  1:4713743/150‑1 (MQ=255)
                                                                                                gcgcATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGa                            <  2:209910/150‑1 (MQ=255)
                                                                                                 cgcATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGAc                           >  2:2878832/1‑150 (MQ=255)
                                                                                                    aTGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCa                        >  1:859597/1‑150 (MQ=255)
                                                                                                      ggAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAgg                      >  2:4471731/1‑150 (MQ=255)
                                                                                                      ggAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAgg                      >  2:4532411/1‑150 (MQ=255)
                                                                                                       gAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAggg                     >  2:749527/1‑150 (MQ=255)
                                                                                                         aaTCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGaa                   <  1:1725705/150‑1 (MQ=255)
                                                                                                         aaTCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGaa                   <  1:1725686/150‑1 (MQ=255)
                                                                                                            ccTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               >  1:4039619/1‑103 (MQ=40)
                                                                                                            ccTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTAt                                                               <  2:4039619/103‑1 (MQ=40)
                                                                                                            ccTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATc                <  2:3449218/150‑1 (MQ=255)
                                                                                                             cTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTCCTGGACGCAGGGAAATCa               <  1:3948349/150‑1 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              <  2:1294830/150‑1 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              <  2:2785872/150‑1 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              >  1:5413434/1‑150 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              >  1:951347/1‑150 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              >  2:1101144/1‑150 (MQ=255)
                                                                                                              tGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGAAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              <  2:4421456/150‑1 (MQ=255)
                                                                                                              caGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATTGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAg              >  1:2326205‑M2/3‑150 (MQ=255)
                                                                                                               aGCATTAAGGGCGATGGTGTCCGCGTATTTTTAAGCGGACGCGATCACCCTTATTGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAgg             >  2:1516576‑M2/2‑150 (MQ=255)
                                                                                                                gCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGTGAAATCATGt            <  2:991581/150‑1 (MQ=255)
                                                                                                                  aaTAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGc          >  2:4341153/1‑150 (MQ=255)
                                                                                                                  aaTAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGc          >  2:2554031/1‑150 (MQ=255)
                                                                                                                  aaTAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGc          <  2:4607581/150‑1 (MQ=255)
                                                                                                                    tAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAg                       <  2:143527/135‑1 (MQ=255)
                                                                                                                    tAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAg                       >  1:143527/1‑135 (MQ=255)
                                                                                                                      aGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGaa      >  2:2539164/1‑150 (MQ=255)
                                                                                                                       gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGGGGAAg     <  2:262944/150‑1 (MQ=255)
                                                                                                                       gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAg     <  1:4461970/150‑1 (MQ=255)
                                                                                                                       gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAg     <  2:4258532/150‑1 (MQ=255)
                                                                                                                       gggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAg     <  1:2035752/150‑1 (MQ=255)
                                                                                                                        ggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGt    >  1:2551517/1‑150 (MQ=255)
                                                                                                                        ggTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGt    >  2:98716/1‑150 (MQ=255)
                                                                                                                         gTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTa   <  2:395627/150‑1 (MQ=255)
                                                                                                                         gTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTa   <  1:4136218/150‑1 (MQ=255)
                                                                                                                          tGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGACATCAGGTGCGGAAGTAc  >  1:4260158/1‑150 (MQ=255)
                                                                                                                          |                                                                                                                                                     
GCCGTGGCTGATGGCGACACTGGACACGGTTGCGCCAGGTTGCAGGCATTCCTGAACGACCTGGGCTTTGAACGGTTTCGGATAAGAGCTTCGTTGGCGCATGGAAATCCTGGCAATAAGGGTGATCGCGTCCGCTTAAAAATACGCGGACACCATCGCCCTTAATGCTGGAGTTCGGAAGGTGAGTTCCCCGGACGCTTACTCGCTGTATGCGTTCGAGACGGGAGCAGAGTTGTACTTTCTGGACGCAGGGAAATCAGGTGCGGAAGTAC  >  NC_002947/6157067‑6157338

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: