Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,960,202 T→A 9.4% intergenic (+40/‑153) PP_1752 → / → PP_1753 hypothetical protein/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,960,2020TA 9.4% 917.6 / 58.5 292intergenic (+40/‑153)PP_1752/PP_1753hypothetical protein/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (124/139);  new base A (12/16);  total (136/156)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

ATTGCCCCTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGA  >  NC_002947/1960052‑1960350
                                                                                                                                                      |                                                                                                                                                    
    ccccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAcc                                                                                                                                                   >  2:2085454/1‑150 (MQ=255)
     cccTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTAAccc                                                                                                                                                  >  2:4546673/1‑150 (MQ=255)
       cTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCcgc                                                                                                                                                >  2:1649085/1‑150 (MQ=255)
        ttGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCcgcg                                                                                                                                               >  2:4018775/1‑150 (MQ=255)
          ggACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGaa                                                                                                                                             >  2:899122/1‑150 (MQ=255)
          ggACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  2:3097569/150‑1 (MQ=255)
          ggACGCTGGTGGAGGCCAATGACAAGCGCTGCGCGCGGGTGAAGGTGGTGCCCACAATCAACCAGGCGCTTGAGGCGGCGATTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  2:5080469/150‑1 (MQ=255)
           gACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGaaa                                                                                                                                            <  1:5042012/150‑1 (MQ=255)
            aCGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAg                                                                                                                                           <  1:681664/150‑1 (MQ=255)
            aCGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaagc                                                                                                                                           <  1:738904/150‑3 (MQ=255)
             cGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  1:4275646/150‑1 (MQ=255)
             cGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  1:4613625/150‑1 (MQ=255)
              gCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAggg                                                                                                                                         >  2:5485275/1‑150 (MQ=255)
               cTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGt                                                                                                                                        <  2:4492289/150‑1 (MQ=255)
                tggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTc                                                                                                                                       >  1:2182438/1‑150 (MQ=255)
                tggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTAACCGCGAAAGGGTc                                                                                                                                       >  1:4937072/1‑150 (MQ=255)
                 ggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCa                                                                                                                                      <  2:3215955/150‑1 (MQ=255)
                 ggtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCa                                                                                                                                      <  2:5030203/150‑1 (MQ=255)
                       ggCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGAAGCGGCCAGAACAg                                                                                                                                <  2:5196832/150‑1 (MQ=37)
                         ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGc                                                                                                                              <  2:2734826/150‑1 (MQ=255)
                         ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGc                                                                                                                              >  1:1868596/1‑150 (MQ=255)
                         ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGAAGCGGCCAGAACAGGc                                                                                                                              <  2:1605770/150‑1 (MQ=25)
                           aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCaa                                                                                                                            >  2:3146902/1‑150 (MQ=255)
                            aTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAac                                                                                                                           <  1:1316998/150‑1 (MQ=255)
                             tGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAaca                                                                                                                          >  2:2880982/1‑150 (MQ=255)
                                 aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTc                                                                                                                      >  1:1434784/1‑150 (MQ=255)
                                   gcgcTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCa                                                                                                                    <  1:5474070/150‑1 (MQ=255)
                                      ctcggCGCGCGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGCCGGCGTGTGCCAACCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTACGCGGGTTAACCCGCGAAAGGCACAGCACAAGCAACACCTCCAccc                                                                                                                 >  2:652552/4‑150 (MQ=255)
                                       tGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCa                                                                                                                <  1:5054394/150‑1 (MQ=255)
                                            gcgGGTGAATGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCAt                                                                                                           <  2:5583613/150‑1 (MQ=255)
                                              gggTGAAGGTGGTGCGCACAATCAACCCGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGACGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGc                                                                                                         <  1:1105689/150‑1 (MQ=255)
                                                 tGAAGGAGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcccc                                                                                                      >  2:1160029/1‑150 (MQ=255)
                                                    aGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaa                                                                                                   <  1:3306409/150‑1 (MQ=255)
                                                     ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                  <  1:1055106/150‑1 (MQ=255)
                                                     ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTAAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                  <  2:39541/150‑1 (MQ=255)
                                                        ggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                               <  1:4151987/150‑1 (MQ=255)
                                                        ggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                               <  1:3835580/150‑1 (MQ=255)
                                                         gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCGCCACCAAGCCAAGCCCCAAAAGAAt                                                                                              >  1:955432/1‑150 (MQ=255)
                                                           gcgcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGa                                                                                            <  2:4144537/150‑1 (MQ=255)
                                                              cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCt                                                                                         <  1:4676091/150‑1 (MQ=255)
                                                              cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCt                                                                                         >  1:25045/1‑150 (MQ=255)
                                                                caATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtga                                                                                       <  2:279684/150‑1 (MQ=255)
                                                                 aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACCCCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                      >  2:4953359/1‑150 (MQ=255)
                                                                 aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                      >  2:2342071/1‑150 (MQ=255)
                                                                   tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgatga                                                                                    >  1:5143372/1‑150 (MQ=255)
                                                                   tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCAAGCCCCAAAAGAATGACCtgatga                                                                                    >  1:1620691/1‑150 (MQ=255)
                                                                   tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCCCGAAAGGGTCAGCACAGGCAACACCACCACCCACCCATGCCCCAAAAGAATGACCtgatga                                                                                    >  1:192241/1‑150 (MQ=255)
                                                                      aCCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTAGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAAtt                                                                                 >  1:1744914/1‑150 (MQ=255)
                                                                         aGGCGCTTGTGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGATTCGCGGGTTTACCCGCGAAAGGGTAAGAACTGGAAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                              <  1:3047487/150‑1 (MQ=255)
                                                                         aGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCtt                                                                              >  2:3445516/1‑150 (MQ=255)
                                                                              cTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTtctc                                                                         <  2:4394719/150‑1 (MQ=255)
                                                                               ttGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGTGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCg                                                                        <  1:2611107/150‑1 (MQ=255)
                                                                                  aGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCa                                                                     <  1:4570642/150‑1 (MQ=255)
                                                                                  aGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCa                                                                     <  1:4025671/150‑1 (MQ=255)
                                                                                      gtcgTTTGCCAAGCACATTAAATAGCCTATGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAATGGTTCAGCACATGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAACTGATGAATTCTTTTCTCGGCActct                                                                 <  1:5340964/148‑1 (MQ=255)
                                                                                       gcgTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACtctc                                                                <  2:1340120/150‑1 (MQ=255)
                                                                                         gTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATTCCCCCAAAGAATGAACTGATGAATTCTTTTCTCGGCACTCTCCa                                                              >  1:4920247/1‑150 (MQ=255)
                                                                                              ccAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtc                                                         >  2:4266155/1‑150 (MQ=255)
                                                                                              ccAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtc                                                         <  2:4876087/150‑1 (MQ=255)
                                                                                               cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtct                                                        >  2:4182586/1‑150 (MQ=255)
                                                                                                aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGcctc                                                       >  2:3878859/1‑150 (MQ=255)
                                                                                                  gCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                      >  2:3562258/1‑149 (MQ=255)
                                                                                                  gCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                      <  1:3562257/149‑1 (MQ=255)
                                                                                                    acaAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTAGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCGGtc                                                   >  1:358026/1‑150 (MQ=255)
                                                                                                     caAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAACGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcc                                                  <  2:2394816/150‑1 (MQ=255)
                                                                                                      aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcct                                                 >  2:4633185/1‑150 (MQ=255)
                                                                                                      aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcct                                                 >  1:5022683/1‑150 (MQ=255)
                                                                                                      aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAACAAAGACCGTGTGGATTATTTTATAGGCACTATCCACCGGCGCTGtcct                                                 >  1:2922024/1‑150 (MQ=255)
                                                                                                       aGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCCGCACAGGCCACACCGCCCCCCAGCCCAGCCCCAAAAGAATGACCTGCTGAATTCTTTTCTCGGCACTCTCCACCGTCTCGGtccgc                                                >  1:370285/1‑148 (MQ=255)
                                                                                                       aGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcctc                                                >  2:4512291/1‑150 (MQ=255)
                                                                                                           aTAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAg                                                                                                                                <  1:4625047/66‑1 (MQ=255)
                                                                                                            tAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcc                                           >  2:772899/1‑150 (MQ=255)
                                                                                                             aGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCa                                                                     <  2:5520220/123‑1 (MQ=255)
                                                                                                              gCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAcccc                                         <  1:2876160/150‑1 (MQ=255)
                                                                                                               ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                        <  1:5222609/150‑1 (MQ=255)
                                                                                                               ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                        <  1:3117622/150‑1 (MQ=255)
                                                                                                                 gAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                      <  1:2715847/150‑1 (MQ=255)
                                                                                                                      tGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAAtt                                 <  1:4182586/150‑1 (MQ=255)
                                                                                                                             gCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCaacaac                          <  2:4789263/150‑1 (MQ=255)
                                                                                                                               cTGTGCTGACCCTTTCGCGGGTAAACCCGCGAAGCGGCCAGAACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                      >  2:3011081/1‑136 (MQ=25)
                                                                                                                                tGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcccc                       >  2:5153651/1‑150 (MQ=255)
                                                                                                                                tGTTCTGGCCGCTTCGCGGGGTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCAGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCCTAGAAATCAAAAAcccc                       >  1:601755/1‑150 (MQ=255)
                                                                                                                                       gCCGCTTCGCGGGTTTCCCCGCGAAAGGGTCAGCACAGGCCACACCTCCCCCCAGCCATGCCCCAAAAGCATGACCGGATGCATTCGTTTCGCGGCCCGCTCCACCGGCTCGGTCCGCCCACCCCGGAGACTTCAACCCCCCCCCCCCCg                >  2:4340686/1‑150 (MQ=255)
                                                                                                                                          gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGc             <  1:182938/150‑1 (MQ=255)
                                                                                                                                           tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCt            >  1:5585079/2‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGACTTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCtt           >  2:2639931/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTg          <  2:1439839/150‑1 (MQ=255)
                                                                                                                                             tCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTg          >  2:2305522/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTg          >  2:1251128/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGAt        >  1:1459562/1‑150 (MQ=255)
                                                                                                                                                cgGGTTTACCCGCGAATGGGTCAGCACAGGCAACACCTACACCCAGCCATGCCCCAAAAGAATGAGATGATGAATTCTTTTCTCGGGACTGTCCAGAGTCTCTGTGCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATc       >  1:2092852/1‑150 (MQ=255)
                                                                                                                                                cgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATc       <  2:1876197/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa     <  2:2183237/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa     <  1:1649085/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa     <  2:1839643/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa     <  1:1267862/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa     <  2:1366672/150‑1 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAg    <  1:2345415/150‑1 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAg    <  1:4526040/150‑1 (MQ=255)
                                                                                                                                                   gTTGACCCGCGCAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAACAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCCCCGTCTCTGTCCTCCACCCCCGTAGAATTCAACAACCCCCCCCCCGGGCTTGATCGAg    <  2:370285/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATTGCCCCTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGA  >  NC_002947/1960052‑1960350

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: