Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,204 G→A 76.8% intergenic (‑105/+128) 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,2040GA76.8% 202.0 / 64.2 112intergenic (‑105/+128)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (15/11);  new base A (45/41);  total (60/52)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 115 total aligned reads displayed.

TCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAT  >  NC_000913/4196107‑4196299
                                                                                                 |                                                                                               
tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTcaca                                                                                               <  1:15958/100‑1 (MQ=255)
tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTcaca                                                                                               <  2:2157885/100‑1 (MQ=255)
tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTcaca                                                                                               <  1:687658/100‑1 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAg                                                                                              >  1:1004735/1‑100 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAg                                                                                              >  1:341332/1‑100 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAg                                                                                              >  1:193240/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAgg                                                                                             <  1:353936/100‑1 (MQ=255)
    aaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTa                                                                                           <  2:1111338/100‑1 (MQ=255)
    aaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTa                                                                                           <  1:1881375/100‑1 (MQ=255)
     aaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTAt                                                                                          <  1:1313360/100‑1 (MQ=255)
      aaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGATGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTAtt                                                                                         <  1:2517407/100‑1 (MQ=255)
        gTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAg                                                                                       >  1:2568636/1‑100 (MQ=255)
        gTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAg                                                                                       <  1:143169/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGcc                                                                                     <  2:1278900/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGcc                                                                                     >  2:1605968/1‑100 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                     >  2:2246277/1‑100 (MQ=255)
           aaGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCt                                                                                    <  2:2785921/100‑1 (MQ=255)
             ggACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGa                                                                                  <  2:1011779/100‑1 (MQ=255)
             ggACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGa                                                                                  >  1:1901276/1‑100 (MQ=255)
              gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGGGATGCGATAGCAGATTAAATCTTGTTCCCTTCACAGGTATTAGCCTGAt                                                                                 >  2:1791866/1‑100 (MQ=255)
              gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTg                                                                                   >  2:3025856/1‑98 (MQ=255)
                cGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCa                                                                               <  1:2015065/100‑1 (MQ=255)
                cGTGGCGTGTCAGACAACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTAGTAGCCTGACCa                                                                               <  2:918389/100‑1 (MQ=255)
                 gTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAg                                                                              <  1:892884/100‑1 (MQ=255)
                   ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGt                                                                            <  1:482688/100‑1 (MQ=255)
                       tGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCcg                                                                        <  2:737172/100‑1 (MQ=255)
                       tGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCcg                                                                        >  2:1434503/1‑100 (MQ=255)
                         tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCcgcg                                                                      <  2:193238/100‑1 (MQ=255)
                          cAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCgg                                                                     <  2:2498869/100‑1 (MQ=255)
                            gacgacGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGAt                                                                   >  2:2030480/1‑100 (MQ=255)
                            gacgacGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGAt                                                                   >  1:2198492/1‑100 (MQ=255)
                              cgacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATcc                                                                 <  2:1431793/100‑1 (MQ=255)
                              cgacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATcc                                                                 <  1:3041590/100‑1 (MQ=255)
                               gactGATGAAACAAGAGACGATCGCCGCGGGGGCGCTGCGATAGCATATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATccc                                                                <  2:3012067/100‑1 (MQ=255)
                                acgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCg                                                               <  1:3076617/100‑1 (MQ=255)
                                 cgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGa                                                              >  1:2674535/1‑100 (MQ=255)
                                  ggATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGaa                                                             <  1:115642/100‑1 (MQ=255)
                                   gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAAt                                                            >  2:1680960/1‑100 (MQ=255)
                                   gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAAt                                                            <  1:1005631/100‑1 (MQ=255)
                                       aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAAc                                                        <  2:1355722/100‑1 (MQ=255)
                                         gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACgg                                                      <  2:2812435/100‑1 (MQ=255)
                                         gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACgg                                                      >  1:2507680/1‑100 (MQ=255)
                                           aaGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCTAATAAACGGtc                                                    <  1:1392968/100‑1 (MQ=255)
                                           aaGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtc                                                    >  1:213255/1‑100 (MQ=255)
                                            agagACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct                                                   >  1:973110/1‑100 (MQ=255)
                                              agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                  <  2:19830/99‑1 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                              >  2:2271326/1‑100 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                              >  2:1841210/1‑100 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCACGAATAAACGGTCTCAGcc                                              >  2:757390/1‑100 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATACCGAATAAACGGTATCa                                                 >  2:2587999/1‑97 (MQ=255)
                                                  gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa                                             <  1:2687580/100‑1 (MQ=255)
                                                  gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa                                             >  2:2681602/1‑100 (MQ=255)
                                                   aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAg                                            >  1:890078/1‑100 (MQ=255)
                                                    tCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGt                                           <  1:2546861/100‑1 (MQ=255)
                                                      gccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTa                                         >  1:1657809/1‑100 (MQ=255)
                                                         gcAGGGGCGATGCGATAGCAGATTAACCCTTGTTCCCTTCGCAGGTATTAGCCTTATCAGGTTCCGCGGTTCCCGAATAAACGGTCTCAGCCAGTTAAgg                                      <  1:2912556/100‑1 (MQ=255)
                                                           aGGGGCGATGCGGTAGCAGATGAACTCTTGTTCCCTTCGCAGGTATTAGCGTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAgg                                      <  2:434845/98‑1 (MQ=255)
                                                            ggggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCAc                                   >  2:2342643/1‑100 (MQ=255)
                                                             gggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACt                                  >  1:1691940/1‑100 (MQ=255)
                                                               gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTcc                                >  2:92881/1‑100 (MQ=255)
                                                               gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTcc                                <  1:2377504/100‑1 (MQ=255)
                                                               gcgatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGACTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTcc                                >  2:1673541/1‑100 (MQ=255)
                                                                 gatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGa                              >  1:2499157/1‑100 (MQ=255)
                                                                   tgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGCATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGAc                             >  1:960299/1‑99 (MQ=255)
                                                                   tgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACa                            <  1:622148/100‑1 (MQ=255)
                                                                    gcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACaa                           <  1:664308/100‑1 (MQ=255)
                                                                     cgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAg                          >  1:163515/1‑100 (MQ=255)
                                                                         ctcAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAt                      <  1:2804154/98‑1 (MQ=255)
                                                                          gCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                     <  1:291378/100‑1 (MQ=255)
                                                                          gCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                     >  2:2132243/1‑100 (MQ=255)
                                                                          gCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                     >  1:423086/1‑100 (MQ=255)
                                                                           cAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtag                    >  2:2804125/1‑98 (MQ=255)
                                                                           cAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTa                    <  1:2805299/100‑1 (MQ=255)
                                                                            aGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                   >  1:700242/1‑100 (MQ=255)
                                                                             gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                  <  2:338178/100‑1 (MQ=255)
                                                                              aTTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcc                 <  1:1764654/100‑1 (MQ=255)
                                                                               ttAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                <  1:2027519/100‑1 (MQ=255)
                                                                               ttAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                <  2:2990834/100‑1 (MQ=255)
                                                                               ttAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                <  2:2627065/100‑1 (MQ=255)
                                                                                tAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcccc               >  1:555351/1‑100 (MQ=255)
                                                                                 aaCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCg              >  1:1965617/1‑100 (MQ=255)
                                                                                 aaCTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCg              <  2:423084/100‑1 (MQ=255)
                                                                                    tctTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaa           >  2:2214041/1‑100 (MQ=255)
                                                                                    tctTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaa           >  2:1710909/1‑100 (MQ=255)
                                                                                       tGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg        >  1:1333067/1‑100 (MQ=255)
                                                                                       tGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg        >  1:783142/1‑100 (MQ=255)
                                                                                        gTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg       >  1:483755/1‑100 (MQ=255)
                                                                                        gTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg       >  2:686485/1‑100 (MQ=255)
                                                                                         ttCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg      >  1:1498831/1‑100 (MQ=255)
                                                                                         ttCCCTTCACAGGTATTAGCCTGATCACGTTCCGCGGATCCCGAATAAACTGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg      <  2:2926438/100‑1 (MQ=255)
                                                                                          tCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg     <  1:2900362/100‑1 (MQ=255)
                                                                                          tCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg     >  2:2937810/1‑100 (MQ=255)
                                                                                          tCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg     <  2:2293739/100‑1 (MQ=255)
                                                                                          tCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg     <  2:154936/100‑1 (MQ=255)
                                                                                           cccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg     >  2:2338362/1‑99 (MQ=255)
                                                                                           cccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc    >  2:1864186/1‑100 (MQ=255)
                                                                                            ccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCCGTTAAGGCACTCCGGCAAGAAATTAGCCCCGCAAATGGGGCa   >  2:2182749/1‑100 (MQ=255)
                                                                                            ccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACCCCGACACGAAATTAGCCCCGCAAATggg      >  2:1471488/1‑97 (MQ=255)
                                                                                             cTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa   >  2:2801627/1‑99 (MQ=255)
                                                                                             cTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAt  >  1:413423/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAT  >  NC_000913/4196107‑4196299

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: