Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,196,250 T→C intergenic (‑151/+82) thiC ← / ← rsd phosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2500TC97.0% 445.2 / ‑0.3 132intergenic (‑151/+82)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (75/53);  total (78/54)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

AGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTT  >  NC_000913/4196153‑4196346
                                                                                                 |                                                                                                
agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCa                                                                                                >  1:1321174/1‑100 (MQ=255)
 gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAg                                                                                               >  1:1526345/1‑100 (MQ=255)
  aCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAATGGTCCCAGc                                                                                              <  1:2064519/100‑1 (MQ=255)
   cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGcc                                                                                             <  2:948410/100‑1 (MQ=255)
   cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGcc                                                                                             >  2:1237849/1‑100 (MQ=255)
   cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGcc                                                                                             >  2:1820428/1‑100 (MQ=255)
    gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCa                                                                                            <  2:1708191/100‑1 (MQ=255)
      tCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGt                                                                                          >  1:1309053/1‑100 (MQ=255)
        gccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTa                                                                                        <  2:2621928/100‑1 (MQ=255)
         ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTaa                                                                                       >  1:2242467/1‑100 (MQ=255)
          cgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAg                                                                                      <  1:399862/100‑1 (MQ=255)
          cgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAg                                                                                      >  1:1499555/1‑100 (MQ=255)
           gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAgg                                                                                     >  2:2170290/1‑100 (MQ=255)
           gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAgg                                                                                     <  1:9428/100‑1 (MQ=255)
           gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAgg                                                                                     <  1:2282672/100‑1 (MQ=255)
           gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAgg                                                                                     <  2:620726/100‑1 (MQ=255)
           gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAgg                                                                                     >  1:1614169/1‑100 (MQ=255)
            cAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGc                                                                                    <  2:2016899/100‑1 (MQ=255)
              ggggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCAc                                                                                  <  2:1862777/100‑1 (MQ=255)
                 gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTcc                                                                               <  2:2340536/100‑1 (MQ=255)
                 gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTcc                                                                               <  1:1737344/100‑1 (MQ=255)
                 gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTcc                                                                               <  2:1526332/100‑1 (MQ=255)
                    atgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGAc                                                                            <  2:1579754/100‑1 (MQ=255)
                     tgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACa                                                                           <  1:280010/100‑1 (MQ=255)
                       cgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAg                                                                         >  2:345560/1‑100 (MQ=255)
                       cgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAg                                                                         >  1:147589/1‑100 (MQ=255)
                         atAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGaa                                                                       >  2:836267/1‑100 (MQ=255)
                         atAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGaa                                                                       >  1:2246816/1‑100 (MQ=255)
                            gCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAAtt                                                                    >  1:2261783/1‑100 (MQ=255)
                            gCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAAtt                                                                    >  2:2022261/1‑100 (MQ=255)
                            gCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAAtt                                                                    <  1:1333659/100‑1 (MQ=255)
                            gCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAAtt                                                                    >  1:916197/1‑100 (MQ=255)
                             cAGATTAACTCTTGTTCCCTTCGCAGGTATTATCCTGATAAGCTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTa                                                                   <  2:2126408/100‑1 (MQ=255)
                             cAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTa                                                                   >  2:515842/1‑100 (MQ=255)
                              aGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                  <  2:1349098/100‑1 (MQ=255)
                              aGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                  >  1:2789412/1‑100 (MQ=255)
                              aGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                  <  2:775320/100‑1 (MQ=255)
                               gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                                                                 >  1:1128400/1‑100 (MQ=255)
                               gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                                                                 >  1:2197469/1‑100 (MQ=255)
                                  tAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcccc                                                              <  2:2180722/100‑1 (MQ=255)
                                   aaCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCg                                                             <  2:916189/100‑1 (MQ=255)
                                   aaCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCg                                                             >  1:450439/1‑100 (MQ=255)
                                    aCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc                                                            <  2:280399/100‑1 (MQ=255)
                                     ctctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTACCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCa                                                           >  1:2134229/1‑100 (MQ=255)
                                      tctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaa                                                          >  1:2934484/1‑100 (MQ=255)
                                       ctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa                                                         >  1:749914/1‑100 (MQ=255)
                                       ctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa                                                         >  2:1104726/1‑100 (MQ=255)
                                         tGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg                                                       >  1:459450/1‑100 (MQ=255)
                                          gTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg                                                      >  2:590643/1‑100 (MQ=255)
                                          gTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg                                                      >  1:669663/1‑100 (MQ=255)
                                             cccTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc                                                   >  1:2095852/1‑100 (MQ=255)
                                             cccTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc                                                   >  2:2865027/1‑100 (MQ=255)
                                              ccTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa                                                  >  2:392867/1‑100 (MQ=255)
                                               cTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAt                                                 >  1:1232303/1‑100 (MQ=255)
                                               cTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAt                                                 <  1:1766391/100‑1 (MQ=255)
                                                 tCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTg                                               <  1:2169369/100‑1 (MQ=255)
                                                   gCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGaa                                             <  1:1920577/100‑1 (MQ=255)
                                                   gCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGaa                                             <  1:66654/100‑1 (MQ=255)
                                                   gCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGaa                                             <  1:761221/100‑1 (MQ=255)
                                                      ggTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGt                                          >  2:2071333/1‑100 (MQ=255)
                                                       gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTa                                         >  2:2877386/1‑100 (MQ=255)
                                                       gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATCGGGCATTGAATGTa                                         >  2:2064414/1‑100 (MQ=255)
                                                        tATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTaa                                        >  1:789468/1‑100 (MQ=255)
                                                          ttAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAAt                                      <  1:680785/100‑1 (MQ=255)
                                                           tAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAAtt                                     >  2:279246/1‑100 (MQ=255)
                                                            aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTa                                    >  1:2526890/1‑100 (MQ=255)
                                                            aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTa                                    >  2:2915661/1‑100 (MQ=255)
                                                            aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTa                                    >  1:828767/1‑100 (MQ=255)
                                                             gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAc                                   <  1:185307/100‑1 (MQ=255)
                                                             gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAc                                   >  1:1386000/1‑100 (MQ=255)
                                                             gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAc                                   >  2:2248719/1‑100 (MQ=255)
                                                               cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgc                                 >  1:2657116/1‑100 (MQ=255)
                                                               cTGATCAGGTTCCGCGGATCCCGAATAAACGGACCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgc                                 >  1:2874479/1‑100 (MQ=255)
                                                                tGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgcg                                <  1:2794675/100‑1 (MQ=255)
                                                                tGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgcg                                >  1:886289/1‑100 (MQ=255)
                                                                     aGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAAc                           <  2:975219/100‑1 (MQ=255)
                                                                      ggTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACa                          >  2:2434842/1‑100 (MQ=255)
                                                                      ggTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACa                          >  2:378850/1‑100 (MQ=255)
                                                                           cgcgGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                     >  2:2172646/1‑100 (MQ=255)
                                                                           cgcgGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                     >  2:1383040/1‑100 (MQ=255)
                                                                           cgcgGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                     <  1:1424772/100‑1 (MQ=255)
                                                                           cgcgGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                     >  2:114059/1‑100 (MQ=255)
                                                                           cgcgGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                     >  1:2222501/1‑100 (MQ=255)
                                                                               gATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAc                 <  2:1499543/100‑1 (MQ=255)
                                                                                aTCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACt                <  1:2064433/100‑1 (MQ=255)
                                                                                 tCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTc               >  2:836215/1‑100 (MQ=255)
                                                                                   ccGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCTGAACTCaa             >  1:2593179/1‑100 (MQ=255)
                                                                                    cGAATAAACGGTCCCAGCCAGTTAAGGCACTCCTACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAg            >  2:2001188/1‑100 (MQ=255)
                                                                                    cGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAg            <  2:95157/100‑1 (MQ=255)
                                                                                    cGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAg            <  1:1102761/100‑1 (MQ=255)
                                                                                     gAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGc           >  2:1132082/1‑100 (MQ=255)
                                                                                     gAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGc           <  2:902646/100‑1 (MQ=255)
                                                                                      aaTAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGTATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCa          <  1:2143642/100‑1 (MQ=255)
                                                                                         aaaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGa       >  1:2912230/1‑100 (MQ=255)
                                                                                           aCGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATg     <  2:907844/100‑1 (MQ=255)
                                                                                            cGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGt    >  1:1901639/1‑100 (MQ=255)
                                                                                            cGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGt    >  1:1387636/1‑100 (MQ=255)
                                                                                            cGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGt    <  1:117124/100‑1 (MQ=255)
                                                                                             ggTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGtt   <  1:836274/100‑1 (MQ=255)
                                                                                              gTCCCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt  >  2:1351501/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTT  >  NC_000913/4196153‑4196346

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