Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,972,637 2 bp→GC coding (430‑431/546 nt) wecD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6370CG99.7% 1444.5 / NA 376R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (184/191);  total (184/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,972,6380GC100.0% 1455.9 / NA 377R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (184/193);  total (184/193)

Only 100 of 390 total aligned reads displayed.

GAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAA  >  NC_000913/3972547‑3972727
                                                                                          ||                                                                                         
gAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgc                                                                                       >  1:9949672/1‑96 (MQ=255)
gAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGAAGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGTGCGGATGCTcgc                                                                                       >  2:4840282/1‑96 (MQ=255)
 aaCTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgcg                                                                                      <  1:10030883/96‑1 (MQ=255)
 aaCTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgcg                                                                                      >  2:3374083/1‑96 (MQ=255)
 aaCTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGGGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgcg                                                                                      >  2:4963329/1‑96 (MQ=255)
   cTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGt                                                                                    >  2:9594334/1‑96 (MQ=255)
    tCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                   <  2:5392464/96‑1 (MQ=255)
    tCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                   >  1:13668860/1‑96 (MQ=255)
     cAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaa                                                                                  >  2:14312231/1‑96 (MQ=255)
       aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaa                                                                                >  1:9040657/1‑96 (MQ=255)
        tGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaac                                                                               <  2:13632130/96‑1 (MQ=255)
         gCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaaca                                                                              >  2:6500559/1‑96 (MQ=255)
          cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaa                                                                             >  1:11183462/1‑96 (MQ=255)
           naCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                                            >  1:1911629/2‑96 (MQ=255)
           gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                                            <  2:2187294/96‑1 (MQ=255)
            aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACt                                                                           >  2:2558466/1‑96 (MQ=255)
             cAGATGCGCGTATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACtt                                                                          <  1:12205737/96‑1 (MQ=255)
             cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGGGCTAAACTGGGCGTATGCTCGCGGTGAAACAACtt                                                                          >  2:13418447/1‑96 (MQ=255)
              aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                                         <  2:14151155/96‑1 (MQ=255)
              aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                                         >  1:5610841/1‑96 (MQ=255)
                  gcgcGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCgg                                                                     <  1:8898479/96‑1 (MQ=255)
                    gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                                                   <  1:4993180/96‑1 (MQ=255)
                    gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGGAAAACAACTTTGCGGGt                                                                   >  1:422905/1‑96 (MQ=255)
                     cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                                  <  2:452228/96‑1 (MQ=255)
                     cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCGAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                                  >  2:5468823/1‑96 (MQ=255)
                      gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTgg                                                                 >  1:6261931/1‑96 (MQ=255)
                      gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTgg                                                                 >  2:12041440/1‑96 (MQ=255)
                      gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCCTCGCGGTAAAACAACTTTGCGGGTgg                                                                 <  2:2781540/96‑1 (MQ=255)
                         ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                                              <  1:8639914/96‑1 (MQ=255)
                         ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                                              <  1:7046764/96‑1 (MQ=255)
                             ccTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCa                                                          >  1:3206412/1‑96 (MQ=255)
                              ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaa                                                         >  2:13041861/1‑96 (MQ=255)
                               tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACTCaaa                                                        <  1:12604683/96‑1 (MQ=255)
                                  tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                     <  2:8237538/96‑1 (MQ=255)
                                   ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                    <  1:7650783/96‑1 (MQ=255)
                                     ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGGGGCGACCCAAATGGGCa                                                  >  2:6756661/1‑96 (MQ=255)
                                     ctggACGCAGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTCTGCGGGTGGCGACCCAAATGGGCa                                                  <  1:4690085/96‑1 (MQ=255)
                                      tggGCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                                 >  2:10123323/1‑96 (MQ=255)
                                      tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                                 >  1:5411117/1‑96 (MQ=255)
                                         aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGcaccac                                              >  2:11015232/1‑96 (MQ=255)
                                         aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                              >  2:624456/1‑96 (MQ=255)
                                          cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                             <  2:11717948/96‑1 (MQ=255)
                                           gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                            <  1:7300176/96‑1 (MQ=255)
                                           gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                            >  1:2457121/1‑96 (MQ=255)
                                             ggTGCAGGTGCTGAGCTGATTCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACACAAATGGGCAACAccgcc                                          >  1:9283948/1‑96 (MQ=255)
                                             ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgcc                                          >  1:1292747/1‑96 (MQ=255)
                                             ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgcc                                          <  1:4377685/96‑1 (MQ=255)
                                               tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                        <  1:13263489/96‑1 (MQ=255)
                                                  aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                     >  2:1468232/1‑96 (MQ=255)
                                                  aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                     >  2:5134470/1‑96 (MQ=255)
                                                   ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                    >  2:5969724/1‑96 (MQ=255)
                                                    gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                   <  2:7302464/96‑1 (MQ=255)
                                                      gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                 >  2:12494269/1‑96 (MQ=255)
                                                      gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                 >  1:130670/1‑96 (MQ=255)
                                                      gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                 >  1:3416433/1‑96 (MQ=255)
                                                       ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                                <  1:8399764/96‑1 (MQ=255)
                                                       ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                                <  1:406736/96‑1 (MQ=255)
                                                        tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                               <  2:3740754/96‑1 (MQ=255)
                                                        tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                               >  1:2660434/1‑96 (MQ=255)
                                                        tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                               >  2:312674/1‑96 (MQ=255)
                                                         gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGa                              <  2:8519039/96‑1 (MQ=255)
                                                           gctgaTGCAAACGGCTCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGATGTCGACCCAAATGGGCAACACAGCTCGGCTGAAACGaga                            >  1:959356/1‑94 (MQ=255)
                                                           gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                            <  1:5078862/96‑1 (MQ=255)
                                                            ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                           >  2:2878284/1‑96 (MQ=255)
                                                              gaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacat                         <  2:356187/96‑1 (MQ=255)
                                                               aTGCAAACGGCGCTAACCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                        <  1:4849006/96‑1 (MQ=255)
                                                                tGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                       <  1:14376372/96‑1 (MQ=255)
                                                                tGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                       <  2:2428208/96‑1 (MQ=255)
                                                                tGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                       >  2:4437938/1‑96 (MQ=255)
                                                                 gCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                      <  2:9774996/96‑1 (MQ=255)
                                                                 gCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                      >  2:11945404/1‑96 (MQ=255)
                                                                  cAAACGGCGCTAAACGGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACCCCGCCGCGCTTAAACGATACATACaa                     <  1:11571502/96‑1 (MQ=255)
                                                                   aaaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                    >  2:13532803/1‑96 (MQ=255)
                                                                   aaaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                    >  2:11271390/1‑96 (MQ=255)
                                                                    aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                   >  2:9266595/1‑96 (MQ=255)
                                                                    aaCGGCGCTAAACTGGGCGTATGCTCGCGGCAAAACAACTTTGCGGGTGGCGCCCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                   >  2:4166871/1‑96 (MQ=255)
                                                                     aCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAATACCGCCGCGCTTAAACGATACATACAAAgt                  >  1:12746923/1‑96 (MQ=255)
                                                                     aCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt                  <  2:13280575/96‑1 (MQ=255)
                                                                     aCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt                  <  1:6682846/96‑1 (MQ=255)
                                                                      cGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg                 >  1:1524074/1‑96 (MQ=255)
                                                                       ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg                >  2:7135633/1‑96 (MQ=255)
                                                                       ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg                >  1:5768701/1‑96 (MQ=255)
                                                                         cgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg              >  1:6124286/1‑96 (MQ=255)
                                                                          gcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAATGGTGc             >  2:1589495/1‑96 (MQ=255)
                                                                           cTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCg            >  1:3763147/1‑96 (MQ=255)
                                                                            tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa           >  1:2720368/1‑96 (MQ=255)
                                                                            tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa           >  2:77782/1‑96 (MQ=255)
                                                                            tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa           <  1:1207597/96‑1 (MQ=255)
                                                                             aaaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAATACCGCCGCGCGTAAACGATACATACAAAGTGGTGCGaa          >  2:3948427/1‑96 (MQ=255)
                                                                              aaCTGGGCGTATGCTCGCGGTAATACAACTTTGCGGGTGGCGACACAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt         <  1:10123323/96‑1 (MQ=255)
                                                                              aaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGCGGTGCGAAt         >  2:8947149/1‑96 (MQ=255)
                                                                               gcTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg        <  2:8091207/95‑1 (MQ=255)
                                                                               aCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg        >  2:12744241/1‑96 (MQ=255)
                                                                                 ttggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa      >  1:10547618/3‑96 (MQ=255)
                                                                                 tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa      <  2:7587568/96‑1 (MQ=255)
                                                                                 tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa      <  1:7903974/96‑1 (MQ=255)
                                                                                 tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa      >  2:3025834/1‑96 (MQ=255)
                                                                                  gggCGTATGCTCGCCGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAg     <  1:10303014/96‑1 (MQ=255)
                                                                                   ggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGa    <  1:13377688/96‑1 (MQ=255)
                                                                                     cGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGaaa  >  1:14164372/1‑96 (MQ=255)
                                                                                          ||                                                                                         
GAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAA  >  NC_000913/3972547‑3972727

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