Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,979 T→C 16.1% T253A (ACG→GCG)  yhbE ← EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,331,9790TC16.1% 262.7 / 41.2 117T253A (ACG→GCG) yhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base C (17/2);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 3.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 125 total aligned reads displayed.

CGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAT  >  NC_000913/3331880‑3332076
                                                                                                   |                                                                                                 
   gaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACGGCGTAAACAGTGGGGTGAGCGtgat                                                                                                <  2:2669228/100‑1 (MQ=255)
   gaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                >  2:405745/1‑100 (MQ=255)
   gaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                >  1:2850844/1‑100 (MQ=255)
    aTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                               <  1:797671/100‑1 (MQ=255)
     tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAACTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                              >  1:367734/1‑100 (MQ=255)
      aaCCTAAAAGGTTTAACATCGGTCTGGAGAAGAAATCGGGCCAGGCCAGTGATAAAAAAAAAAAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                <  1:1004856/97‑1 (MQ=255)
      aaCCTAAAAGGTTAAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTgatgat                                                                                             >  2:1961229/1‑100 (MQ=255)
        ccTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcg                                                                                           <  2:2852724/100‑1 (MQ=255)
         cTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                          >  2:2169703/1‑100 (MQ=255)
           aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                        <  2:2428182/100‑1 (MQ=255)
           aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                        >  2:2821864/1‑100 (MQ=255)
           aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                        <  1:497875/100‑1 (MQ=255)
               ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAc                                                                                    <  1:1934506/100‑1 (MQ=255)
               ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAc                                                                                    <  1:1356597/100‑1 (MQ=255)
                gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                   >  2:1731958/1‑100 (MQ=255)
                  ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                 <  2:1413571/100‑1 (MQ=255)
                    aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcg                                                                               >  1:2420336/1‑100 (MQ=255)
                    aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCACTCACCTgcg                                                                               <  1:1176884/100‑1 (MQ=255)
                         cGGTCTGGCGAAGAAATCGGGCCAGGCCAGTTATAAACGATCTGAAAGAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                          <  1:1969932/100‑1 (MQ=255)
                         cGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                          <  2:942748/100‑1 (MQ=255)
                          ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcct                                                                         >  2:2442787/1‑100 (MQ=255)
                               ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                    >  2:2377136/1‑100 (MQ=255)
                               ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                    <  2:120952/100‑1 (MQ=255)
                               ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                    >  1:2292392/1‑100 (MQ=255)
                               ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                    >  1:962877/1‑100 (MQ=255)
                                gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGc                                                                   >  1:1810231/1‑100 (MQ=255)
                                gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGc                                                                   <  1:2888377/100‑1 (MQ=255)
                                gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGc                                                                   <  1:1593945/100‑1 (MQ=255)
                                 cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGCAAAAAACACCGTACACCGTGGGGTGAGCGTGATACTCGCGCTCACCTGCGCTGCCTGCCAGc                                                                   >  1:1449539/1‑99 (MQ=255)
                                 cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCg                                                                  >  2:1788212/1‑100 (MQ=255)
                                  gaagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGa                                                                 >  1:2584130/1‑100 (MQ=255)
                                   aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                <  1:1728716/100‑1 (MQ=255)
                                   aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                >  2:487248/1‑100 (MQ=255)
                                      aaaTAGGGCCAGGCCAGTGATAAAAGTTCTGTAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCTCCTGCGCTGGCAGGCAGCgcagg                                                              <  1:854103/99‑5 (MQ=255)
                                         tCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCg                                                          <  2:1032662/100‑1 (MQ=255)
                                           gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                        >  2:766468/1‑100 (MQ=255)
                                              ccaggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTcc                                                     >  2:1446235/1‑100 (MQ=255)
                                                   ccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAg                                                <  2:842698/100‑1 (MQ=255)
                                                    cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                               >  2:1267931/1‑100 (MQ=255)
                                                    cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                               >  2:2696294/1‑100 (MQ=255)
                                                      gTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGTGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgc                                             <  2:2050308/100‑1 (MQ=255)
                                                       tGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgcg                                            <  1:1639447/100‑1 (MQ=255)
                                                       tGATAAAAGATCTCAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgcg                                            >  2:1982167/1‑100 (MQ=255)
                                                        gATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAACTGCGCTGCCAGACAGCGAGCCATCGCTTCCGCCAg                                                >  1:2906263/1‑95 (MQ=255)
                                                         aTAAAAGGTCTGAAAACAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                          <  1:2802173/100‑1 (MQ=255)
                                                         aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                          <  2:1834832/100‑1 (MQ=255)
                                                         aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                          >  2:1907835/1‑100 (MQ=255)
                                                         aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                          <  1:9964/100‑1 (MQ=255)
                                                         aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                          >  1:1503735/1‑100 (MQ=255)
                                                          tAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                         <  1:1811505/100‑1 (MQ=255)
                                                           aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                        >  2:2243771/1‑100 (MQ=255)
                                                           aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                        >  2:553202/1‑100 (MQ=255)
                                                              aGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATc                                     <  1:2811273/100‑1 (MQ=255)
                                                               gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATcc                                    <  1:1214686/100‑1 (MQ=255)
                                                                  cTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                 <  1:325197/100‑1 (MQ=255)
                                                                   tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                <  2:1567942/100‑1 (MQ=255)
                                                                    gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAGCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                <  1:2370510/99‑1 (MQ=255)
                                                                    gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                               >  2:1899623/1‑100 (MQ=255)
                                                                    gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAGCGCTTCCCCCAGGGCGCCATATCCt                                   >  1:283309/1‑96 (MQ=255)
                                                                     aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                              >  1:628839/1‑100 (MQ=255)
                                                                      aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                             <  1:1483912/100‑1 (MQ=255)
                                                                      aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                             <  2:1290101/100‑1 (MQ=255)
                                                                       aaaaaCAGCGTGAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                            <  1:1405438/100‑1 (MQ=255)
                                                                        aaaaCAGCGTAAACAGTGGGGTGAGCGTGAGGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGt                           <  2:1761234/100‑1 (MQ=255)
                                                                          aaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAt                         >  2:941130/1‑100 (MQ=255)
                                                                          aaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAt                         <  2:1329603/100‑1 (MQ=255)
                                                                           aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAtt                        >  1:219281/1‑100 (MQ=255)
                                                                           aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAtt                        >  1:1851335/1‑100 (MQ=255)
                                                                            cAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                       >  2:2662059/1‑100 (MQ=255)
                                                                            cAGCGTAAACAGTGGGGTCAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTACCGCCAGGGCGCCATATCCTACCAAGGTATTc                       <  2:209207/100‑1 (MQ=255)
                                                                              gCGTAAACAGTTGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                     <  2:127957/100‑1 (MQ=255)
                                                                              gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                     <  2:2584145/100‑1 (MQ=255)
                                                                              gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                     <  1:267161/100‑1 (MQ=255)
                                                                                 tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                  >  1:262421/1‑100 (MQ=255)
                                                                                 tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                  <  1:1893396/100‑1 (MQ=255)
                                                                                 tAAACAGTGGGGTGAGCGTGATGATCGCGCGCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                  <  1:2468073/100‑1 (MQ=255)
                                                                                  aaaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCTCCATATCCTACCAAGGTATTCAGTcc                  <  1:2188989/99‑1 (MQ=255)
                                                                                     cAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa              >  2:144578/1‑100 (MQ=255)
                                                                                     agcgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa              >  1:2082572/7‑100 (MQ=38)
                                                                                      gcgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa             >  2:844806/6‑100 (MQ=38)
                                                                                      gcgcaggTGAGCGCGATCATCACGCTCACCCGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa             >  2:2243623/6‑100 (MQ=37)
                                                                                      aGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa             <  2:2245987/100‑1 (MQ=255)
                                                                                       ggggggTAGGCGTGTCGATCTCGCTCACCTGCGCTGCCGGCTAGCGAGGCATCGACTCCGCGAGGGAGCCACATCCTAGCAAGTTATTCAGTcc                  <  1:2191674/92‑1 (MQ=255)
                                                                                       cgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa            >  1:803190/5‑100 (MQ=38)
                                                                                       cgcaggTGAGCGCGATCATCACGCTCACCCGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa            >  1:166262/5‑100 (MQ=37)
                                                                                        gcaggagAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCGTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCacga            >  1:2248335/7‑96 (MQ=38)
                                                                                        gcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTGTTCAGTCCGCaaaa            >  2:2867605/4‑99 (MQ=38)
                                                                                         ggggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGTCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa          >  1:2648350/1‑100 (MQ=255)
                                                                                         caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa          >  2:211666/3‑100 (MQ=38)
                                                                                         caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa          >  1:1495260/3‑100 (MQ=39)
                                                                                         caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa          >  1:2089989/3‑100 (MQ=38)
                                                                                          aggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt         >  2:613578/2‑100 (MQ=39)
                                                                                          aggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt         >  2:1537353/2‑100 (MQ=39)
                                                                                           ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt        >  2:2492212/1‑100 (MQ=255)
                                                                                           ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt        >  2:1748974/1‑100 (MQ=255)
                                                                                           ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt        >  1:1164598/1‑100 (MQ=255)
                                                                                           ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt         >  2:1698092/1‑99 (MQ=39)
                                                                                              gAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa     <  1:1015042/100‑1 (MQ=255)
                                                                                               aGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc    <  1:1267926/100‑1 (MQ=255)
                                                                                                 cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt  >  2:391172/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAT  >  NC_000913/3331880‑3332076

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