Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,404,477 A→G 5.1% intergenic (+36/‑116) PP_3010 → / → PP_3011 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,404,4770AG 5.1% 500.5 / 9.9 155intergenic (+36/‑116)PP_3010/PP_3011hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (86/61);  new base G (3/5);  total (89/66)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 161 total aligned reads displayed.

CTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAA  >  NC_002947/3404328‑3404626
                                                                                                                                                     |                                                                                                                                                     
cTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCgaga                                                                                                                                                       <  2:3594406/150‑1 (MQ=255)
 ttGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGagag                                                                                                                                                      >  1:5033603/1‑150 (MQ=255)
    tAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAATAGCCCGCCGTCGAGAGCgg                                                                                                                                                   >  2:3473407/1‑150 (MQ=255)
         ttGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttt                                                                                                                                              <  1:4772377/150‑1 (MQ=255)
          tGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCtatt                                                                                                                                             <  1:5076688/150‑1 (MQ=255)
           gCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttttt                                                                                                                                            <  2:3078852/150‑1 (MQ=255)
             ggggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTCCGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgt                                                                                                                                          >  2:2825360/1‑150 (MQ=255)
                gCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgtg                                                                                                                                       >  1:3358291/1‑150 (MQ=255)
                gCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgtg                                                                                                                                       <  1:2351313/150‑1 (MQ=255)
                gCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgtg                                                                                                                                       <  1:4956896/150‑1 (MQ=38)
                  ccGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCAGGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAAACGCCCGCTCTCGACGGCGGGCTATTggtgtggg                                                                                                                                     <  2:1811924/149‑3 (MQ=38)
                    gggAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCa                                                                                                                                   >  2:823365/1‑150 (MQ=255)
                    gggAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCa                                                                                                                                   >  2:288747/1‑150 (MQ=255)
                    gggAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCa                                                                                                                                   >  2:4912307/1‑150 (MQ=255)
                      gAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATc                                                                                                                                 <  1:1134984/150‑1 (MQ=255)
                       aaGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgtgggggttg                                                                                                                                <  1:4569690/150‑8 (MQ=37)
                       aaGCGAGCGCCACTTCACGGCCTACCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCg                                                                                                                                <  2:3619259/150‑1 (MQ=255)
                         gcgagcgCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGcc                                                                                                                              >  2:1996440/1‑150 (MQ=255)
                               gcCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCg                                                                                                                        >  1:1565171/1‑150 (MQ=255)
                               gcCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCg                                                                                                                        >  1:1545344/1‑150 (MQ=255)
                                ccACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                                       >  1:842690/1‑150 (MQ=255)
                                ccACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                                       >  1:4454607/1‑150 (MQ=255)
                                 cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGc                                                                                                                      >  2:1287776/1‑150 (MQ=255)
                                 cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGc                                                                                                                      >  2:498928/1‑150 (MQ=255)
                                   cTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCaa                                                                                                                    >  1:363291/1‑150 (MQ=255)
                                       aCGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGtt                                                                                                                <  1:1479818/150‑1 (MQ=255)
                                         ggCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATCGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCa                                                                                                              <  2:4613234/150‑1 (MQ=255)
                                           ccTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaa                                                                                                            <  1:3866672/150‑1 (MQ=255)
                                            cTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaa                                                                                                           >  1:3013531/1‑150 (MQ=255)
                                              aTCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaa                                                                                                         >  2:1833191/1‑150 (MQ=255)
                                              aTCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaa                                                                                                         >  2:4257262/1‑150 (MQ=255)
                                               tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaaa                                                                                                        >  2:4213226/1‑150 (MQ=255)
                                                  aCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAggg                                                                                                     >  2:153347/1‑150 (MQ=255)
                                                     tGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCgg                                                                                                  >  2:2566758/1‑150 (MQ=255)
                                                      gATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCggg                                                                                                 <  1:3783456/150‑1 (MQ=255)
                                                        tGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTg                                                                                               <  2:1698499/150‑1 (MQ=255)
                                                         gCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTgg                                                                                              >  1:4651586/1‑150 (MQ=255)
                                                          ccTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGa                                                                                             <  1:913627/150‑1 (MQ=255)
                                                           cTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAg                                                                                            >  1:2866451/1‑150 (MQ=255)
                                                            tGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGc                                                                                           >  1:2106485/1‑150 (MQ=255)
                                                              gTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                                         >  2:3934528/1‑150 (MQ=255)
                                                              gTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                                         >  2:3934543/1‑150 (MQ=255)
                                                                     cATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                  >  2:4258251/1‑150 (MQ=255)
                                                                     cATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                  >  1:3642630/1‑150 (MQ=255)
                                                                     cATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTGTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                  >  1:4578855/1‑150 (MQ=255)
                                                                      ctGGCAGGTTGCTGGATGCTGCACTGAGGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCATTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGcc                                                                                 <  1:5043375/149‑1 (MQ=255)
                                                                                 cTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGCGCAGCGtt                                                                      <  2:615963/150‑1 (MQ=255)
                                                                                  tGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTg                                                                     >  1:125455/1‑150 (MQ=255)
                                                                                  tGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTg                                                                     >  1:2678675/1‑150 (MQ=255)
                                                                                      tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCACTAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCAGCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                 <  2:2209594/150‑1 (MQ=255)
                                                                                      tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                 <  2:1171470/150‑1 (MQ=255)
                                                                                      tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                 <  2:1636436/150‑1 (MQ=255)
                                                                                      tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                 <  1:4186840/150‑1 (MQ=255)
                                                                                      agctgcACAGATGTGTTCAACAGCGTAATCAGCCCACCCCCCCCCCCATCGCCCGCCGTCGAGAGCGGGCTATTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                 <  2:5023883/149‑1 (MQ=255)
                                                                                         tgcACTGATGTGTTCAAAATCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                              <  2:3997497/150‑1 (MQ=255)
                                                                                         tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAATGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                              >  1:1760349/1‑150 (MQ=255)
                                                                                         tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                              >  2:2225427/1‑150 (MQ=255)
                                                                                         tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                              >  1:769355/1‑150 (MQ=255)
                                                                                              tGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGc                                                         <  1:2533864/150‑1 (MQ=255)
                                                                                              tGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGc                                                         >  1:602912/1‑150 (MQ=255)
                                                                                                aTGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCg                                                       >  2:30688/1‑150 (MQ=255)
                                                                                                   tgtTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcc                                                    >  2:1001443/1‑150 (MQ=255)
                                                                                                     ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                  >  2:325807/1‑150 (MQ=255)
                                                                                                     ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                  >  1:868348/1‑150 (MQ=255)
                                                                                                     ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                  >  1:389139/1‑150 (MQ=255)
                                                                                                     ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                  >  1:4788176/1‑150 (MQ=255)
                                                                                                       cAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgc                                                <  2:2283358/150‑1 (MQ=255)
                                                                                                       cAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgc                                                >  2:1445212/1‑150 (MQ=255)
                                                                                                        aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgcc                                               >  2:4240150/1‑150 (MQ=255)
                                                                                                        aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgcc                                               >  1:4471872/1‑150 (MQ=255)
                                                                                                        aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgcc                                               >  2:1571167/1‑150 (MQ=255)
                                                                                                           aGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGt                                            >  1:798611/1‑150 (MQ=255)
                                                                                                              gTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAat                                         >  2:2656896/1‑150 (MQ=255)
                                                                                                                aaTCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAatat                                       <  1:823365/150‑1 (MQ=255)
                                                                                                                  tcatcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTc                                     <  1:2621066/150‑1 (MQ=255)
                                                                                                                    atcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGt                                   <  2:326173/150‑1 (MQ=255)
                                                                                                                    atcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGt                                   >  2:3420875/1‑150 (MQ=255)
                                                                                                                     tcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                                  <  2:602571/150‑1 (MQ=255)
                                                                                                                     tcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                                  <  2:4193621/150‑1 (MQ=255)
                                                                                                                        aaCCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTa                               >  1:2090529/1‑150 (MQ=255)
                                                                                                                           ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                            >  2:4756265/1‑150 (MQ=255)
                                                                                                                           ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                            >  1:4696500/1‑150 (MQ=255)
                                                                                                                           ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                            >  1:4696493/1‑150 (MQ=255)
                                                                                                                             ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtgg                          >  2:598108/1‑150 (MQ=255)
                                                                                                                             ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtgg                          >  2:1119676/1‑150 (MQ=255)
                                                                                                                               acacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggat                        <  1:1736725/150‑1 (MQ=255)
                                                                                                                                 acaaaaaGCCCGCTCTCGACGGCGGGCTATTGGTGTGGCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgg                      >  2:2599236/7‑150 (MQ=38)
                                                                                                                                  ccactcgCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAACAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga                     <  2:5033603/144‑1 (MQ=255)
                                                                                                                                   cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGgag                    >  2:755391/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCCGCTCTCGACGGCGGGCTATTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGgag                    >  1:1148956/5‑150 (MQ=255)
                                                                                                                                    aaTAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGgaga                   >  1:1097637/1‑150 (MQ=255)
                                                                                                                                     aTAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGagag                  >  1:3158854/1‑150 (MQ=255)
                                                                                                                                        gCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCg               >  1:558855/1‑150 (MQ=255)
                                                                                                                                          ccgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCa             <  2:3193938/150‑1 (MQ=255)
                                                                                                                                           cgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAg            <  2:968733/150‑1 (MQ=255)
                                                                                                                                           cgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAg            >  1:1023395/1‑150 (MQ=255)
                                                                                                                                             ccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGt          <  2:3670639/150‑1 (MQ=255)
                                                                                                                                               gTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTgg        <  1:1001443/150‑1 (MQ=255)
                                                                                                                                                   agagCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGc    <  1:1833191/150‑1 (MQ=255)
                                                                                                                                                     agCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCaa  <  2:4454607/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAA  >  NC_002947/3404328‑3404626

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: