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

ACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGA  >  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)
          atgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCgg‑g                                                                                      <  2:1222116/100‑1 (MQ=255)
           tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGACGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑Gt                                                                                     >  1:538413/1‑100 (MQ=255)
            gatgGCGCCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTg                                                                                    >  2:1353043/1‑100 (MQ=255)
            gatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTg                                                                                    >  1:922823/1‑100 (MQ=255)
            gatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTg                                                                                    >  1:1688887/1‑100 (MQ=255)
              tgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTg                                                                                    <  1:2175432/98‑1 (MQ=255)
               ggCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGc                                                                                 <  1:1030418/100‑1 (MQ=255)
                gCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCt                                                                                <  2:314871/100‑1 (MQ=255)
                 cGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAATGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑TTGAg                                                                                  <  1:2158212/97‑1 (MQ=255)
                 cGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACTCTCGCGG‑GTGAGCTg                                                                               >  2:2558950/1‑100 (MQ=255)
                  gACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGccag                                                                              <  2:2129847/100‑4 (MQ=255)
                   aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGc                                                                             <  2:1906543/100‑1 (MQ=255)
                    ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCt                                                                            >  2:346564/1‑100 (MQ=255)
                       cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAc                                                                         >  1:1538631/1‑100 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgcc                                                                      >  1:1717192/1‑100 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCCAGCTCAcccg                                                                      <  1:1532041/100‑4 (MQ=255)
                          tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCCAGCTCAcccg                                                                      <  1:2104730/100‑4 (MQ=255)
                           aTCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgccg                                                                     <  2:222155/100‑1 (MQ=255)
                            tCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgc                                                                     <  2:1473287/100‑5 (MQ=255)
                              gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCtt                                                                  >  2:883387/1‑100 (MQ=255)
                              gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCtt                                                                  >  2:1022608/1‑100 (MQ=255)
                               cTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCttt                                                                 <  2:1523717/100‑1 (MQ=255)
                                 ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCg                                                               <  1:32041/100‑1 (MQ=255)
                                  tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGt                                                              <  2:2543061/100‑1 (MQ=255)
                                   gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGtt                                                             >  2:704053/1‑100 (MQ=255)
                                   gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGtt                                                             >  2:704043/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCCTCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  1:1062853/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  2:672431/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  2:1619361/1‑100 (MQ=255)
                                    cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTa                                                            >  1:497011/1‑100 (MQ=255)
                                     ccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCCCTTTCGttct                                                           >  1:2015594/1‑98 (MQ=255)
                                      cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                          >  2:1167632/1‑100 (MQ=255)
                                      cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGCGCTGGCTCACGCCGCTTTCGTTAtt                                                          >  1:2594454/1‑100 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGTTTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:160000/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:500777/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:808553/100‑1 (MQ=255)
                                       aCGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                         <  1:1821992/100‑1 (MQ=255)
                                         gCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAg                                                       <  2:1731496/100‑1 (MQ=255)
                                             aGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAcc                                                   >  2:904594/1‑100 (MQ=255)
                                               cgcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCg                                                 >  2:1153246/1‑100 (MQ=255)
                                                gcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCgg                                                >  2:2408057/1‑100 (MQ=255)
                                                gcccgcCGAGCTGGAAGGGGAAGATGCATGCTCATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAcc                                                   <  2:2345536/94‑1 (MQ=255)
                                                  gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGa                                              <  1:10502/100‑1 (MQ=255)
                                                  gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGa                                              <  1:926318/100‑1 (MQ=255)
                                                     cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                           <  2:519570/100‑1 (MQ=255)
                                                     cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                           >  2:2156538/1‑100 (MQ=255)
                                                       gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                         >  2:358321/1‑100 (MQ=255)
                                                       gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCAGCTTTCGTTATTCAGCACCCGGGAAACTg                                         >  1:1560542/1‑100 (MQ=255)
                                                        agCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGt                                        <  1:996523/100‑1 (MQ=255)
                                                         gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                       <  2:708487/100‑1 (MQ=255)
                                                         gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                       <  1:2372289/100‑1 (MQ=255)
                                                          cGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAg                                      <  1:1831777/100‑1 (MQ=255)
                                                           gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                     >  2:1375726/1‑100 (MQ=255)
                                                            ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                    >  2:2451167/1‑100 (MQ=255)
                                                            ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                    >  2:2045520/1‑100 (MQ=255)
                                                                 gggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTccc                               <  2:1372678/100‑1 (MQ=255)
                                                                  ctaaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcccc                              <  2:123485/98‑1 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCag                            >  1:123484/1‑99 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAc                            >  1:284291/1‑100 (MQ=255)
                                                                    aaaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTCTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAc                            <  1:1590307/100‑1 (MQ=255)
                                                                     aaaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcg                           >  2:2141947/1‑100 (MQ=255)
                                                                      aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                          >  2:1284006/1‑100 (MQ=255)
                                                                       aaccACGCTGATTAATCACTGCTTTAACACGCTCGCGG‑GTGACCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgcc                         >  2:1755411/3‑100 (MQ=255)
                                                                       aTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgcc                         >  2:2031694/1‑100 (MQ=255)
                                                                          cATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCt                      <  1:1316543/100‑1 (MQ=255)
                                                                          cATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCt                      <  1:2558277/100‑1 (MQ=255)
                                                                           aTGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCtt                     <  1:672428/100‑1 (MQ=255)
                                                                             gCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGGAAATGTAGATTTCCCCACGCCGCTTa                    >  1:2495542/1‑99 (MQ=255)
                                                                                gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGACACTGTAGATTTCCCCACGCCGCTTAAgcg                >  1:1123471/1‑100 (MQ=255)
                                                                                gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg                <  1:1641931/100‑1 (MQ=255)
                                                                                  ttCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcg              >  1:1448523/1‑100 (MQ=255)
                                                                                   tcatcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             <  2:782346/100‑1 (MQ=255)
                                                                                   tcatcaCTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             >  1:372456/1‑100 (MQ=255)
                                                                                   tcatcaCTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcgcgc             <  1:693044/100‑1 (MQ=255)
                                                                                    catcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGcgcgcg            >  1:210178/1‑100 (MQ=255)
                                                                                    catcaCTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGcgcgcg            <  1:306656/100‑1 (MQ=255)
                                                                                     atcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGa           >  1:1287104/1‑100 (MQ=255)
                                                                                       caCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGata         >  1:2273893/1‑100 (MQ=255)
                                                                                          tGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCt      <  1:1534323/100‑1 (MQ=255)
                                                                                           gCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  1:1752466/1‑100 (MQ=255)
                                                                                           gCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:1097955/1‑100 (MQ=255)
                                                                                           gCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:1560162/1‑100 (MQ=255)
                                                                                           gCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt     >  2:2023083/1‑100 (MQ=255)
                                                                                            cTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCttt    >  2:2560844/1‑100 (MQ=255)
                                                                                              tCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGa  >  1:1690793/1‑100 (MQ=255)
                                                                                               |                                                                                                   
ACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTTTGA  >  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: ‑