Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,331 C→G 100% intergenic (‑232/+1) 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,3310CG95.6% 627.5 / 13.2 180intergenic (‑232/+1)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base G (82/90);  total (86/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 181 total aligned reads displayed.

GATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAG  >  NC_000913/4196232‑4196426
                                                                                                   |                                                                                               
gATCCCGAATAAACGGTCTCAGCCAGTTAAGGTACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAc                                                                                                 <  1:53104/100‑1 (MQ=255)
   cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCa                                                                                              <  2:3707594/100‑1 (MQ=255)
    ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCaa                                                                                             >  1:1544689/1‑100 (MQ=255)
    ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCaa                                                                                             >  1:1498988/1‑100 (MQ=255)
     cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAg                                                                                            <  1:2784978/100‑1 (MQ=255)
     cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAg                                                                                            >  1:1129949/1‑100 (MQ=255)
     cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAg                                                                                            <  2:1722565/100‑1 (MQ=255)
          aaaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGa                                                                                       >  2:1990260/1‑100 (MQ=255)
          aaaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGa                                                                                       <  2:85090/100‑1 (MQ=255)
           aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGAt                                                                                      >  1:2756734/1‑100 (MQ=255)
           aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGAt                                                                                      >  1:3641550/1‑100 (MQ=255)
             cGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGt                                                                                    <  1:2346722/100‑1 (MQ=255)
             cGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGt                                                                                    <  2:3776963/100‑1 (MQ=255)
               gTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGttt                                                                                  <  1:3650400/100‑1 (MQ=255)
                 ctcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGa                                                                                <  1:1268519/100‑1 (MQ=255)
                  tcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAc                                                                               >  1:3126014/1‑100 (MQ=255)
                   cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAcg                                                                              >  2:1438923/1‑100 (MQ=255)
                     gCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGAcgcg                                                                            >  2:2952216/1‑100 (MQ=255)
                        aGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATGACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgggc                                                                         <  2:1691015/100‑1 (MQ=255)
                        aGTTAAGACAGTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGCAGTCAAGCAGGATGTTTGACgcggg                                                                          >  2:2776491/1‑99 (MQ=255)
                         gTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgggcg                                                                        <  1:2989424/100‑1 (MQ=255)
                         gTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgggcg                                                                        <  1:2330858/100‑1 (MQ=255)
                          ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACgcgggcgg                                                                       <  2:3601688/100‑1 (MQ=255)
                           tAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGggcggc                                                                      >  2:1649082/1‑100 (MQ=255)
                            ctggCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                     <  1:2296386/98‑1 (MQ=255)
                            aaGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                     <  2:1879550/100‑1 (MQ=255)
                             aGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGgcggcg                                                                     >  2:2230120/1‑99 (MQ=255)
                             aGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGt                                                                    >  2:1255181/1‑100 (MQ=255)
                             aGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGt                                                                    <  1:2237675/100‑1 (MQ=255)
                             aGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGt                                                                    <  1:1361555/100‑1 (MQ=255)
                              ggCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGt                                                                    >  1:3734541/1‑99 (MQ=255)
                              ggCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTc                                                                   >  2:3196504/1‑100 (MQ=255)
                               gCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                  <  1:2892152/100‑1 (MQ=255)
                               gCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCa                                                                  >  1:2176985/1‑100 (MQ=255)
                                cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCaa                                                                 >  1:1159133/1‑100 (MQ=255)
                                cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCaa                                                                 <  1:3327968/100‑1 (MQ=255)
                                 aCTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCaaa                                                                <  1:349206/100‑2 (MQ=255)
                                  cTCCGACAAGAAATTAGCCCCGCAAATGGGGCATCGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGTCGGTGTCAAGc                                                               <  2:2068576/100‑1 (MQ=255)
                                   tCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                                                              <  1:2952234/100‑1 (MQ=255)
                                   tCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCAAACa                                                              <  1:1416937/100‑1 (MQ=255)
                                    ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc                                                             >  1:2308076/1‑100 (MQ=255)
                                    ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCAAACat                                                             <  1:2013730/100‑2 (MQ=255)
                                    ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCAAACat                                                             <  1:3457805/100‑2 (MQ=255)
                                    ccGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGCTTGACGCCGCCCGCGTCAAACat                                                             <  2:3347065/100‑2 (MQ=38)
                                     cGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCGAGCAcc                                                            <  2:2850780/100‑1 (MQ=255)
                                      gACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca                                                           >  1:3597371/1‑100 (MQ=255)
                                       aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                          <  2:3156309/100‑1 (MQ=255)
                                       aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCCGCCCGCGTCaaacatcct                                                          <  1:1214745/100‑8 (MQ=37)
                                        cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc                                                         >  1:3757256/1‑100 (MQ=255)
                                        cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc                                                         >  2:1500156/1‑100 (MQ=255)
                                         aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                                                        <  2:907379/100‑1 (MQ=255)
                                              nttAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATc                                                   >  1:1895048/2‑100 (MQ=255)
                                               ttAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCa                                                  <  1:2975338/100‑1 (MQ=255)
                                                 aGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGc                                                <  1:2204763/100‑1 (MQ=255)
                                                  gCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGc                                                <  2:3397825/99‑1 (MQ=255)
                                                  gCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                                               >  1:1107230/1‑100 (MQ=255)
                                                  gCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                                               <  1:2494088/100‑1 (MQ=255)
                                                  gCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                                               >  1:1424071/1‑100 (MQ=255)
                                                    cccGCAAATGGGGGATTGAATGTAAATTACGCGTTAACCGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa                                             >  2:2432585/1‑100 (MQ=255)
                                                    cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa                                             >  1:272285/1‑100 (MQ=255)
                                                    cccGCAAATGGGGCATTGAATGTAAATTACGCGTAAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa                                             >  1:1488077/1‑100 (MQ=255)
                                                       gCAAATGGGGCGTTGAATGTAAATTACGAGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATcg                                          <  2:882436/100‑2 (MQ=255)
                                                         aaaTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTc                                        <  1:1698295/100‑1 (MQ=255)
                                                             ggggCATTGAATGTAACTTACGCGTTAACAGCGCAGAAGTAGAGCAGGATGTTTGAAGCGGGCCGCGTCAAGCCCCAACAGAATCAGCTTATCTTACAGc                                    <  2:2299773/100‑1 (MQ=255)
                                                             ggggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCatc                                    <  1:2659001/100‑3 (MQ=255)
                                                              gggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGca                                   <  1:2085059/100‑1 (MQ=255)
                                                               ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcac                                  >  2:232499/1‑100 (MQ=255)
                                                                 cATTGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGcaca                                 >  1:1729369/1‑99 (MQ=255)
                                                                    tGAATGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgc                             <  1:1070889/100‑1 (MQ=255)
                                                                      aaTGTAAATTACGTGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc                           <  2:3188247/100‑1 (MQ=255)
                                                                      aaTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc                           >  2:3358798/1‑100 (MQ=255)
                                                                       aTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                          >  2:816620/1‑100 (MQ=255)
                                                                       aTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                          >  1:3813643/1‑100 (MQ=255)
                                                                       aTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcg                          >  2:1712597/1‑100 (MQ=255)
                                                                       aTGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc                           >  2:2996125/1‑99 (MQ=255)
                                                                        tGTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgcgc                         >  1:795763/1‑100 (MQ=255)
                                                                         gTAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                        >  2:974561/1‑100 (MQ=255)
                                                                          tAAATTACGCGTTAACAGCGCATAAGTCAAGCAGGATTTTTGACGCGGACGGCGTCAAGCACCAGGAGAATGAGCTTATCTTCCAGCACAAAGCGCGCtt                       >  2:1928909/1‑100 (MQ=255)
                                                                          tAAATTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCtt                       >  2:1904769/1‑100 (MQ=255)
                                                                          tAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCtt                       >  2:74713/1‑100 (MQ=255)
                                                                            aaTTACGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                     <  2:187244/100‑1 (MQ=255)
                                                                                aCGCGTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAgcg                 <  1:761466/100‑1 (MQ=255)
                                                                                   cgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCt              >  2:1771983/1‑100 (MQ=255)
                                                                                   cgTTAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCt              >  1:2994144/1‑100 (MQ=255)
                                                                                     ttAACAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCg            <  1:2869462/100‑1 (MQ=255)
                                                                                       aaCAGCGCAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGcc          <  2:1725930/100‑1 (MQ=255)
                                                                                        aCAGCGCAGAAGTCAAGCAGGATGTTTAACGCGGGAGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAACGAGCTTCCAGCGCCTCGc           >  2:43379/1‑98 (MQ=255)
                                                                                          aGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt       >  1:1310954/1‑100 (MQ=255)
                                                                                           gcgcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg      <  2:809017/100‑1 (MQ=255)
                                                                                           gcgcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg      >  1:2126213/1‑100 (MQ=255)
                                                                                           gagcAGAACTCCAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg      <  1:3547495/98‑1 (MQ=255)
                                                                                            cgcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGt     <  1:957908/100‑1 (MQ=255)
                                                                                             gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc    >  2:3521158/1‑100 (MQ=255)
                                                                                             gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc    <  1:3385075/100‑1 (MQ=255)
                                                                                             gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc    <  2:3676085/100‑1 (MQ=255)
                                                                                             gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc    >  1:3255796/1‑100 (MQ=255)
                                                                                             gcAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc    <  2:1786700/100‑1 (MQ=255)
                                                                                              cAGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCa   <  1:15837/100‑1 (MQ=255)
                                                                                               aGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg  <  1:722570/100‑1 (MQ=255)
                                                                                               aGAAGTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAg  <  2:652337/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAG  >  NC_000913/4196232‑4196426

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