Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,120 C→T 8.3% intergenic (+34/+131) 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,1200CT 8.3% 212.9 / 7.3 73intergenic (+34/+131)ybcW/ylcIDLP12 prophage; uncharacterized protein/DUF3950 family protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base T (3/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

TGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTT  >  NC_000913/579972‑580266
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                           <  3:75139/125‑1 (MQ=255)
                                                 gCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTct                                                                                                                           >  4:75139/1‑125 (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                                                          >  8:100004/1‑130 (MQ=255)
                                                                                                             gACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAAt                                                          <  7:100004/130‑1 (MQ=255)
                                                                                                                  aTCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGAAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                             >  6:26694/1‑138 (MQ=255)
                                                                                                                  aTCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGAAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTc                                             <  5:26694/138‑1 (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                                                                                                           <  7:194729/56‑1 (MQ=255)
                                                                                                                                      ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAg                                                                                                           >  8:194729/1‑56 (MQ=255)
                                                                                                                                       tCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAtt             >  5:116733/1‑149 (MQ=255)
                                                                                                                                        cAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAtt             <  6:241254/148‑1 (MQ=255)
                                                                                                                                         aaaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                                                                 >  4:130215/1‑79 (MQ=255)
                                                                                                                                         aaaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGc                                                                                 <  3:130215/79‑1 (MQ=255)
                                                                                                                                          aaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttctt          <  7:31769/149‑1 (MQ=255)
                                                                                                                                           actactTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAttcttc         <  1:201100/149‑1 (MQ=255)
                                                                                                                                                ttAACCTGTGACAGAGCCAACACGCAGTCTGTCAATGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCATCGTAATTATGTGAATTTACAATATCGTCGTGTTCGGAGGTAAGAACGCGGGATGTTCATTCTTCATCAc    >  4:230197/1‑149 (MQ=255)
                                                                                                                                                  aaCCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACtt  >  7:231790/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTT  >  NC_000913/579972‑580266

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: