Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,467,129 C→T 84.1% E40K (GAA→AAA)  treR ← trehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,467,1290CT84.1% 219.6 / 31.0 94E40K (GAA→AAA) treRtrehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (10/5);  new base T (44/35);  total (54/40)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGA  >  NC_000913/4467034‑4467227
                                                                                               |                                                                                                   
aCGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAaca                                                                                                 <  2:1898870/100‑1 (MQ=255)
aCGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAaca                                                                                                 >  1:1866104/1‑100 (MQ=255)
 cGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAacac                                                                                                >  2:1015226/1‑100 (MQ=255)
   ggTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACCCGc                                                                                              >  1:1634222/1‑100 (MQ=255)
   ggTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGc                                                                                              >  2:498285/1‑100 (MQ=255)
    gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCt                                                                                             >  2:926321/1‑100 (MQ=255)
      aaCAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcg                                                                                           >  2:1869911/1‑100 (MQ=255)
       aCAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgc                                                                                          <  2:862496/100‑1 (MQ=255)
         aaTGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCACAGCGGGAAGGGGAAAATCCATACTGATTCATCACTGCTTTAACACGCTCACgg                                                                                        >  1:98666/1‑100 (MQ=255)
          atgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCggg                                                                                      <  2:1222116/100‑1 (MQ=255)
           tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGACGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGt                                                                                     >  1:538413/1‑100 (MQ=255)
            gatgGCGCCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTg                                                                                    >  2:1353043/1‑100 (MQ=255)
            gatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTg                                                                                    >  1:922823/1‑100 (MQ=255)
            gatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTg                                                                                    >  1:1688887/1‑100 (MQ=255)
              tgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTg                                                                                    <  1:2175432/98‑1 (MQ=255)
               ggCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGc                                                                                 <  1:1030418/100‑1 (MQ=255)
                gCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCt                                                                                <  2:314871/100‑1 (MQ=255)
                 cGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAATGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGTTGAg                                                                                  <  1:2158212/97‑1 (MQ=255)
                 cGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACTCTCGCGGGTGAGCTg                                                                               >  2:2558950/1‑100 (MQ=255)
                  gACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGccag                                                                              <  2:2129847/100‑4 (MQ=255)
                   aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGc                                                                             <  2:1906543/100‑1 (MQ=255)
                    ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCt                                                                            >  2:346564/1‑100 (MQ=255)
                       cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAc                                                                         >  1:1538631/1‑100 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgcc                                                                      >  1:1717192/1‑100 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCCAGCTCAcccg                                                                      <  1:1532041/100‑4 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCCAGCTCAcccg                                                                      <  1:2104730/100‑4 (MQ=255)
                           aTCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgccg                                                                     <  2:222155/100‑1 (MQ=255)
                            tCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgc                                                                     <  2:1473287/100‑5 (MQ=255)
                              gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                                  >  2:883387/1‑100 (MQ=255)
                              gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                                  >  2:1022608/1‑100 (MQ=255)
                               cTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCttt                                                                 <  2:1523717/100‑1 (MQ=255)
                                 ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCg                                                               <  1:32041/100‑1 (MQ=255)
                                  tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGt                                                              <  2:2543061/100‑1 (MQ=255)
                                   gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGtt                                                             >  2:704053/1‑100 (MQ=255)
                                   gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGtt                                                             >  2:704043/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCCTCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  1:1062853/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  2:672431/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  2:1619361/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  1:497011/1‑100 (MQ=255)
                                     ccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCCCTTTCGttct                                                           >  1:2015594/1‑98 (MQ=255)
                                      cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                          >  2:1167632/1‑100 (MQ=255)
                                      cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGCGCTGGCTCACGCCGCTTTCGTTAtt                                                          >  1:2594454/1‑100 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGTTTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:160000/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:500777/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:808553/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:1821992/100‑1 (MQ=255)
                                         gCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAg                                                       <  2:1731496/100‑1 (MQ=255)
                                             aGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAcc                                                   >  2:904594/1‑100 (MQ=255)
                                               cgcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCg                                                 >  2:1153246/1‑100 (MQ=255)
                                                gcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCgg                                                >  2:2408057/1‑100 (MQ=255)
                                                gcccgcCGAGCTGGAAGGGGAAGATGCATGCTCATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAcc                                                   <  2:2345536/94‑1 (MQ=255)
                                                  gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGa                                              <  1:10502/100‑1 (MQ=255)
                                                  gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGa                                              <  1:926318/100‑1 (MQ=255)
                                                     cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                           <  2:519570/100‑1 (MQ=255)
                                                     cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                           >  2:2156538/1‑100 (MQ=255)
                                                       gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                         >  2:358321/1‑100 (MQ=255)
                                                       gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCAGCTTTCGTTATTCAGCACCCGGGAAACTg                                         >  1:1560542/1‑100 (MQ=255)
                                                        agCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGt                                        <  1:996523/100‑1 (MQ=255)
                                                         gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                       <  2:708487/100‑1 (MQ=255)
                                                         gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                       <  1:2372289/100‑1 (MQ=255)
                                                          cGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAg                                      <  1:1831777/100‑1 (MQ=255)
                                                           gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                     >  2:1375726/1‑100 (MQ=255)
                                                            ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                    >  2:2451167/1‑100 (MQ=255)
                                                            ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                    >  2:2045520/1‑100 (MQ=255)
                                                                 gggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTccc                               <  2:1372678/100‑1 (MQ=255)
                                                                  ctaaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcccc                              <  2:123485/98‑1 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCag                            >  1:123484/1‑99 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAc                            >  1:284291/1‑100 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTCTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAc                            <  1:1590307/100‑1 (MQ=255)
                                                                     aaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcg                           >  2:2141947/1‑100 (MQ=255)
                                                                      aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                          >  2:1284006/1‑100 (MQ=255)
                                                                       aaccACGCTGATTAATCACTGCTTTAACACGCTCGCGGGTGACCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgcc                         >  2:1755411/3‑100 (MQ=255)
                                                                       aTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgcc                         >  2:2031694/1‑100 (MQ=255)
                                                                          cATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCt                      <  1:1316543/100‑1 (MQ=255)
                                                                          cATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCt                      <  1:2558277/100‑1 (MQ=255)
                                                                           aTGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCtt                     <  1:672428/100‑1 (MQ=255)
                                                                             gCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGGAAATGTAGATTTCCCCACGCCGCTTa                    >  1:2495542/1‑99 (MQ=255)
                                                                                gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGACACTGTAGATTTCCCCACGCCGCTTAAgcg                >  1:1123471/1‑100 (MQ=255)
                                                                                gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg                <  1:1641931/100‑1 (MQ=255)
                                                                                  ttCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcg              >  1:1448523/1‑100 (MQ=255)
                                                                                   tcatcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             <  2:782346/100‑1 (MQ=255)
                                                                                   tcatcaCTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             >  1:372456/1‑100 (MQ=255)
                                                                                   tcatcaCTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             <  1:693044/100‑1 (MQ=255)
                                                                                    catcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGcgcgcg            >  1:210178/1‑100 (MQ=255)
                                                                                    catcaCTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGcgcgcg            <  1:306656/100‑1 (MQ=255)
                                                                                     atcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGa           >  1:1287104/1‑100 (MQ=255)
                                                                                       caCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGata         >  1:2273893/1‑100 (MQ=255)
                                                                                          tGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCt      <  1:1534323/100‑1 (MQ=255)
                                                                                           gCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  1:1752466/1‑100 (MQ=255)
                                                                                           gCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:1097955/1‑100 (MQ=255)
                                                                                           gCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:1560162/1‑100 (MQ=255)
                                                                                           gCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:2023083/1‑100 (MQ=255)
                                                                                            cTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCttt    >  2:2560844/1‑100 (MQ=255)
                                                                                              tCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGa  >  1:1690793/1‑100 (MQ=255)
                                                                                               |                                                                                                   
ACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGA  >  NC_000913/4467034‑4467227

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