Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,753,703 A→G 13.4% intergenic (‑163/‑308) thiI ← / → glnA tRNA sulfurtransferase/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,753,7030AG13.4% 611.4 / 71.3 249intergenic (‑163/‑308)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (99/107);  new base G (14/19);  total (114/135)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-02

Only 100 of 251 total aligned reads displayed.

GCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCA  >  NC_002947/5753554‑5753852
                                                                                                                                                     |                                                                                                                                                     
  gcAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAgg                                                                                                                                                     >  2:7434926/1‑150 (MQ=255)
  gcAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAgg                                                                                                                                                     <  1:2664763/150‑1 (MQ=255)
      ggCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGcttt                                                                                                                                                 <  2:8155443/150‑4 (MQ=39)
      ggCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGt                                                                                                                                                 <  2:1517322/150‑1 (MQ=255)
       gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACGGGctttg                                                                                                                                                <  2:4451571/150‑5 (MQ=255)
       gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTg                                                                                                                                                <  1:7874066/150‑1 (MQ=255)
       gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTg                                                                                                                                                <  2:4125128/150‑1 (MQ=255)
       gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCGACAGGCGGTg                                                                                                                                                >  1:4173863/1‑150 (MQ=255)
        ccTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGc                                                                                                                                               <  1:2412526/150‑1 (MQ=255)
         cTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCa                                                                                                                                              <  2:2895430/150‑1 (MQ=39)
           gCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCaaa                                                                                                                                            <  2:4080681/150‑1 (MQ=255)
           gCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCaaa                                                                                                                                            <  2:8718846/150‑1 (MQ=255)
             cGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGGGCAAAGc                                                                                                                                          >  1:1651203/1‑150 (MQ=255)
               aaaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                                        <  1:8545220/150‑1 (MQ=255)
               aaaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCATTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                                        <  1:6579510/150‑1 (MQ=255)
                 aaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAAGGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGt                                                                                                                                      <  2:1561481/150‑1 (MQ=255)
                  aTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTa                                                                                                                                     <  1:3694531/150‑1 (MQ=255)
                    gAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACGGGCTTTGCACCGCCTGTAgg                                                                                                                                   <  2:7910067/150‑1 (MQ=38)
                       gggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGc                                                                                                                                >  2:6667112/1‑150 (MQ=255)
                       gggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGc                                                                                                                                >  1:5543371/1‑150 (MQ=255)
                        ggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGcgcg                                                                                                                               <  2:397008/150‑1 (MQ=25)
                         gcgcGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                              <  1:6093353/150‑1 (MQ=38)
                         gcgcGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                              <  1:270312/150‑1 (MQ=255)
                          cgcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                             <  1:6986028/150‑1 (MQ=255)
                          cgcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                             <  1:1721257/150‑1 (MQ=255)
                           gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGt                                                                                                                            <  1:644255/150‑1 (MQ=255)
                               aTTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAc                                                                                                                        <  2:7004241/150‑1 (MQ=255)
                                ttATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAcc                                                                                                                       <  1:4414641/150‑1 (MQ=255)
                                 tataTCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAccc                                                                                                                      <  1:5499510/150‑1 (MQ=255)
                                  atatCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcg                                                                                                                     <  1:5263448/150‑1 (MQ=255)
                                     tCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                              >  2:8584429/1‑138 (MQ=37)
                                     tCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:190167/1‑150 (MQ=255)
                                     tCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:1133449/1‑150 (MQ=255)
                                       ggAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGcc                                                                                                                                         >  2:7283607/1‑125 (MQ=255)
                                       ggAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGcc                                                                                                                                         <  1:7283607/125‑1 (MQ=255)
                                        gAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                               <  1:8628262/150‑1 (MQ=255)
                                               cTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGAGAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTa                                                                                                        <  2:8444643/150‑1 (MQ=255)
                                                tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGGAGGAGCGGGTTTACCCGCGAAGAGGCCCTAg                                                                                                       >  1:4855267/1‑150 (MQ=255)
                                                tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCGACAGGCGGTGCAAAGCCCGGAGGAGCGGGTTTAACCGCGAAGAGGCCCAAg                                                                                                       >  1:8599262/1‑150 (MQ=255)
                                                  aGGATTTGACCAACTTTTTATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGTTAAACCCGCGCCTACAGGGGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGaa                                                                                                     <  1:7981713/150‑1 (MQ=255)
                                                     aTTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCt                                                                                                  >  1:1399535/1‑150 (MQ=255)
                                                      tttGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTg                                                                                                 <  2:6186702/150‑1 (MQ=255)
                                                      tttGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTg                                                                                                 <  2:2192102/150‑1 (MQ=255)
                                                      tttGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTg                                                                                                 <  1:7478713/150‑1 (MQ=255)
                                                       ttGCCCAACTTTTGATCAGACGATTTGATATAAATGCAGAGGCCAGTTTTTGCAGCGCCTGTTCGGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                                <  1:1232684/150‑1 (MQ=37)
                                                       ttGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                                <  1:2240949/150‑1 (MQ=255)
                                                          aCCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCGCCAACAGGCGGTGCAAAGCCCGTAGGAGCGGGGTTACCCCCGAAGAGGCCCTAGAACCAGGcac                                                                                             >  1:5753837/1‑150 (MQ=38)
                                                            cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCacac                                                                                           <  2:3621620/150‑1 (MQ=255)
                                                            cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCacac                                                                                           <  1:8345442/150‑1 (MQ=255)
                                                            cAACTCCTGATCAGAAGCGTTGAGAGAAATGCCGAGGCCAGGATTTGCAGCGCCTGGTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCacac                                                                                           <  2:3621548/150‑1 (MQ=37)
                                                             aaCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCAcacc                                                                                          >  2:5755889/1‑150 (MQ=255)
                                                              aCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAATGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCAcacca                                                                                         <  1:2174947/150‑1 (MQ=255)
                                                               cTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCAGGCAcaccac                                                                                        >  2:3082311/1‑150 (MQ=255)
                                                                ttttGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTACAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCAGGCAcaccaca                                                                                       >  1:8831310/1‑149 (MQ=255)
                                                                   tGATCATAAGCTTTGATATAAATGCAGTGGCCTGGTTTTGCAGCGCGTGTTCGGGCCCTTTCGGGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGTGCGGTTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcac                                                                                    <  1:1176329/150‑1 (MQ=255)
                                                                   tGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcac                                                                                    <  2:5230195/150‑1 (MQ=255)
                                                                   tGAACAGACGCTTTGAGATAAATGCAGCGGCCAGGTTTGGCAGCGCCTGTTCCGGCCCTTCCGCGGTTAACCCCGCGCCTCCCGGCGGTGCAAAGCCCGTCGGCGGGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcac                                                                                    <  1:2521156/150‑1 (MQ=37)
                                                                     aTCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgc                                                                                  <  2:542648/150‑1 (MQ=255)
                                                                      tCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                              <  2:5238707/105‑1 (MQ=21)
                                                                      tCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                              >  1:5238707/1‑105 (MQ=21)
                                                                        aGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgccac                                                                               >  2:2094048/1‑150 (MQ=255)
                                                                        aGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTAAAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTTACCCCGAAGAGGCCCTAGAACCAGGCACACCACGcaagccac                                                                               >  1:6364298/1‑150 (MQ=255)
                                                                         gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTctcgcc                                                                                 <  2:8953189/147‑1 (MQ=255)
                                                                         gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgcc                                                                                 >  1:8953189/1‑147 (MQ=255)
                                                                             cTTTGACATAAATGCAGAGGCCAGTTTTTGCAGCGCCTGTTCCGGCCCTTACGCGGGTAAACCCGCGCCTACACGCGGTGCAAAGCCCGTAGCAGCGGGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCAcc                                                                          <  1:2052167/150‑1 (MQ=255)
                                                                              tttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAAGCTGGCACACCACTCACGCCACGCAccc                                                                         >  2:8791791/1‑150 (MQ=255)
                                                                                  agaTAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTAtt                                                                     <  2:670886/150‑1 (MQ=255)
                                                                                  acaTAAATGCAGAGGCAAGGTTTTGCAGCGCCTGTACCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTAtt                                                                     <  2:2748293/148‑1 (MQ=255)
                                                                                   gaTAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTg                                                                    <  2:3711298/150‑1 (MQ=255)
                                                                                    aTAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCCCGCCTACACGCGGTGCAAAGCCCGTAGGAGCGGGTTTAACCGCGCAGAGGCCCTAGAAACTGGCACAACAAAAACCCCACCCACCCAATTgg                                                                   >  1:3347064/1‑150 (MQ=255)
                                                                                            agagGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGacac                                                           <  2:7625501/150‑1 (MQ=255)
                                                                                              agGCCAGGTGTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGGAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCCCGCCCCCGATTGGTGCGACACGa                                                         >  2:6873931/1‑150 (MQ=38)
                                                                                                 ccAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCCCGCCTACAGGCGGTGCAAAGCCCGGAGGAGCGGGTTTTCCCGCGAAGAGGCCCTAGAAACAGGCACACCACTCACCCCACCCACCCCATTGGGGCGaccccaaca                                                      >  1:6162322/1‑143 (MQ=255)
                                                                                                  cAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTCGGAGCGGGGTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCCCGCACCCGATTGGGGCGACACGATCaa                                                     >  2:2265962/1‑150 (MQ=39)
                                                                                                  cAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCaa                                                     >  2:4964568/1‑150 (MQ=255)
                                                                                                      ttttGCAGCGCGTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCaaaaaa                                                 <  1:2007306/150‑1 (MQ=255)
                                                                                                          gCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCa                                             <  1:5677282/150‑1 (MQ=255)
                                                                                                          gCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCa                                             <  2:6332862/150‑1 (MQ=255)
                                                                                                          gCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCCCCGACAGGCGGTGCAAAGCCCGTAGGGGCCGGGTTACCCGCGAAGAGGCCCGAGAACCAGGCACACCCCACACGCCACGCACCCAATTGGTGCGACACGATCaaaaaaacca                                             >  1:6356177/1‑146 (MQ=255)
                                                                                                            aGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATc                                           <  2:7743103/150‑1 (MQ=255)
                                                                                                             gcgcCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCGACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCAGGCACACCACTCACGCCACGCACCCGATTGGTGCGACACGATCAAAAACGCATCGt                                         >  2:2406015/1‑150 (MQ=255)
                                                                                                              cgcCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGt                                         >  1:691674/1‑150 (MQ=255)
                                                                                                                ccTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCATACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTaa                                       >  2:4355953/1‑150 (MQ=255)
                                                                                                                    ttCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCGTCGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGt                                   <  1:994871/150‑1 (MQ=255)
                                                                                                                    ttCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCGATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGt                                   >  1:3791323/1‑150 (MQ=255)
                                                                                                                        ggCCCTTTCGCGGGTAAAACCGCGCCTACAGGCGGTGCAAAGCCCCTAGGAGCGGGTTTAACCGCGAAGAGGCCCTAGAACCAGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATAGTAAGGGTGCAc                               >  1:5117764/1‑150 (MQ=39)
                                                                                                                                   gggTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCt                    <  2:4717822/150‑1 (MQ=255)
                                                                                                                                    ggTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTa                   >  2:7499968/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCGCCTACAGGGGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTaa                  <  1:3211763/150‑1 (MQ=255)
                                                                                                                                       aaaCCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAcc                <  2:5892062/150‑1 (MQ=255)
                                                                                                                                          cccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa             <  2:5543371/150‑1 (MQ=255)
                                                                                                                                          cccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa             <  1:517855/150‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACTTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc            >  2:3210099/1‑150 (MQ=37)
                                                                                                                                           ccGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc            >  1:1306155/1‑150 (MQ=37)
                                                                                                                                           ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc            >  1:3073170/1‑150 (MQ=255)
                                                                                                                                           ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc            >  1:5485989/1‑150 (MQ=255)
                                                                                                                                           ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc            >  1:565254/1‑150 (MQ=255)
                                                                                                                                                 cTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCaaa      >  1:2963842/1‑150 (MQ=255)
                                                                                                                                                  tACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAAt     >  2:7174192/1‑150 (MQ=255)
                                                                                                                                                  tACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAAt     >  1:4831707/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCA  >  NC_002947/5753554‑5753852

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: