Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,554,884 T→A 5.9% intergenic (+577/+84) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,554,8840TA 5.9% 693.7 / 24.0 230intergenic (+577/+84)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (84/131);  new base A (11/2);  total (95/135)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 238 total aligned reads displayed.

TCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGT  >  NC_002947/5554735‑5555033
                                                                                                                                                     |                                                                                                                                                     
   gTAATCATTTTCGCGCTCTCATAAGTCACTTCGCATCAGTTGGTAAGTCTGGTGTACAGCCTGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGACAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                    <  2:1511244/150‑1 (MQ=18)
    tACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCa                                                                                                                                                   <  2:8982196/150‑1 (MQ=255)
     aCTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaa                                                                                                                                                  >  2:3897970/1‑150 (MQ=255)
      cTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                                 <  1:160623/150‑1 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTGTTAGCAAAg                                                                                                                                                >  1:281978/1‑150 (MQ=21)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                >  2:4680692/1‑150 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                >  2:8696750/1‑150 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                >  2:6725888/1‑150 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                >  2:4922267/1‑150 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGATTCGCAAAg                                                                                                                                                >  2:4922079/1‑150 (MQ=255)
       tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGATTAGCaaaa                                                                                                                                                >  1:7745207/1‑149 (MQ=255)
           tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCg                                                                                                                                            >  1:2038735/1‑150 (MQ=255)
             tCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGa                                                                                                                                          <  2:2465351/150‑1 (MQ=255)
             tCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGa                                                                                                                                          <  1:2107569/150‑1 (MQ=255)
                  ctctCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCCGCAAcgg                                                                                                                                     >  2:6116220/1‑148 (MQ=255)
                   tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCGGaaaaagg                                                                                                                                    >  1:8943671/1‑146 (MQ=255)
                   tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCCAAGCCCCGAAACAgg                                                                                                                                    >  1:4496876/1‑150 (MQ=255)
                   tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgg                                                                                                                                    <  2:1119547/150‑1 (MQ=255)
                    ctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                   <  1:2461092/150‑1 (MQ=255)
                     tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                  <  1:8244269/150‑1 (MQ=255)
                        tAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAATCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAATCCCCGAAACAGGGGttt                                                                                                                               >  1:5934555/1‑150 (MQ=17)
                        tAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCGGAAACTGGGGttt                                                                                                                               >  2:2545019/1‑150 (MQ=255)
                        tAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTAGCAAAGACAGGAAACAGGGGttt                                                                                                                               >  2:1852210/1‑150 (MQ=17)
                               cTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGATTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCATCtgg                                                                                                                        <  1:8597009/150‑1 (MQ=255)
                                tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                       >  2:4635255/1‑150 (MQ=255)
                                  cgcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACACGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggtgg                                                                                                                     <  1:5483785/150‑1 (MQ=255)
                                  cgcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggtgg                                                                                                                     <  1:6335282/150‑1 (MQ=255)
                                   gcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGACCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                    <  2:1661232/150‑1 (MQ=255)
                                   gcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTCTTGGGAAAGCCCCGAAAAAGGGGTTTAAGCagggggg                                                                                                                    >  1:3353801/1‑149 (MQ=17)
                                   gcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                    <  2:6261888/150‑1 (MQ=255)
                                      tCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctgc                                                                                                                 <  2:5642183/150‑1 (MQ=255)
                                        aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                               <  2:4150198/150‑1 (MQ=255)
                                        aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                               <  2:9061620/150‑1 (MQ=255)
                                        aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                               <  1:5767051/150‑1 (MQ=255)
                                         gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCGGAAACGGGGGTTTCAGCGGGTGGCGGCGTg                                                                                                              >  2:1190408/1‑150 (MQ=17)
                                          ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                             <  1:998466/150‑1 (MQ=255)
                                          ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                             <  2:7330097/150‑1 (MQ=255)
                                          ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                             <  2:4230977/150‑1 (MQ=255)
                                                cGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAACGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGttt                                                                                                       <  2:4973544/150‑1 (MQ=255)
                                                cGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGttt                                                                                                       <  1:4550011/150‑1 (MQ=255)
                                                 gTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGATTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                      <  1:82810/150‑1 (MQ=255)
                                                  tCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAg                                                                                                     <  1:4226181/150‑1 (MQ=255)
                                                  tCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCATCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCGTGAGTGTTTat                                                                                                     >  2:7999136/1‑149 (MQ=255)
                                                   cTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGa                                                                                                    <  1:8386406/150‑1 (MQ=255)
                                                         tACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCATGGCAAAGCCCGGAAACTGGGGTTTCAGCGGGTGGCTGCGTGAGTGTTTAGAAAACaa                                                                                              >  1:708903/1‑150 (MQ=17)
                                                          aCAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                             <  1:3873955/150‑1 (MQ=255)
                                                             gCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                          <  2:5958186/150‑1 (MQ=255)
                                                               cGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACGGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACCATGGGtt                                                                                        >  2:5134951/1‑150 (MQ=255)
                                                                ggTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTg                                                                                       <  1:4293594/150‑1 (MQ=255)
                                                                       cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCGTCGCCAAGCCCGGAAACTGGGGTTTGGGCGGGTGGGTGGTTGGGTGTTTTTTAAAAAATGGGTTTCagtgtg                                                                                >  2:3139696/1‑145 (MQ=17)
                                                                       cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                                >  2:3509837/1‑150 (MQ=255)
                                                                        aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                               <  2:3177035/150‑1 (MQ=255)
                                                                        aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                               >  2:4547369/1‑150 (MQ=255)
                                                                        aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTCGCAAAGCCCCGAAACAGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAAAAAGGGTTGCAATATGt                                                                               >  1:3751090/1‑150 (MQ=17)
                                                                         agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                              >  2:8295405/1‑150 (MQ=255)
                                                                          gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAACGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                             <  1:6844226/150‑1 (MQ=255)
                                                                          gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                             <  1:3672532/150‑1 (MQ=255)
                                                                               aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCa                                                                        <  2:7457955/150‑1 (MQ=255)
                                                                                 gCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTGGGTCGCAAAACCCCGAAACAGGGGTTTGAGCTGGTGGGTGCTTGTTTGTTTAGAAAAAAAAGTGTTTCAATCTGttttaacacg                                                                      >  1:8842605/1‑143 (MQ=17)
                                                                                    tATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCt                                                                   >  2:6836804/1‑150 (MQ=255)
                                                                                       gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAc                                                                <  1:421598/150‑1 (MQ=255)
                                                                                       gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAc                                                                <  1:300954/150‑1 (MQ=255)
                                                                                       gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTCGCAAAGCCCGGGAACGGGGGTTTCAGCTGGTGGCGGCGTGAGTGGTTAGGAAACAAAGGGTTGCATTGTGTTTGACCCGGGCTTAc                                                                >  2:6939207/1‑150 (MQ=17)
                                                                                         tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCGTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAACt                                                              >  2:5554597/1‑150 (MQ=255)
                                                                                         tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACATGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTTGAAAAAAAAGGGTTGCATTGTGTTTGAACAGGGCTTAcca                                                              >  2:2773410/1‑149 (MQ=17)
                                                                                            aCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                           <  1:3250371/150‑1 (MQ=255)
                                                                                            aCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                           <  1:761424/150‑1 (MQ=255)
                                                                                             cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                          <  1:6591130/150‑1 (MQ=255)
                                                                                                    ttGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTa                                                   <  2:3955603/150‑1 (MQ=255)
                                                                                                       aGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAATCCCCGAAACAGGGGTTTCAGCTGGTGGGTGCGTGAGTGTTTATAAAACAAAGGGTTGCATTCTGTTTGGCCAGGGCTTAACTCCAGTACTCTATcc                                                >  2:5369824/1‑150 (MQ=18)
                                                                                                               cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCCGAAAATGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAAAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTAAAGTACTCTATACAAAGaaaa                                        >  2:6480055/1‑150 (MQ=21)
                                                                                                                gAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaa                                       <  2:628499/150‑1 (MQ=255)
                                                                                                                 aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaa                                      <  2:8786013/150‑1 (MQ=255)
                                                                                                                 aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACGGGGGTTTCAGCTGGTGGCTGCTTGAGTGGTTAGAAAACAATGGGGTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaa                                      >  2:2948558/1‑150 (MQ=255)
                                                                                                                  aCCAGAAGCCAAGACCGTCTGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTTTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaa                                     <  2:819607/150‑1 (MQ=255)
                                                                                                                   ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                    <  2:8837144/150‑1 (MQ=255)
                                                                                                                      gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                                 <  1:1589277/150‑1 (MQ=255)
                                                                                                                      gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                                 <  2:3353801/150‑1 (MQ=255)
                                                                                                                            aaGACCGGCAGGTCCGGGCTTTGCTTCGCAAACCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGATTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCa                           <  2:4247213/150‑1 (MQ=255)
                                                                                                                                       ttCCGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGc                >  1:7907515/2‑150 (MQ=255)
                                                                                                                                        tCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCa               <  1:6140097/150‑1 (MQ=255)
                                                                                                                                         ccGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaa              >  1:6903857/1‑150 (MQ=255)
                                                                                                                                          cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaa             >  1:1318329/1‑150 (MQ=255)
                                                                                                                                          cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTAGAGTGTTTAGAAAACAATGGGTTGCATTCTGATTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaa             <  2:5529918/150‑1 (MQ=255)
                                                                                                                                          cGGGCTTTGCGAAGCAAAGCCCGGACACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGTCCAGGGCTTACCTCAAGTACTCTATCCAAAGAAATAAAAGTACTATCAGTGCTATCAGCaat             >  2:5183161/1‑149 (MQ=25)
                                                                                                                                           gggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa            >  2:2362979/1‑150 (MQ=255)
                                                                                                                                           gggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa            >  1:5045784/1‑150 (MQ=255)
                                                                                                                                           gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa            >  1:7379527/1‑150 (MQ=255)
                                                                                                                                           gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa            >  1:2760898/1‑150 (MQ=255)
                                                                                                                                            ggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCGATCAGCAAAAt           <  1:3143935/150‑1 (MQ=255)
                                                                                                                                            ggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAt           >  1:4716319/1‑150 (MQ=255)
                                                                                                                                              cTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt         <  1:5162692/150‑1 (MQ=255)
                                                                                                                                              cTTTGCTCCTCCAAGCCCGGAAACTGGGGTTTCCGCTGGGGGCTGCTTGAGTGTTTAGAAAACAAGGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt         <  2:6163219/150‑1 (MQ=255)
                                                                                                                                               tttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTg        <  2:700233/150‑1 (MQ=255)
                                                                                                                                                ttGCTTCGCAAAGCCCTGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAATAAAGAACAATCAGTGCAATCAGCAAAATTTGt       <  2:7745207/150‑1 (MQ=255)
                                                                                                                                                 tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa      >  2:7138952/1‑150 (MQ=255)
                                                                                                                                                 tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa      <  1:3509837/150‑1 (MQ=255)
                                                                                                                                                 tGCTCCGCAACGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACTATCAGTGCAATCAGCAAAATTTGTa      <  1:4680724/150‑1 (MQ=255)
                                                                                                                                                  gCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAc     >  2:2082312/1‑150 (MQ=255)
                                                                                                                                                    ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAccg   <  1:4922079/150‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGT  >  NC_002947/5554735‑5555033

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