Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,169 A→T 34.6% I128N (ATC→AAC)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,1690AT34.6% 144.9 / 126.4 104I128N (ATC→AAC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (39/29);  new base T (15/21);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 105 total aligned reads displayed.

TTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCG  >  NC_000913/2867073‑2867267
                                                                                                |                                                                                                  
ttCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTc                                                                                                  >  2:1677532/1‑99 (MQ=255)
   aTAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGc                                                                                               <  2:150712/99‑1 (MQ=255)
      aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcccc                                                                                            >  1:1634996/1‑99 (MQ=255)
      aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcccc                                                                                            >  2:1712774/1‑99 (MQ=255)
       tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCa                                                                                           >  2:469786/1‑99 (MQ=255)
       tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCa                                                                                           >  2:2418040/1‑99 (MQ=255)
        cGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAg                                                                                          <  1:512330/99‑1 (MQ=255)
         gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAgg                                                                                         <  1:2331325/99‑1 (MQ=255)
          cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                        <  1:2146782/99‑1 (MQ=255)
          cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                        <  2:206061/99‑1 (MQ=255)
            cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTg                                                                                      >  1:2181652/1‑99 (MQ=255)
            cGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGCCAAACTTCCCTACGGCGCGGATCAGCCCCGGGtgg                                                                                      >  2:1339715/1‑97 (MQ=255)
              ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGcc                                                                                    <  1:1907318/99‑1 (MQ=255)
              ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGcc                                                                                    >  1:847145/1‑99 (MQ=255)
              ttCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGcc                                                                                    >  1:311128/1‑99 (MQ=255)
               tCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGccc                                                                                   >  1:1027690/1‑99 (MQ=255)
                  aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCctct                                                                                <  2:2241887/99‑1 (MQ=255)
                  aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctct                                                                                <  2:1496829/99‑1 (MQ=255)
                   tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCtt                                                                               <  2:2093539/99‑1 (MQ=255)
                   tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCtt                                                                               <  2:930444/99‑1 (MQ=255)
                      tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGTAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                            >  2:618115/1‑99 (MQ=255)
                      tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGa                                                                            >  1:700031/1‑99 (MQ=255)
                      tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                            >  1:1137724/1‑99 (MQ=255)
                       cTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGAt                                                                           <  2:1772627/99‑1 (MQ=255)
                         ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATaa                                                                         >  2:1072011/1‑99 (MQ=255)
                          gCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAg                                                                        <  1:932667/99‑1 (MQ=255)
                            gAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                                      <  1:1725437/99‑1 (MQ=255)
                            gAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGt                                                                      >  1:896582/1‑99 (MQ=255)
                              aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                                    >  2:758367/1‑99 (MQ=255)
                              aTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTcc                                                                    <  2:1492376/99‑1 (MQ=255)
                               tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                                   <  1:1547831/99‑1 (MQ=255)
                                ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                                  <  2:1009845/99‑1 (MQ=255)
                                 caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                                 <  2:2264848/99‑1 (MQ=255)
                                  accaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTCCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCa                                                                >  1:1485347/1‑99 (MQ=255)
                                   ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                               >  1:29010/1‑99 (MQ=255)
                                   ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                               >  2:1351727/1‑99 (MQ=255)
                                    caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAAc                                                              <  1:431205/99‑1 (MQ=255)
                                      ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                            >  2:1268391/1‑99 (MQ=255)
                                      ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                            <  1:1696699/99‑1 (MQ=255)
                                        ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                          <  2:1947534/99‑1 (MQ=255)
                                         tGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAg                                                         <  2:1348267/99‑1 (MQ=255)
                                             tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCa                                                     >  1:978814/1‑99 (MQ=255)
                                             tATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCa                                                     <  2:1513952/99‑1 (MQ=255)
                                                gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                  <  1:1669150/99‑1 (MQ=255)
                                                gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                  <  2:1784901/99‑1 (MQ=255)
                                                gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                  <  1:1070031/99‑1 (MQ=255)
                                                gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                  <  2:2188426/99‑1 (MQ=255)
                                                 ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                                 <  1:1456955/99‑1 (MQ=255)
                                                    agaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                              >  1:2449418/1‑99 (MQ=255)
                                                    agaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGc                                              >  1:299880/1‑99 (MQ=255)
                                                     gaAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGcc                                             >  2:2236818/1‑99 (MQ=255)
                                                       aGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCAt                                           >  1:805496/1‑99 (MQ=255)
                                                        gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                                          >  2:1770930/1‑99 (MQ=255)
                                                        gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATa                                          >  2:1316356/1‑99 (MQ=255)
                                                         cGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                                         <  2:804266/99‑1 (MQ=255)
                                                         cGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATaa                                         <  2:1336988/99‑1 (MQ=255)
                                                           gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                       >  2:703908/1‑99 (MQ=255)
                                                           gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                       >  2:911389/1‑99 (MQ=255)
                                                           gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                       >  1:2365555/1‑99 (MQ=255)
                                                            aaaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcgg                                      >  1:1101809/1‑99 (MQ=255)
                                                             aaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                                     >  1:1756424/1‑99 (MQ=255)
                                                             aaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                                     <  1:771687/99‑1 (MQ=255)
                                                              aCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                                    <  2:2175837/99‑1 (MQ=255)
                                                                  cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCa                                >  1:241095/1‑99 (MQ=255)
                                                                  cGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCa                                >  1:1704813/1‑99 (MQ=255)
                                                                   gTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                               <  2:1792754/99‑1 (MQ=255)
                                                                   gTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                               >  2:665155/1‑99 (MQ=255)
                                                                   gTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                               <  1:392537/99‑1 (MQ=255)
                                                                    ttCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                              >  2:891148/1‑99 (MQ=255)
                                                                    ttCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                              >  2:690068/1‑99 (MQ=255)
                                                                       cGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtttt                           <  2:535189/99‑1 (MQ=255)
                                                                        gggTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAttttt                          <  1:317190/99‑1 (MQ=255)
                                                                         ggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                         >  1:1238864/1‑99 (MQ=255)
                                                                         ggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                         <  2:288433/99‑1 (MQ=255)
                                                                          gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAGGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                        >  2:749433/1‑99 (MQ=255)
                                                                          gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                        >  1:1841228/1‑99 (MQ=255)
                                                                          gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                        >  1:309267/1‑99 (MQ=255)
                                                                           tCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                       >  2:879540/1‑99 (MQ=255)
                                                                           tCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                       <  1:966474/99‑1 (MQ=255)
                                                                            cAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca                      <  1:1351727/99‑1 (MQ=255)
                                                                              aaCTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                    >  1:2133026/1‑99 (MQ=255)
                                                                              aaCTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                    >  2:1290571/1‑99 (MQ=255)
                                                                              aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                    >  2:1274738/1‑99 (MQ=255)
                                                                               aCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca                   >  2:2217630/1‑99 (MQ=255)
                                                                                 ttCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa                 <  2:859756/99‑1 (MQ=255)
                                                                                  tctcTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc                <  2:1095083/99‑1 (MQ=255)
                                                                                      tACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa            >  1:1848908/1‑99 (MQ=255)
                                                                                      tACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa            >  2:778894/1‑99 (MQ=255)
                                                                                       aCCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg           >  1:2247461/1‑99 (MQ=255)
                                                                                        ccGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt          <  2:1906710/99‑1 (MQ=255)
                                                                                        ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt          <  1:2001641/99‑1 (MQ=255)
                                                                                          gcgcgGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTa        <  2:501463/99‑1 (MQ=255)
                                                                                          gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTa        <  1:2271784/99‑1 (MQ=255)
                                                                                           cgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAc       >  2:1149377/1‑99 (MQ=255)
                                                                                           cgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAc       <  1:1836092/99‑1 (MQ=255)
                                                                                            gcgGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAct      <  1:2081359/99‑1 (MQ=255)
                                                                                            gcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAct      <  2:69351/99‑1 (MQ=255)
                                                                                             cgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActc     >  2:788358/1‑99 (MQ=255)
                                                                                               gATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACtctc   <  1:1328955/99‑1 (MQ=255)
                                                                                                aTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCg  >  1:1403800/1‑99 (MQ=255)
                                                                                                |                                                                                                  
TTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCG  >  NC_000913/2867073‑2867267

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