Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,330 G→T 18.7% C25F (TGC→TTC)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3300GT18.7% 449.6 / 79.8 175C25F (TGC→TTC) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (62/80);  new base T (18/15);  total (80/95)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 175 total aligned reads displayed.

ATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAA  >  NC_000913/2560233‑2560428
                                                                                                 |                                                                                                  
attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                                   >  2:3850128/1‑99 (MQ=255)
 ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                 >  2:918946/1‑100 (MQ=255)
 ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                 >  2:1476275/1‑100 (MQ=255)
  tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgg                                                                                                >  2:190921/1‑100 (MQ=255)
   cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                               <  1:3322778/100‑1 (MQ=38)
   cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                               <  2:3213649/100‑1 (MQ=255)
   cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                               >  1:44308/1‑100 (MQ=255)
   cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                               >  1:3131084/1‑100 (MQ=255)
   cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                               >  1:2251819/1‑100 (MQ=255)
   cAATAACAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCTGCATGTTGCAGCATTCCggg                                                                                               <  2:1443723/100‑1 (MQ=255)
    aaTATCAAGAACTTAAATCATGTTTGCAACGAAAGATCTAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCTGGCATGTTGCAGCATTCCggg                                                                                               <  2:3626762/98‑1 (MQ=255)
     atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                             <  2:1561246/100‑1 (MQ=38)
     atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                             <  1:1840529/100‑1 (MQ=38)
     atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                             <  2:4170251/100‑1 (MQ=255)
      tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                            <  1:3458822/100‑1 (MQ=38)
      tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                            >  1:701514/1‑100 (MQ=255)
       atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCAGAATGCTGCAACATGCCGGGGGCg                                                                                           <  2:3133341/100‑1 (MQ=255)
        tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                          <  1:2013486/100‑1 (MQ=255)
        tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                          <  1:1334504/100‑1 (MQ=255)
        tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                          <  2:2356829/100‑1 (MQ=255)
          aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                        <  2:3159288/100‑1 (MQ=255)
          aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                        <  2:1589072/100‑1 (MQ=255)
          aaGAACTTAAATCCATGAATGCAACGAAAGATCCAGAATTCGAAAAAAGGAACAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGCCGGTa                                                                                        <  2:1345058/100‑1 (MQ=255)
           aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                                       <  1:820006/100‑1 (MQ=255)
            gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                      <  1:277904/100‑1 (MQ=255)
             aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                     <  2:2859203/100‑1 (MQ=255)
               cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAACAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                   <  1:2628035/100‑1 (MQ=255)
                ttAAATCCATGTTTGCAACGAAAGAACCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttttatttgt                                                                                  <  2:569061/100‑11 (MQ=21)
                  aaaTCCGTGTTTGCAACGAAAGATCCATAGCTCGAAAATAGGATCAACACAAATAAAAGACCCCGCAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                  <  2:3955685/98‑1 (MQ=255)
                     tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                             >  2:3065487/1‑100 (MQ=255)
                      ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                            <  2:3798181/100‑1 (MQ=255)
                          gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                        >  2:2984606/1‑100 (MQ=255)
                          gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                        <  2:550868/100‑1 (MQ=255)
                             tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGga                                                                     <  2:3020006/100‑1 (MQ=255)
                               cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                                   >  1:3760211/1‑100 (MQ=255)
                                aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                  >  1:3921014/1‑100 (MQ=255)
                                aaCGAAAGATCCAGAGTTCGAAAAGAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                                   <  2:1851726/99‑1 (MQ=255)
                                 aCGAAAGAGCCAGAGTAAGAAAAGAGGATCAACACAAGTAAAAGCCCCCAGAATGCTGCAACATGCCGGGGGCAGTATGAGCAACCAGCTAAAGGAGAAt                                                                 <  2:1211081/100‑1 (MQ=255)
                                     aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                            <  1:1485468/100‑1 (MQ=255)
                                      aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                            <  2:3523741/100‑1 (MQ=255)
                                      aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                            <  2:1715043/100‑1 (MQ=255)
                                        aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                          >  1:2381181/1‑100 (MQ=255)
                                        aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                          >  2:775454/1‑100 (MQ=255)
                                        aTCCAAAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                          <  1:1721496/100‑1 (MQ=255)
                                         tCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATg                                                         <  1:2691366/100‑1 (MQ=255)
                                           cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGt                                                        >  2:1799944/1‑99 (MQ=255)
                                              agTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                                    <  2:3423700/100‑1 (MQ=255)
                                               gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                   >  1:850309/1‑100 (MQ=255)
                                                ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                                  <  1:3820613/100‑1 (MQ=255)
                                                ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCAGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGGCTGata                                                  >  2:1953331/1‑100 (MQ=255)
                                                 tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                                 <  1:961578/100‑1 (MQ=255)
                                                   gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                               <  1:3986597/100‑1 (MQ=255)
                                                   gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                               <  1:589885/100‑1 (MQ=255)
                                                    aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                              >  1:922835/1‑100 (MQ=255)
                                                       aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                           <  2:1353024/100‑1 (MQ=255)
                                                        tAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                                          <  1:3060385/100‑1 (MQ=255)
                                                         aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                         <  1:1410203/100‑1 (MQ=255)
                                                          ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                        <  2:1975799/100‑1 (MQ=255)
                                                          agctCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                                          <  2:1640643/95‑1 (MQ=255)
                                                            aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                                      >  2:600792/1‑100 (MQ=255)
                                                               aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                   <  2:2215957/100‑1 (MQ=255)
                                                                acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGATATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                                  >  2:277019/1‑100 (MQ=255)
                                                                  ttctcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTATGATATGCTAGCAGTAGAACAAg                                >  2:2419915/11‑100 (MQ=21)
                                                                  acaAATAAAAGCCCCCGGAATGTTGCAACATGCCGGGGGCGGTATGAGTAACAAGCTAACGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                                >  2:1871148/1‑100 (MQ=255)
                                                                  acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                                <  2:2973752/100‑1 (MQ=255)
                                                                  acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                                >  1:2216984/1‑100 (MQ=255)
                                                                   caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                               <  1:1945377/100‑1 (MQ=255)
                                                                     aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                             >  1:380101/1‑100 (MQ=255)
                                                                      aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                            <  1:1776482/100‑1 (MQ=255)
                                                                      aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                            <  2:2564082/100‑1 (MQ=255)
                                                                      aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                            <  2:1328246/100‑1 (MQ=255)
                                                                       tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                           <  1:1672787/100‑1 (MQ=255)
                                                                         accgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa                         >  1:2804482/4‑100 (MQ=37)
                                                                          ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                        >  1:1631100/3‑100 (MQ=37)
                                                                            gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                      >  2:4021680/1‑100 (MQ=38)
                                                                            gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                      >  2:269247/1‑100 (MQ=38)
                                                                            gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                      >  1:2889047/1‑100 (MQ=38)
                                                                            gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCGAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                      >  2:3344710/1‑100 (MQ=38)
                                                                             cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGACTTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                     >  2:4108477/1‑100 (MQ=38)
                                                                             cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                     >  1:3941717/1‑100 (MQ=38)
                                                                             cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                     >  2:2698318/1‑100 (MQ=38)
                                                                              ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGACTTTCTCATGTCCGATATGCTAACAGTAGAACAAGCAACAAACAatg                    >  2:881362/1‑100 (MQ=25)
                                                                              ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg                    >  2:3459536/1‑100 (MQ=255)
                                                                               cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                   >  2:4115009/1‑100 (MQ=38)
                                                                                ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat                  >  1:1093381/1‑100 (MQ=255)
                                                                                ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat                  >  1:1174558/1‑100 (MQ=255)
                                                                                 cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg                 <  1:201562/100‑1 (MQ=255)
                                                                                 cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGGCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg                 <  1:773332/100‑1 (MQ=255)
                                                                                   gAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt               >  1:2637777/1‑100 (MQ=255)
                                                                                      tgctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGt            <  2:1649722/100‑1 (MQ=255)
                                                                                       gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc           >  1:3961108/1‑100 (MQ=255)
                                                                                        ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa          >  2:651952/1‑100 (MQ=255)
                                                                                         tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa         >  1:1080520/1‑100 (MQ=255)
                                                                                           cAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTGTGATATGCTAGCAATAGAACAAGAAACAAACAATGATGTTCGTCAAt        >  2:2023815/1‑99 (MQ=255)
                                                                                            aaCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttt      <  2:537026/100‑1 (MQ=255)
                                                                                             aCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtttt     >  2:495597/1‑100 (MQ=255)
                                                                                              cATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttttt    >  2:1974572/1‑100 (MQ=255)
                                                                                               aTGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa   <  2:2476722/100‑1 (MQ=255)
                                                                                               aTGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa   <  1:598079/100‑1 (MQ=255)
                                                                                                tGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaa  <  1:2984815/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
ATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAA  >  NC_000913/2560233‑2560428

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