Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,703 A→T 16.3% intergenic (+54/+48) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,7030AT16.3% 993.3 / 65.9 391intergenic (+54/+48)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (231/77);  new base T (2/68);  total (233/158)
Fisher's exact test for biased strand distribution p-value = 2.52e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 409 total aligned reads displayed.

TTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGG  >  NC_000913/3608554‑3608851
                                                                                                                                                     |                                                                                                                                                    
  aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                    <  2:4344997/150‑1 (MQ=255)
  aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                    <  2:6413104/150‑1 (MQ=255)
  aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTCTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                                    <  2:703617/150‑1 (MQ=255)
     gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                 <  1:8118907/150‑1 (MQ=255)
     gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                 >  1:8220282/1‑150 (MQ=255)
     gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                 <  1:8638461/150‑1 (MQ=255)
      aTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                <  1:4556632/150‑1 (MQ=255)
        gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                              <  2:1865106/150‑1 (MQ=255)
         aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGGTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                             >  2:5757638/1‑150 (MQ=255)
           cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                           >  1:3480165/1‑150 (MQ=255)
            gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                          >  1:7421410/1‑150 (MQ=255)
              gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                        <  1:6004308/150‑1 (MQ=255)
                 gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                                     >  2:2761810/1‑150 (MQ=255)
                    gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                  >  1:6680561/1‑150 (MQ=255)
                      tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                >  2:7752275/1‑150 (MQ=255)
                        gctgTGCTGGCTGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTATGATTGTTGCGCAGGTTATAAGGAAAACCTATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGTAAAAggcg                                                                                                                              <  2:1595147/147‑1 (MQ=255)
                        gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                              >  1:4866352/1‑150 (MQ=255)
                        gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                              >  1:587765/1‑150 (MQ=255)
                          aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                            >  2:5968149/1‑150 (MQ=255)
                                ggCTGATACGGGGGCGAGGGTGCTGGTGACAGCGATTGCGTTAAGTTTGTTGCGCAGGAGATACGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCGCTGGGGTGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                                      <  2:8196489/150‑1 (MQ=255)
                                ggCAGATACGGGGGCGCCGGTGCTGTTGACAGCGAATGCGTTAAGTTTGTTGCGCAGGAGTTAAGGCAAACCGATCGCACCATTGCGCGCGATGGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGaa                                                                                                                      <  2:1182457/150‑1 (MQ=38)
                                 gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                     >  1:2968296/1‑150 (MQ=255)
                                 gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                                     >  2:7132119/1‑150 (MQ=255)
                                           gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                           >  2:7272900/1‑150 (MQ=255)
                                            ggggccGGTGCGGGGCACCGCTAATGGGTTAATATTTTGGCGCCGGATATACGCCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGTGCGGGTTCGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                          <  2:1701635/145‑1 (MQ=255)
                                            ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                          >  1:2253651/1‑150 (MQ=255)
                                            ggCGACGGTGCTGGGGACCGCGATTGCGTTATGCTTGTGGCGCAGGACATATCGCAAACCGATCGCAAAATTGAGCGCTATCGGCCCCATCGCCCCGCGGGGGCGCGGGTTCGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGccc                                                                                                          <  2:6075720/150‑1 (MQ=255)
                                               gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                       >  1:1001974/1‑150 (MQ=255)
                                               gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                       >  1:538168/1‑150 (MQ=255)
                                                ccGGTGCTGGTGCCAGGGATGGCGTTAAGATGGTGGCGCAGGTGTTAGGGCAACCCGAGGGCACCTTTGCGCGCGATGGGTCCCCTCGCCCCTGTGGGGAGCGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                                      <  2:5128863/149‑1 (MQ=255)
                                                  ggagCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCACAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                                    >  2:4749468/4‑150 (MQ=255)
                                                   gtggtggtgACCGGGATTGGGTTATGTTTGTTGCGCTGGTGTTACGGCAACCCGTTGGCAACTTTGAGCGGGATGGGTCCCCTCGCCCCTGTGGGGTGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaa                                                                                                   <  2:8355093/150‑1 (MQ=255)
                                                         gTGACCGCTACTGTGTTATGTTGGTGGCGCTGTAGTTACGGCAACCCGACCGCAAATTGGCGCGCGCTCGGCCCCCGCGCCCCTGTGGGGCGTGGGTTGGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                             <  1:5536744/150‑1 (MQ=255)
                                                          tGACAGCGAATGCGTGAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtg                                                                                            >  2:2326189/1‑149 (MQ=255)
                                                              cgCGTATGCGTTATTTTTGTGGCGCAGGCTTTACGGAAAACCTATCGAAACTTTGCGCGCGATGGGGCCCCTCGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                        <  1:8280235/149‑1 (MQ=255)
                                                                 gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                                     >  2:6557795/1‑150 (MQ=255)
                                                                    tGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAATGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTAACGAATCaa                                                                                  >  1:2119737/1‑150 (MQ=255)
                                                                       gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                               >  2:3738581/1‑150 (MQ=255)
                                                                        ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                              >  1:6348027/1‑150 (MQ=255)
                                                                         tatgATTGTGGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGGCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCCATCAACCGCCTTTACGAATCAAATAAc                                                                             <  2:2822008/147‑1 (MQ=39)
                                                                               tgttgCGCAGGAGATATGGCAAACCGATCGAAACTTTGAGCGCGATGGGTCCCCTCGCCCCTGTGGGTAGTGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                       <  2:4598555/150‑1 (MQ=255)
                                                                                gttcCGCAGGACATAAGCCAAACCGATCGCAACATTGACCGCGATCGGACCCCTCGCCCCTCTGGGCAGAGGGTTAGGGTCAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                      <  2:2908902/150‑1 (MQ=255)
                                                                                 ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGc                                                                     >  2:4287229/1‑150 (MQ=255)
                                                                                  tgggcAGGAGTTAAGGAAAACGGATCGCAACATTGCGCGCTTGCGGACCCCCCGCCCCTGTGGGGTGTGTGTTTGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                    <  2:2397290/147‑1 (MQ=255)
                                                                                    ctctgtcgaTACGGAAACCCTATCTCAACTTTGCGCGGTTTCGGCCCCCGCGCCCCTGTGGGTTTTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                  <  2:3891099/143‑1 (MQ=255)
                                                                                    cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATAAGGCAGt                                                                  >  1:4936370/1‑150 (MQ=255)
                                                                                      cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                                >  1:7194427/1‑150 (MQ=255)
                                                                                             tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                         >  1:1192106/1‑150 (MQ=255)
                                                                                               aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                       >  1:6210046/1‑150 (MQ=255)
                                                                                               aGGAAACCCGATCGCAACTTTGAGCGCTAGCGGTCCCCCCGCCCCTCTGGGTTGTGGGTTTGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                       <  1:1385103/150‑1 (MQ=255)
                                                                                                ggCACACCGAGCGCAACTTGGCGCGCGCTGGGCCCCCTCGCCCGTGTGGGGAGTGGGTTGGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                      <  2:7609849/150‑1 (MQ=255)
                                                                                                ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                      >  1:4351733/1‑150 (MQ=255)
                                                                                                 gCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                     >  2:451391/1‑150 (MQ=255)
                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAGGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                     >  2:8069632/1‑150 (MQ=255)
                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                     >  1:6958945/1‑150 (MQ=255)
                                                                                                 gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCGTGACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                     >  2:2044985/1‑150 (MQ=255)
                                                                                                     aCCGATCGAAACATTGAGCGCGATCGCACCCCACCCCCCTCTGGGGAGAGGTTTAGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                                 <  1:2756420/150‑1 (MQ=255)
                                                                                                        gttCGCAACATTTGGGGGGTTGGGGCCCCGCCCCCCGGGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                              <  1:6849223/148‑1 (MQ=255)
                                                                                                        gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                              >  1:6396760/1‑150 (MQ=255)
                                                                                                        gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                              >  1:2152498/1‑150 (MQ=255)
                                                                                                          tCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                            <  1:5635592/150‑1 (MQ=255)
                                                                                                           cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                           >  2:3578953/1‑150 (MQ=255)
                                                                                                             cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                         >  2:5544976/1‑150 (MQ=255)
                                                                                                              aaCTTTGCGCGCGATCGGGCCCCCCGCCCCGGTGGGGTGCGGGTTTGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                        <  1:3819286/150‑1 (MQ=255)
                                                                                                              aaCATGGCGCGCGATGGGCCCCCGCGCCCCTGGGGGGTGGGGTTGCGGGTGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                        <  1:2437032/150‑1 (MQ=255)
                                                                                                                 aTTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                     <  2:1361019/150‑1 (MQ=255)
                                                                                                                  ttGAGCGCTTTCGGTCCCCTCGCCCCTGTGGGGAGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                    <  1:8443908/150‑1 (MQ=255)
                                                                                                                  ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                    >  1:2863807/1‑150 (MQ=255)
                                                                                                                    gcgcgcgcTCGGTCCCCGCGCCCGGCTGGGGTGTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                                  <  2:3372599/148‑1 (MQ=255)
                                                                                                                     tgcgcTATGGGCCCCCTCGCCCCTCGGGGGAGAGGGTTTGGGGGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                 <  2:131178/149‑1 (MQ=255)
                                                                                                                     aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                 >  2:2531612/1‑150 (MQ=255)
                                                                                                                          gagcGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGTTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                            <  2:3222749/147‑1 (MQ=255)
                                                                                                                           ttgggcccccccGCCCCTGTGGTGAGTGGTTTTGGGTGGGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                           <  1:3029441/144‑1 (MQ=255)
                                                                                                                           accctaacccTCTCCCCAGAGGGGCGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCCACCAGTTCGTGTTCCATAAAg                           >  2:8508899/8‑150 (MQ=255)
                                                                                                                           aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                           >  2:6993873/1‑150 (MQ=255)
                                                                                                                           aTCGGTCCCCTCGCCCCTCGGGGGAGCGGGTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                           <  2:6263764/150‑1 (MQ=255)
                                                                                                                              ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                        >  1:7558134/1‑150 (MQ=255)
                                                                                                                               gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCACATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                       >  1:5912674/1‑150 (MQ=255)
                                                                                                                                tCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                      >  2:2106271/1‑150 (MQ=255)
                                                                                                                                  cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                    >  1:4182671/1‑150 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                  >  1:5135794/1‑150 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                  >  2:7497300/1‑150 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                  >  2:3171604/1‑150 (MQ=255)
                                                                                                                                      cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTCCAAAAccc                >  1:1549116/1‑150 (MQ=255)
                                                                                                                                       gCCCGTCTGGGGTGTGGTTTTGGGGGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa               <  1:7352/150‑1 (MQ=255)
                                                                                                                                       gCCCCTGTGGGTTGTGTGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa               <  2:419664/150‑1 (MQ=255)
                                                                                                                                       gCCCCGCGGGGGTGTGGTTTTGGGGGTGGCGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa               <  1:8027164/150‑1 (MQ=255)
                                                                                                                                       gCCCATCTGGGGAGTGGTTTTGGTTGAGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa               <  2:5342043/150‑1 (MQ=255)
                                                                                                                                         cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg             >  2:4035336/1‑150 (MQ=255)
                                                                                                                                         cccTCGGGGTTGTGGGTTCGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg             <  2:6672497/150‑1 (MQ=255)
                                                                                                                                           ctctGGGGAGCGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa           <  1:4585402/150‑1 (MQ=255)
                                                                                                                                            tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat          >  1:4724/1‑150 (MQ=255)
                                                                                                                                            tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat          >  2:299559/1‑150 (MQ=255)
                                                                                                                                              tGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAatat        >  1:903122/1‑150 (MQ=255)
                                                                                                                                                gggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcg      >  1:1066519/1‑150 (MQ=255)
                                                                                                                                                 ggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc     >  1:6758488/1‑150 (MQ=255)
                                                                                                                                                  gtgagGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg    <  1:8559802/148‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                       <  1:1838894/83‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgcg    >  1:7795587/1‑150 (MQ=255)
                                                                                                                                                   tgagGGTTTGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg   <  1:316129/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGG  >  NC_000913/3608554‑3608851

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: