Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,475 G→T 5.9% intergenic (‑229/‑451) PP_5406 ← / → PP_5407 transposase/TnsD‑like transposition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,171,4750GT 5.9% 742.9 / 26.1 221intergenic (‑229/‑451)PP_5406/PP_5407transposase/TnsD‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base G (109/99);  new base T (5/8);  total (114/107)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 223 total aligned reads displayed.

CTGCGGGCAGGTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCTGACGC  >  NC_002947/6171327‑6171624
                                                                                                                                                    |                                                                                                                                                     
cTGCGGGCAGGTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGt                                                                                                                                                      >  1:980951/1‑150 (MQ=255)
 tGCGGGCAGGTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTc                                                                                                                                                     >  2:539620/1‑150 (MQ=255)
      gCAGGTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACt                                                                                                                                                >  1:3822250/1‑150 (MQ=255)
         ggTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTgg                                                                                                                                             <  1:1850340/150‑1 (MQ=255)
         ggTAACCGCGTTGCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTgg                                                                                                                                             <  2:4274074/150‑1 (MQ=255)
          gTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTggg                                                                                                                                            >  1:1347631/1‑150 (MQ=255)
             aCCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGTGCGACCTTGGGccc                                                                                                                                         <  2:1896497/150‑2 (MQ=255)
             aCCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCt                                                                                                                                         <  2:2232636/150‑1 (MQ=255)
              ccGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCtt                                                                                                                                        <  1:4127513/150‑1 (MQ=255)
               cgcgTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCttt                                                                                                                                       <  2:966607/150‑1 (MQ=255)
               cgcgTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCttt                                                                                                                                       <  1:2963188/150‑1 (MQ=255)
               cgcgTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAAGGCCCAAGTGCGACCTTGGGCCctt                                                                                                                                       <  2:5107281/150‑1 (MQ=255)
                 cgTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAAGGCCCAAGTGCGACCTTGGGCCCTTGt                                                                                                                                     <  1:147149/150‑1 (MQ=255)
                 cgTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAAGGCCCAAGTGCGACCTTGGGCCCTTGt                                                                                                                                     <  2:4633772/150‑1 (MQ=255)
                  gTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTc                                                                                                                                    <  2:2769698/150‑1 (MQ=255)
                  gTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTc                                                                                                                                    >  2:5461878/1‑150 (MQ=255)
                         gTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGa                                                                                                                             <  2:1371967/150‑1 (MQ=255)
                              ggATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTcc                                                                                                                        >  1:968212/1‑150 (MQ=255)
                               gATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCa                                                                                                                       <  2:2497928/150‑1 (MQ=255)
                                  aaCGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAgcg                                                                                                                    >  2:618580/1‑150 (MQ=255)
                                     gcgGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGctc                                                                                                                 >  2:857834/1‑150 (MQ=255)
                                      cgGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGctct                                                                                                                >  2:5190637/1‑150 (MQ=255)
                                       ggAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCtctc                                                                                                               >  2:1257385/1‑150 (MQ=255)
                                           aTTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGttt                                                                                                           <  2:5322652/150‑1 (MQ=255)
                                              cTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTtagc                                                                                                        <  1:4578300/150‑1 (MQ=255)
                                               ttGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTtagct                                                                                                       >  2:5574707/1‑150 (MQ=255)
                                               ttGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTtagct                                                                                                       >  2:4380319/1‑150 (MQ=255)
                                                tgtgGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTgg                                                                                                                                             >  1:2777085/1‑111 (MQ=255)
                                                tgtgGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTgg                                                                                                                                             <  2:2777085/111‑1 (MQ=255)
                                                tgtgGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTtagcta                                                                                                      >  1:2341888/1‑150 (MQ=255)
                                                  tgGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTtagctagc                                                                                                    <  2:2914888/150‑1 (MQ=255)
                                                      tGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCCCTCGTTTAGCTAGCTgtc                                                                                                >  2:2703431/1‑150 (MQ=255)
                                                       gAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTgtcg                                                                                               >  2:2071527/1‑150 (MQ=255)
                                                         gCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTgtcgtc                                                                                             >  2:303168/1‑150 (MQ=255)
                                                          cGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCt                                                                                            <  1:1232502/150‑1 (MQ=255)
                                                           gAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTg                                                                                           >  2:5296932/1‑150 (MQ=255)
                                                                cGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAAttt                                                                                      >  1:5548718/1‑150 (MQ=255)
                                                                    cGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAAtttcttt                                                                                  <  2:1175008/150‑1 (MQ=255)
                                                                    cGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAAtttcttt                                                                                  <  2:3039944/150‑1 (MQ=255)
                                                                      aaaCTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCg                                                                                >  1:1431874/1‑150 (MQ=255)
                                                                      aaaCTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCg                                                                                >  1:4376909/1‑150 (MQ=255)
                                                                       ggcatGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGt                                                                               <  2:5150658/146‑1 (MQ=255)
                                                                       aaCTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGt                                                                               >  1:775494/1‑150 (MQ=255)
                                                                        aCTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGtt                                                                              <  2:2369797/150‑1 (MQ=255)
                                                                        aCTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCCCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGtt                                                                              <  1:2703431/150‑1 (MQ=255)
                                                                         cTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTa                                                                             <  1:4209391/150‑1 (MQ=255)
                                                                          ttGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAg                                                                            >  1:3358727/1‑150 (MQ=255)
                                                                          ttGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAg                                                                            >  2:3787847/1‑150 (MQ=255)
                                                                               gTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACCTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTa                                                                       >  1:4785602/1‑150 (MQ=255)
                                                                                tttCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTAc                                                                      >  1:3896327/1‑150 (MQ=255)
                                                                                       tGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGc                                                               <  1:612251/150‑1 (MQ=255)
                                                                                          aaGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGtc                                                            <  2:3884384/150‑1 (MQ=255)
                                                                                          aaGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCAAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGtc                                                            >  1:4684087/1‑150 (MQ=255)
                                                                                          aaGGGATAGACAATGGAGGAGCCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGtc                                                            >  2:3611814/1‑150 (MQ=255)
                                                                                             ggATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACACGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTcc                                                         <  1:1907549/150‑1 (MQ=255)
                                                                                               aTAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTccgc                                                       <  1:3599891/150‑1 (MQ=255)
                                                                                                tAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTccgcc                                                      <  2:2661228/150‑1 (MQ=255)
                                                                                                  gACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGa                                                    >  1:3457276/1‑150 (MQ=255)
                                                                                                   aCAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGAt                                                   >  1:2455802/1‑150 (MQ=255)
                                                                                                    cAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGAtt                                                  <  2:1310735/150‑1 (MQ=255)
                                                                                                    cAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGAtt                                                  <  1:1034450/150‑1 (MQ=255)
                                                                                                       tGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATtgtg                                               >  1:1326212/1‑150 (MQ=255)
                                                                                                        ggaggaGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGa                                              <  1:1196423/150‑1 (MQ=255)
                                                                                                          aggagTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATc                                            >  1:1645680/1‑150 (MQ=255)
                                                                                                          aggagTCAGGTGCCTTTTCGGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATc                                            <  2:3996229/150‑1 (MQ=255)
                                                                                                              gTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCACGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGc                                        <  1:326652/150‑1 (MQ=255)
                                                                                                               tCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGcc                                       <  2:3284608/150‑1 (MQ=255)
                                                                                                                 aGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAgg                                                                           <  1:316861/112‑1 (MQ=255)
                                                                                                                 aGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAgg                                                                           >  2:316861/1‑112 (MQ=255)
                                                                                                                  ggTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCAc                                    >  1:505552/1‑150 (MQ=255)
                                                                                                                      ccTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTc                                >  2:4936251/1‑150 (MQ=255)
                                                                                                                        ttttCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCAc                              >  1:158428/1‑150 (MQ=255)
                                                                                                                            cTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGc                          >  2:761649/1‑150 (MQ=255)
                                                                                                                             tGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGcc                         <  2:5605054/150‑1 (MQ=255)
                                                                                                                              gACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTCTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGccc                        >  2:3511809/1‑150 (MQ=255)
                                                                                                                               ccgcagccACACGGGCCCACGGTCGCCCTTGGGCCTTTGTCGATCTGCTTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGACTTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCt                       <  2:4841706/144‑1 (MQ=255)
                                                                                                                                 ttATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTtc                     >  1:5125101/1‑150 (MQ=255)
                                                                                                                                 ttATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTtc                     >  1:2494671/1‑150 (MQ=255)
                                                                                                                                 ttATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTtc                     >  2:18042/1‑150 (MQ=255)
                                                                                                                                 ttATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTtc                     >  1:4310814/1‑150 (MQ=255)
                                                                                                                                    tCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTtctct                  <  1:3033523/150‑1 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  1:4917066/150‑1 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  2:4545677/150‑1 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  1:5420385/150‑1 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  1:67615/150‑1 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                >  1:3748302/1‑150 (MQ=255)
                                                                                                                                      cacaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  1:59197/150‑1 (MQ=255)
                                                                                                                                       acaAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGTCCACGCTCACGAGCCCTTCTCTCAg               >  2:2012690/1‑150 (MQ=255)
                                                                                                                                         aaaggCCCAAGTGCGACCTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCa                <  1:3889162/144‑1 (MQ=255)
                                                                                                                                         aaaggCCCAAGTGCGACCTTGGGCCCTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGc             >  1:1998514/4‑150 (MQ=255)
                                                                                                                                         aaaggCCCAAGTGCGACCTTGGGCCCTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGc             >  1:2175288/4‑150 (MQ=255)
                                                                                                                                         aaaggCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGGCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGc             >  2:3058030/4‑150 (MQ=255)
                                                                                                                                           gggCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTGTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCGGATTGTGATCTTGCCCACGCTCACGTGCCCTTCTCTCAGGCtt           >  2:1679822/1‑150 (MQ=255)
                                                                                                                                           aggCCCAAGTGCGACCTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCtt           >  1:4439129/2‑150 (MQ=255)
                                                                                                                                           aggCCCAAGTGCGACCTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCtt           >  2:3570606/2‑150 (MQ=255)
                                                                                                                                               ccAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCt       <  2:1518199/150‑1 (MQ=255)
                                                                                                                                               ccAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGAGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCt       <  1:1370615/150‑1 (MQ=255)
                                                                                                                                                 aaGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCTGa     <  1:618580/150‑1 (MQ=255)
                                                                                                                                                   ggTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCTGACg   <  1:1383655/150‑1 (MQ=255)
                                                                                                                                                    gTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCTGACGc  >  1:822705/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCGGGCAGGTAACCGCGTTTCCGGTTGTGGATAACGCGGAGATTCTTGTGGATGAGCGAATCCGTGCGAAACTTGACGTTTCACATGCAAGGGATAGACAATGGAGGAGTCAGGTGCCTTTTCTGACTTATCCACAAGGGCCCAAGGTCGCACTTGGGCCTTTGTCTATCTGATTTCCAGCGCTCTCGTTTAGCTAGCTGTCGTCTGAATTTCTTTCGTTAGGTCTACCTTAGGCGTCTCCGCCGATTGTGATCTTGCCCACGCTCACGAGCCCTTCTCTCAGGCTTCGCTGACGC  >  NC_002947/6171327‑6171624

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