Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,169 A→T 15.2% 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,1690AT15.2% 355.3 / 48.7 113I128N (ATC→AAC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base A (49/46);  new base T (10/6);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 115 total aligned reads displayed.

GGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTC  >  NC_000913/2867071‑2867266
                                                                                                  |                                                                                                 
ggTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGa                                                                                                   >  2:1192390/1‑99 (MQ=255)
   tCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCa                                                                                                <  1:2629228/99‑1 (MQ=255)
      tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcc                                                                                             >  1:3024512/1‑99 (MQ=255)
      tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcc                                                                                             >  1:466312/1‑99 (MQ=255)
      tAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcc                                                                                             >  2:2558856/1‑99 (MQ=255)
       aaTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGccc                                                                                            <  2:1356432/99‑1 (MQ=255)
        aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcccc                                                                                           >  1:2980345/1‑99 (MQ=255)
        aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcccc                                                                                           <  1:1412801/99‑1 (MQ=255)
        aTCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGcccc                                                                                           >  1:702211/1‑99 (MQ=255)
         tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCa                                                                                          >  2:2410210/1‑99 (MQ=255)
         tCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCa                                                                                          >  2:1105380/1‑99 (MQ=255)
           gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAgg                                                                                        <  2:128240/99‑1 (MQ=255)
           gCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAgg                                                                                        <  1:3014143/99‑1 (MQ=255)
            cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                       >  2:2723844/1‑99 (MQ=255)
            cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                       >  1:627949/1‑99 (MQ=255)
            cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                       >  1:2139069/1‑99 (MQ=255)
            cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                       >  2:735373/1‑99 (MQ=255)
            cccGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGGCAAACTTCTCTACCGCGCGGATCAGCCCCAGGt                                                                                       <  1:2800970/99‑1 (MQ=255)
               gTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGc                                                                                    >  2:1944906/1‑99 (MQ=255)
                  cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCct                                                                                 >  2:1712099/1‑99 (MQ=255)
                  cAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCct                                                                                 <  2:221898/99‑1 (MQ=255)
                   aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGCTCAGCCCCAGGTTGCCctc                                                                                >  2:102429/1‑99 (MQ=255)
                   aaTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctc                                                                                >  1:203576/1‑99 (MQ=255)
                    aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctct                                                                               >  2:2247609/1‑99 (MQ=255)
                    aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctct                                                                               <  1:630192/99‑1 (MQ=255)
                    aTCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCctct                                                                               <  1:2347656/99‑1 (MQ=255)
                     tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCtt                                                                              >  2:3060330/1‑99 (MQ=255)
                     tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCtt                                                                              >  1:2004582/1‑99 (MQ=255)
                     tCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCtt                                                                              <  1:1699967/99‑1 (MQ=255)
                        tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                           >  2:1749522/1‑99 (MQ=255)
                        tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                           >  1:2244253/1‑99 (MQ=255)
                        tCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGGTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGa                                                                           <  2:606485/99‑1 (MQ=255)
                          tGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATa                                                                         >  2:2618986/1‑99 (MQ=255)
                           ggCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATaa                                                                        <  2:387204/99‑1 (MQ=255)
                             cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                      <  1:609735/99‑1 (MQ=255)
                             cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                      <  1:2472421/99‑1 (MQ=255)
                             cGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAgg                                                                      <  2:192582/99‑1 (MQ=255)
                               aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc                                                                    >  1:2968521/1‑99 (MQ=255)
                               aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc                                                                    <  1:592677/99‑1 (MQ=255)
                               aaTCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTc                                                                    >  1:2145796/1‑99 (MQ=255)
                                 tCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCa                                                                  >  2:2482309/1‑99 (MQ=255)
                                  ccaccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAg                                                                 >  1:761971/1‑99 (MQ=255)
                                   caccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGc                                                                >  1:1017016/1‑99 (MQ=255)
                                     ccaGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCaa                                                              >  2:3010931/1‑99 (MQ=255)
                                      caGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAAc                                                             >  2:2937418/1‑99 (MQ=255)
                                       aGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACg                                                            >  2:1141595/1‑99 (MQ=255)
                                        ggTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGc                                                           <  2:145125/99‑1 (MQ=255)
                                         gTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGcc                                                          >  2:2985858/1‑99 (MQ=255)
                                          ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                         >  1:2511247/1‑99 (MQ=255)
                                          ttGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCa                                                         >  1:2947737/1‑99 (MQ=255)
                                                  gTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGa                                                 <  1:1231950/99‑1 (MQ=255)
                                                   ttGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGAt                                                <  2:2918631/99‑1 (MQ=255)
                                                        aaGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCa                                           <  1:2011425/99‑1 (MQ=255)
                                                          gCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTTCCATa                                         >  2:2290146/1‑99 (MQ=255)
                                                             gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                      >  1:221384/1‑99 (MQ=255)
                                                             gAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcg                                      >  2:2914522/1‑99 (MQ=255)
                                                              aaaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcgg                                     >  1:937191/1‑99 (MQ=255)
                                                               aaCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggc                                    <  2:1955785/99‑1 (MQ=255)
                                                                aCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                                   >  1:2357995/1‑99 (MQ=255)
                                                                aCCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAAcggcg                                   <  2:1739622/99‑1 (MQ=255)
                                                                   aCGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGc                                >  1:626828/1‑99 (MQ=255)
                                                                     gTTCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCaa                              >  1:1807643/1‑99 (MQ=255)
                                                                      ttCCGGGTCAAACTTCTCTACCGCGCGGCTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                             <  1:2572097/99‑1 (MQ=255)
                                                                      ttCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                             >  2:1993006/1‑99 (MQ=255)
                                                                      ttCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAt                             >  1:46484/1‑99 (MQ=255)
                                                                       tCCGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtt                            <  1:1607365/99‑1 (MQ=255)
                                                                         cGGGTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAtttt                          <  2:2961349/99‑1 (MQ=255)
                                                                           ggTCAAACTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                        >  1:1401484/1‑99 (MQ=255)
                                                                           ggTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTa                        <  2:1493821/99‑1 (MQ=255)
                                                                            gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                       <  2:2861327/99‑1 (MQ=255)
                                                                            gTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTac                       <  1:1819274/99‑1 (MQ=255)
                                                                             tCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                      >  2:1443854/1‑99 (MQ=255)
                                                                             tCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacc                      >  2:498727/1‑99 (MQ=255)
                                                                              cAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca                     <  2:1504882/99‑1 (MQ=255)
                                                                              cAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTacca                     <  2:774162/99‑1 (MQ=255)
                                                                               aaaCTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccac                    <  1:959627/99‑1 (MQ=255)
                                                                                aaCTTCTCTACCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                   >  2:2009502/1‑99 (MQ=255)
                                                                                aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                   >  1:2274669/1‑99 (MQ=255)
                                                                                aaCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                   >  1:2502908/1‑99 (MQ=255)
                                                                                 aCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca                  <  2:1569828/99‑1 (MQ=255)
                                                                                 aCTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTAccacca                  >  1:2307567/1‑99 (MQ=255)
                                                                                  cTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAg                 >  1:1613106/1‑99 (MQ=255)
                                                                                   ttCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa                <  2:2306610/99‑1 (MQ=255)
                                                                                   ttCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa                <  1:1463134/99‑1 (MQ=255)
                                                                                   ttCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGa                <  1:1951338/99‑1 (MQ=255)
                                                                                    tctcTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGAc               <  2:2759629/99‑1 (MQ=255)
                                                                                        tACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCaa           <  2:203576/99‑1 (MQ=255)
                                                                                         aCCGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAg          >  1:1294114/1‑99 (MQ=255)
                                                                                          ccGCGCGGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt         <  1:2539858/99‑1 (MQ=255)
                                                                                          ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt         <  1:203008/99‑1 (MQ=255)
                                                                                          ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt         <  1:2388353/99‑1 (MQ=255)
                                                                                           cgcgcgGTTCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGtt        >  2:843228/1‑99 (MQ=255)
                                                                                               cgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActc    >  1:2473722/1‑99 (MQ=255)
                                                                                               cgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActc    >  2:1523737/1‑99 (MQ=255)
                                                                                                ggATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActct   <  1:1653331/99‑1 (MQ=255)
                                                                                                ggATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActct   >  2:1335166/1‑99 (MQ=255)
                                                                                                ggATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActct   <  1:1417185/99‑1 (MQ=255)
                                                                                                ggATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActct   >  2:2611458/1‑99 (MQ=255)
                                                                                                 gATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACtctc  <  2:2271594/99‑1 (MQ=255)
                                                                                                 gATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACtctc  >  2:2084893/1‑99 (MQ=255)
                                                                                                  |                                                                                                 
GGTTCATAATCGCCCGTTCAATCGTCTGGCGAATCCACCAGGTTGCGTATGTTGAGAAGCGGAAACCACGTTCCGGGTCAAACTTCTCTACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTC  >  NC_000913/2867071‑2867266

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