Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,115 C→T 8.3% intergenic (+29/+136) ybcW → / ← ylcI DLP12 prophage; uncharacterized protein/DUF3950 family protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,1150CT 8.3% 219.5 / 7.4 72intergenic (+29/+136)ybcW/ylcIDLP12 prophage; uncharacterized protein/DUF3950 family protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base T (3/3);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTC  >  NC_000913/579967‑580259
                                                                                                                                                    |                                                                                                                                                
ccTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAactac                                                                                                                                                  >  7:94537/1‑149 (MQ=255)
ccTGATGTGGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAactac                                                                                                                                                  <  6:181993/149‑1 (MQ=255)
     tgttgtAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAAc                                                                                                                                             <  8:206224/149‑1 (MQ=255)
      gttgtAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAAc                                                                                                                                             <  5:83614/148‑1 (MQ=255)
         gtAAATGCAACTATAGAACATAAGGTGACTCTGGAAGCATACAGCGCAATTGAGGCAGCGTTGGTGCAGCACGATAATAATATGAAGGATAATTCCCAGGTCGTTGACTGATCACCATAACTGATAATCATTCAAACTACTTAACCagt                                                                                                                                         >  8:121275/1‑149 (MQ=255)
          tAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCtgtg                                                                                                                                        <  4:115070/149‑1 (MQ=255)
           aaTTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCtgtg                                                                                                                                        >  3:68465/1‑148 (MQ=255)
                 aTGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGc                                                                                                                                 <  8:178838/149‑1 (MQ=255)
                  tGTATATAACATAAGGTGTCTCTTGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAATGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGcc                                                                                                                                <  5:183148/149‑1 (MQ=255)
                   gTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGcc                                                                                                                                >  3:2132/1‑148 (MQ=255)
                     ataGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCaa                                                                                                                              >  3:206798/1‑148 (MQ=255)
                      taGATCATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAaca                                                                                                                            >  6:48911/1‑149 (MQ=255)
                      taGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGTCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCATGGAGGTTGAATGATCACAATAAATGCTAATCATTCAAACTACTTACACTGTGACAGAGacaaca                                                                                                                            >  8:62033/1‑149 (MQ=255)
                      taGAACATAAGGTGTCTCTGGAAGCATACACGCCAAATGACGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTCGTCGTTGACTGATCACCATAACTCCTAATCATTCAAACTACATAACCTGTGACAGAGCCAaca                                                                                                                            >  8:65547/1‑149 (MQ=255)
                          aCATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACg                                                                                                                          >  7:182798/1‑147 (MQ=255)
                             tAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTc                                                                                                                     >  7:206899/1‑149 (MQ=255)
                                 gtgtCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgt                                                                                                                  >  4:82827/1‑148 (MQ=255)
                                  tgtCTCTGGAAGCATTCACGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTAATCCCTGGTGGTTTACTCATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACACAGCCAACACGCAGTctgtca                                                                                                                <  4:526/149‑1 (MQ=255)
                                    tctcTTGAAGCATTCAGGGCAATTGAGGCATCGTTGGTGAAGCACGATAATAATATGAATGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcact                                                                                                              <  8:56641/149‑1 (MQ=255)
                                    tctcTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcact                                                                                                              <  5:202715/149‑1 (MQ=255)
                                      tctGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcactgt                                                                                                            >  6:51254/1‑149 (MQ=255)
                                      tctGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcactgt                                                                                                            >  4:178978/1‑149 (MQ=255)
                                      tctGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTctgtcactgt                                                                                                            >  5:210847/1‑149 (MQ=255)
                                        tGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACAATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcactgt                                                                                                            <  4:68465/147‑1 (MQ=255)
                                          gAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtcactgtca                                                                                                          <  8:41449/147‑1 (MQ=255)
                                             gCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGaaa                                                                                                     <  4:2132/149‑1 (MQ=255)
                                                ttCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTg                                                                                                  <  7:180839/149‑1 (MQ=255)
                                                ttCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTg                                                                                                  <  8:94537/149‑1 (MQ=255)
                                                  cAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTgg                                                                                                 <  4:112781/148‑1 (MQ=255)
                                                      gCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTct                                                                                                                    >  4:75139/1‑125 (MQ=255)
                                                      gCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTct                                                                                                                    <  3:75139/125‑1 (MQ=255)
                                                      gCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCACTGGTGGTTGACTGATCACCATAACTGCTAATCATTCGAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTaaaa                                                                                            <  3:72847/149‑1 (MQ=255)
                                                          ttGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTg                                                                                         >  5:206048/1‑148 (MQ=255)
                                                             atgCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTTGTGGGTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCaa                                                                                      <  3:146548/146‑1 (MQ=255)
                                                             aGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCaa                                                                                      <  5:48911/148‑1 (MQ=255)
                                                               gCAGCGTTGGTGAAGCAAGAAAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTc                                                                                   <  4:46827/149‑1 (MQ=255)
                                                                       ggTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACt                                                                            <  1:74092/148‑1 (MQ=255)
                                                                        gTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTg                                                                           <  6:206048/148‑1 (MQ=255)
                                                                               aCGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGccc                                                                   >  5:241254/1‑149 (MQ=255)
                                                                                  ataataataTGAAGGATTATTCCCTGGTTGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCg                                                                <  7:178737/149‑1 (MQ=255)
                                                                                  ataataataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCg                                                                >  2:201100/1‑149 (MQ=255)
                                                                                      taataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGaaa                                                                                                     >  5:240413/1‑108 (MQ=255)
                                                                                      taataTGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGaaa                                                                                                     <  6:240413/108‑1 (MQ=255)
                                                                                         tatGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTaa                                                         <  8:206899/149‑1 (MQ=255)
                                                                                           tGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGt                                                       <  6:210847/149‑1 (MQ=255)
                                                                                             aaGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                                      >  2:115815/1‑148 (MQ=255)
                                                                                             aaGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGa                                                     <  3:82827/149‑1 (MQ=255)
                                                                                                     ttCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAt                                                   <  2:26706/143‑1 (MQ=255)
                                                                                                     ttCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAt                                                   >  1:26706/1‑143 (MQ=255)
                                                                                                     ttCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACaa                                             <  1:115815/149‑1 (MQ=255)
                                                                                                     ttCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCATCGTAATTAAGTGAATTTACaa                                             >  3:156093/1‑149 (MQ=255)
                                                                                                        ccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAatat                                          >  4:16623/1‑149 (MQ=255)
                                                                                                            gtggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                      <  3:16623/149‑1 (MQ=255)
                                                                                                             tggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTcc                                     <  8:182798/149‑1 (MQ=255)
                                                                                                                ttGACTGATCACCATAACTGATAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCCACTCAATTACTGCAATGCCCGCGAAATTAAGTGAATTTACAATATCGTCCTGt                                  >  2:118111/1‑149 (MQ=255)
                                                                                                                  gACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAt                                                   <  7:100004/130‑1 (MQ=255)
                                                                                                                  gACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAt                                                   >  8:100004/1‑130 (MQ=255)
                                                                                                                       aTCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGAAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                      <  5:26694/138‑1 (MQ=255)
                                                                                                                       aTCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGAAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                      >  6:26694/1‑138 (MQ=255)
                                                                                                                          aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtca                                                                                                                >  7:155600/1‑61 (MQ=255)
                                                                                                                          aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTctgtca                                                                                                                <  8:155600/61‑1 (MQ=255)
                                                                                                                          aCCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGgaa                        <  3:118669/149‑1 (MQ=255)
                                                                                                                              tAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAc                    >  3:36574/1‑149 (MQ=255)
                                                                                                                                aCTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAcgc                  <  5:51254/149‑1 (MQ=255)
                                                                                                                                     tAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTAACTGTAAGGAAATTGGTAAATCTGCAACTCAATTACTGCAATGCCCTCGGAATTAAGTGAATTTACAATATCGTCCTGTTCGAAGGGAAGAAAGCGGGAt             >  8:226875/1‑149 (MQ=255)
                                                                                                                                       aTCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGt           >  4:88654/1‑149 (MQ=255)
                                                                                                                                          aTTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCa        >  2:36558/1‑149 (MQ=255)
                                                                                                                                           ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAg                                                                                                    >  8:194729/1‑56 (MQ=255)
                                                                                                                                           ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAg                                                                                                    <  7:194729/56‑1 (MQ=255)
                                                                                                                                            tCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAtt      >  5:116733/1‑149 (MQ=255)
                                                                                                                                             cAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAtt      <  6:241254/148‑1 (MQ=255)
                                                                                                                                              aaaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                                                          <  3:130215/79‑1 (MQ=255)
                                                                                                                                              aaaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                                                          >  4:130215/1‑79 (MQ=255)
                                                                                                                                               aaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttctt   <  7:31769/149‑1 (MQ=255)
                                                                                                                                                actactTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttcttc  <  1:201100/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTC  >  NC_000913/579967‑580259

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