Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,467,129 C→T 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,1290CT99.2% 418.4 / ‑4.9 122E40K (GAA→AAA) treRtrehalose 6‑phosphate‑inducible trehalose regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (60/61);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 122 total aligned reads displayed.

GACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTT  >  NC_000913/4467033‑4467224
                                                                                                |                                                                                                 
gACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAac                                                                                                <  1:821297/100‑1 (MQ=255)
gACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAac                                                                                                <  1:547073/100‑1 (MQ=255)
 aCGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAaca                                                                                               <  1:156077/100‑1 (MQ=255)
  cGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAacac                                                                                              >  2:220330/1‑100 (MQ=255)
  cGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAacac                                                                                              <  1:1648138/100‑1 (MQ=255)
    ggTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGc                                                                                            >  2:615311/1‑100 (MQ=255)
    ggTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGc                                                                                            >  2:952151/1‑100 (MQ=255)
     gTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCt                                                                                           >  2:2104828/1‑100 (MQ=255)
         gatctatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgcg                                                                                       <  2:913095/95‑1 (MQ=255)
         cAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgcg                                                                                       <  1:2341812/100‑1 (MQ=255)
         cAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgc                                                                                        >  2:859574/1‑99 (MQ=255)
          aaTGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCGATGCTGATTCATCACTGCTTTAACACGCTCGCgg                                                                                      <  2:956700/100‑1 (MQ=255)
          aaTGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCgg                                                                                      <  2:1694219/100‑1 (MQ=255)
          aaTGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCgg                                                                                      <  2:1960690/100‑1 (MQ=255)
           atgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAACATCCATGCTGATTCATCACTGCTTTAACACGCTCGCggg                                                                                    >  2:2639194/1‑100 (MQ=255)
            tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGt                                                                                   >  1:1343168/1‑100 (MQ=255)
            tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGt                                                                                   >  2:1239836/1‑100 (MQ=255)
              atgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgcgagatc                                                                                  >  1:913097/1‑95 (MQ=255)
                   gACCACTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGc                                                                           <  2:1884009/100‑1 (MQ=255)
                    aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGc                                                                           <  1:1689975/100‑1 (MQ=255)
                    aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGc                                                                           <  1:42320/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCt                                                                          <  1:1507292/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCt                                                                          <  1:1232248/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCt                                                                          <  2:515407/100‑1 (MQ=255)
                      cACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTc                                                                         <  1:848590/100‑1 (MQ=255)
                        cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAc                                                                       >  2:2713408/1‑100 (MQ=255)
                        cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAc                                                                       >  1:999590/1‑100 (MQ=255)
                          ttATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgc                                                                     <  1:742043/100‑1 (MQ=255)
                          ttATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgc                                                                     <  2:1152253/100‑1 (MQ=255)
                           tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgcc                                                                    >  2:2768347/1‑100 (MQ=255)
                           tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCAcgcc                                                                    <  1:1493157/100‑1 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                                >  2:1141557/1‑100 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                                >  1:1504914/1‑100 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCtt                                                                >  1:792708/1‑100 (MQ=255)
                                 tttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTc                                                              <  1:1899429/100‑1 (MQ=255)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgcgagcg                                                              <  2:959148/100‑10 (MQ=37)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgcgagcg                                                              <  1:1272463/100‑10 (MQ=37)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCg                                                             <  1:344470/100‑1 (MQ=255)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGAAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCg                                                             <  1:20820/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGt                                                            <  2:2455653/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGt                                                            <  2:1527060/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGt                                                            >  2:37546/1‑100 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGt                                                            >  2:2681823/1‑100 (MQ=255)
                                    gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGtt                                                           >  2:2724365/1‑100 (MQ=255)
                                    gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGtt                                                           >  2:582896/1‑100 (MQ=255)
                                     cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTa                                                          >  2:2575089/1‑100 (MQ=255)
                                      ccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTAt                                                         <  2:2870082/100‑1 (MQ=255)
                                       cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                        >  1:27511/1‑100 (MQ=255)
                                       cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                        >  2:2835907/1‑100 (MQ=255)
                                        tctcATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                       <  2:1437794/97‑1 (MQ=255)
                                         cGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCa                                                      <  2:2469958/100‑1 (MQ=255)
                                           cATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAga                                                    >  1:1437800/1‑99 (MQ=255)
                                           cATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGc                                                    <  1:1995602/100‑1 (MQ=255)
                                              aGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTAGCTCACGCCGCTTTCGTTATTCAGCAcc                                                 >  1:2515845/1‑100 (MQ=255)
                                               gcgcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAccc                                                >  1:2379865/1‑100 (MQ=255)
                                                  cgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCggg                                             >  2:653503/1‑100 (MQ=255)
                                                  cgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCggg                                             <  2:362973/100‑1 (MQ=255)
                                                   gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTTACGCCGCTTTCGTTATTCAGCACCCGGGa                                            <  2:2364205/100‑1 (MQ=255)
                                                      cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                         >  1:676277/1‑100 (MQ=255)
                                                       agagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACt                                        >  2:390117/1‑100 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       <  2:262098/100‑1 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       >  2:1912630/1‑100 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       >  1:532224/1‑100 (MQ=255)
                                                         agCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGt                                      <  1:615313/100‑1 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  2:1426097/1‑100 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  1:1476363/1‑100 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  1:292476/1‑100 (MQ=255)
                                                            gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                   >  1:719827/1‑100 (MQ=255)
                                                            gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                   <  2:2331717/100‑1 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCGGGGAAACTGTAGAt                                  >  2:1967811/1‑100 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                  >  2:744396/1‑100 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                  <  2:2679219/100‑1 (MQ=255)
                                                              gAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAtt                                 <  2:2744676/100‑1 (MQ=255)
                                                              gAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAtt                                 <  1:1996373/100‑1 (MQ=255)
                                                                 ggggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              >  2:2848487/1‑100 (MQ=255)
                                                                 ggggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              <  1:2768368/100‑1 (MQ=255)
                                                                  gggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTccc                             <  2:2236704/100‑1 (MQ=255)
                                                                  gggAAAATCCATGCTGATTCATCACGGCTTCAACACGCTCTCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              <  2:2646121/99‑1 (MQ=255)
                                                                    gAAAATCCATGCTGGTTCCTCGCTGATTTAACACGCTTGCGGGTGGGCTGGCTCACACCGCTTTCTTTCTTCAGCACCCGGGGAACTGTAGATTTcccc                            <  2:1244669/99‑1 (MQ=255)
                                                                       aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTTGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                        >  2:333648/1‑100 (MQ=255)
                                                                       aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                        >  2:867032/1‑100 (MQ=255)
                                                                        aTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgcc                      <  1:942215/100‑1 (MQ=255)
                                                                         tCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgccg                     >  2:811003/1‑100 (MQ=255)
                                                                           cATGCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCt                   <  1:2614123/100‑1 (MQ=255)
                                                                              gCTGATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTaa                >  1:2069351/1‑100 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg             >  2:2754231/1‑100 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg             <  1:2575107/100‑1 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg             <  2:1972594/100‑1 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAgcg             >  2:2384076/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGat       >  1:2616678/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGat       >  1:1712486/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGat       >  2:2123185/1‑100 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCGCGCTTAAgcggcgtgg       <  1:1098018/100‑7 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGata      >  2:1247080/1‑100 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGata      >  2:1498202/1‑100 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTGGATTTCCCCACGCCGCTTCAGCGCGCGATATc    >  1:1115825/1‑100 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATc    <  1:2500339/100‑1 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATc    >  2:1291735/1‑100 (MQ=255)
                                                                                           tGCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCt   >  1:783662/1‑100 (MQ=255)
                                                                                            gCTTTAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAATGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCtt  >  2:736219/1‑100 (MQ=255)
                                                                                                |                                                                                                 
GACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGGGTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCCGCTTAAGCGCGCGATATCTT  >  NC_000913/4467033‑4467224

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: