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.

GACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATCTT  >  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)
           atgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAACATCCATGCTGATTCATCACTGCTTTAACACGCTCGCgg‑g                                                                                    >  2:2639194/1‑100 (MQ=255)
            tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑Gt                                                                                   >  1:1343168/1‑100 (MQ=255)
            tgatgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑Gt                                                                                   >  2:1239836/1‑100 (MQ=255)
              atgGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTcgcgagatc                                                                                  >  1:913097/1‑95 (MQ=255)
                   gACCAC‑TTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGc                                                                           <  2:1884009/100‑1 (MQ=255)
                    aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGc                                                                           <  1:1689975/100‑1 (MQ=255)
                    aCCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGc                                                                           <  1:42320/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCt                                                                          <  1:1507292/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCt                                                                          <  1:1232248/100‑1 (MQ=255)
                     ccACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCt                                                                          <  2:515407/100‑1 (MQ=255)
                      cACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTc                                                                         <  1:848590/100‑1 (MQ=255)
                        cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAc                                                                       >  2:2713408/1‑100 (MQ=255)
                        cTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAc                                                                       >  1:999590/1‑100 (MQ=255)
                          ttATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgc                                                                     <  1:742043/100‑1 (MQ=255)
                          ttATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgc                                                                     <  2:1152253/100‑1 (MQ=255)
                           tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgcc                                                                    >  2:2768347/1‑100 (MQ=255)
                           tATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCAcgcc                                                                    <  1:1493157/100‑1 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCtt                                                                >  2:1141557/1‑100 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCtt                                                                >  1:1504914/1‑100 (MQ=255)
                               gCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCtt                                                                >  1:792708/1‑100 (MQ=255)
                                 tttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTc                                                              <  1:1899429/100‑1 (MQ=255)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgcgagcg                                                              <  2:959148/100‑10 (MQ=37)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGGCGTGAGCCAGCTCAcccgcgagcg                                                              <  1:1272463/100‑10 (MQ=37)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCg                                                             <  1:344470/100‑1 (MQ=255)
                                  ttGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGAAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCg                                                             <  1:20820/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGt                                                            <  2:2455653/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGt                                                            <  2:1527060/100‑1 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGt                                                            >  2:37546/1‑100 (MQ=255)
                                   tGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGt                                                            >  2:2681823/1‑100 (MQ=255)
                                    gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGtt                                                           >  2:2724365/1‑100 (MQ=255)
                                    gCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGtt                                                           >  2:582896/1‑100 (MQ=255)
                                     cccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTa                                                          >  2:2575089/1‑100 (MQ=255)
                                      ccACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTAt                                                         <  2:2870082/100‑1 (MQ=255)
                                       cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                        >  1:27511/1‑100 (MQ=255)
                                       cACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTAtt                                                        >  2:2835907/1‑100 (MQ=255)
                                        tctcATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTc                                                       <  2:1437794/97‑1 (MQ=255)
                                         cGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCa                                                      <  2:2469958/100‑1 (MQ=255)
                                           cATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAga                                                    >  1:1437800/1‑99 (MQ=255)
                                           cATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGc                                                    <  1:1995602/100‑1 (MQ=255)
                                              aGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTAGCTCACGCCGCTTTCGTTATTCAGCAcc                                                 >  1:2515845/1‑100 (MQ=255)
                                               gcgcgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCAccc                                                >  1:2379865/1‑100 (MQ=255)
                                                  cgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCggg                                             >  2:653503/1‑100 (MQ=255)
                                                  cgcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCggg                                             <  2:362973/100‑1 (MQ=255)
                                                   gcgcAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTTACGCCGCTTTCGTTATTCAGCACCCGGGa                                            <  2:2364205/100‑1 (MQ=255)
                                                      cAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAc                                         >  1:676277/1‑100 (MQ=255)
                                                       agagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACt                                        >  2:390117/1‑100 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       <  2:262098/100‑1 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       >  2:1912630/1‑100 (MQ=255)
                                                        gagCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTg                                       >  1:532224/1‑100 (MQ=255)
                                                         agCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGt                                      <  1:615313/100‑1 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  2:1426097/1‑100 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  1:1476363/1‑100 (MQ=255)
                                                          gCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTa                                     >  1:292476/1‑100 (MQ=255)
                                                            gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                   >  1:719827/1‑100 (MQ=255)
                                                            gggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGa                                   <  2:2331717/100‑1 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCGGGGAAACTGTAGAt                                  >  2:1967811/1‑100 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                  >  2:744396/1‑100 (MQ=255)
                                                             ggAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAt                                  <  2:2679219/100‑1 (MQ=255)
                                                              gAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAtt                                 <  2:2744676/100‑1 (MQ=255)
                                                              gAAGGGGAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGAtt                                 <  1:1996373/100‑1 (MQ=255)
                                                                 ggggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              >  2:2848487/1‑100 (MQ=255)
                                                                 ggggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              <  1:2768368/100‑1 (MQ=255)
                                                                  gggAAAATCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTccc                             <  2:2236704/100‑1 (MQ=255)
                                                                  gggAAAATCCATGCTGATTCATCACGGCTTCAACACGCTCTCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTcc                              <  2:2646121/99‑1 (MQ=255)
                                                                    gAAAATCCATGCTGGTTCCTCGCTGATTTAACACGCTTGCGG‑GTGGGCTGGCTCACACCGCTTTCTTTCTTCAGCACCCGGGGAACTGTAGATTTcccc                            <  2:1244669/99‑1 (MQ=255)
                                                                       aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTTGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                        >  2:333648/1‑100 (MQ=255)
                                                                       aaTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc                        >  2:867032/1‑100 (MQ=255)
                                                                        aTCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc‑c                      <  1:942215/100‑1 (MQ=255)
                                                                         tCCATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCAcgc‑cg                     >  2:811003/1‑100 (MQ=255)
                                                                           cATGCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCt                   <  1:2614123/100‑1 (MQ=255)
                                                                              gCTGATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTaa                >  1:2069351/1‑100 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAgcg             >  2:2754231/1‑100 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAgcg             <  1:2575107/100‑1 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAgcg             <  2:1972594/100‑1 (MQ=255)
                                                                                 gATTCATCACTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAgcg             >  2:2384076/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGat       >  1:2616678/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGat       >  1:1712486/1‑100 (MQ=255)
                                                                                       tcaCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGat       >  2:2123185/1‑100 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGCGCGCTTAAgcggcgtgg       <  1:1098018/100‑7 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGata      >  2:1247080/1‑100 (MQ=255)
                                                                                        caCTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGata      >  2:1498202/1‑100 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTGGATTTCCCCACGC‑CGCTTCAGCGCGCGATATc    >  1:1115825/1‑100 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATc    <  1:2500339/100‑1 (MQ=255)
                                                                                          cTGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATc    >  2:1291735/1‑100 (MQ=255)
                                                                                           tGCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATCt   >  1:783662/1‑100 (MQ=255)
                                                                                            gCTTTAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAAATGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATCtt  >  2:736219/1‑100 (MQ=255)
                                                                                                |                                                                                                 
GACGGGTAACAATGATGGCGACCACTTTATCGCTTTGCCCACGCATAGCGCGCGCAGAGCGGGAAGGGGAAAATCCATGCTGATTCATCACTGCTTCAACACGCTCGCGG‑GTGAGCTGGCTCACGCCGCTTTCGTTATTCAGCACCCGGGAAACTGTAGATTTCCCCACGC‑CGCTTAAGCGCGCGATATCTT  >  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: ‑