Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,285 A→G 15.3% intergenic (+24/+20) tyrR → / ← tpx aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,388,2850AG15.3% 323.8 / 44.4 137intergenic (+24/+20)tyrR/tpxaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base A (53/63);  new base G (8/13);  total (61/76)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

CTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAT  >  NC_000913/1388187‑1388384
                                                                                                  |                                                                                                   
ctcGGCGTTTCACATACCGCGCTTGCCAATACGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGAGAATATGCCTGATGGTGCAac                                                                                                    >  1:1697318/1‑100 (MQ=255)
ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                    <  1:1595923/100‑1 (MQ=255)
ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                    >  1:2566337/1‑100 (MQ=255)
 tcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAaca                                                                                                   <  1:920855/100‑1 (MQ=255)
 tcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAaca                                                                                                   <  2:1721909/100‑1 (MQ=255)
   ggCGTTTAACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                  <  1:2702147/99‑1 (MQ=255)
      gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                              <  2:1532576/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGCATATGCCTCATGGTGTTGCACCATCAgg                                                                                           <  1:3772853/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:1107274/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:1325559/100‑1 (MQ=39)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:2487199/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:2776041/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                           <  1:1072105/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                           <  1:2035007/100‑1 (MQ=255)
          nacaTACCACGATTGCCAATAAGATGCGGGAAGATGGTATGAGTCAGAAGAAAAGCGAAGAGAAAGCGCGAATATGTCTTATTGTGTAACACCAtaaagc                                                                                          >  1:773223/2‑95 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  1:2263435/99‑1 (MQ=38)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                          <  2:570607/100‑1 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGc                                                                                          <  1:827990/100‑1 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGc                                                                                          <  2:823284/100‑1 (MQ=255)
           acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  1:609966/98‑1 (MQ=38)
           acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGCa                                                                                         <  2:3319780/100‑1 (MQ=38)
                ccGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTa                                                                                    >  2:810683/1‑100 (MQ=255)
                     aTTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTa                                                                               >  2:1945755/1‑100 (MQ=255)
                        gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGc                                                                            >  1:1689408/1‑100 (MQ=255)
                           aaTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCttt                                                                         <  2:664477/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                        <  2:359480/100‑1 (MQ=255)
                             tAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCa                                                                       >  1:1298721/1‑100 (MQ=255)
                                gTTGCGGGCATATGGGCTAAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGCGCCACACCAGCAGGCCTATTCACTAATCCCTTCAGt                                                                     >  1:3041172/1‑99 (MQ=255)
                                gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                    <  2:499025/100‑1 (MQ=255)
                                gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                    >  1:1868744/1‑100 (MQ=255)
                                 ctgagggAATATGGTCTGACTCAAAAGAATAACCACGAGTAAGCGCAAATATGCTAGATGGTGCAACACCATCATTCATACTAAATGCTGCGTTc                                                                        <  1:1639823/91‑1 (MQ=255)
                                  tGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACa                                                                  >  1:2290604/1‑100 (MQ=255)
                                   gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                 <  1:2769118/100‑1 (MQ=255)
                                   gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                 <  1:1008571/100‑1 (MQ=255)
                                    cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                <  2:1790116/100‑1 (MQ=255)
                                    cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                >  2:1713157/1‑100 (MQ=255)
                                    cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCGTTCAGTACAGc                                                                >  2:3432811/1‑100 (MQ=255)
                                    cGGGAATATGGTCTGAGGCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                <  1:1733388/100‑1 (MQ=255)
                                     gggAATATGGTCTGAGTCAGAAGAAGAGCGAAGAGTAAGCGCTAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                               <  1:1916793/100‑1 (MQ=255)
                                        aaTATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                            >  2:3187929/1‑100 (MQ=255)
                                         atatGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCagag                                                           <  2:1755810/100‑1 (MQ=255)
                                          tatGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgc                                                          >  2:3662191/1‑100 (MQ=255)
                                           atGGTGTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgct                                                         <  1:2969902/100‑1 (MQ=255)
                                           atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgct                                                         >  1:1034834/1‑100 (MQ=255)
                                              gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCCACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgct                                                      <  1:3295737/100‑1 (MQ=255)
                                                 tGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                   <  2:1283764/100‑1 (MQ=255)
                                                 tGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                   <  2:1051360/100‑1 (MQ=255)
                                                  gAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                  <  1:2466990/100‑1 (MQ=255)
                                                    gTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                                >  2:2641787/1‑100 (MQ=255)
                                                      cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATc                                              >  1:3032059/1‑100 (MQ=255)
                                                      cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATc                                              >  2:1103506/1‑100 (MQ=255)
                                                      cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCCTCGTAATc                                              >  2:2919118/1‑100 (MQ=255)
                                                       agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                             <  2:1372062/100‑1 (MQ=255)
                                                         aagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCgg                                           >  1:253198/1‑100 (MQ=255)
                                                          agaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGc                                          <  2:2005131/100‑1 (MQ=255)
                                                           gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCt                                         >  2:1007595/1‑100 (MQ=255)
                                                            aagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                        >  2:2213247/1‑100 (MQ=255)
                                                             agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCCGAGCTGCTTCGTAATCCGGCTCg                                       >  1:835241/1‑100 (MQ=255)
                                                              gaaggaagAGTAAGCGCGAATATGCCTGATGGTGCAACACCTTCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                      <  1:1575372/100‑1 (MQ=255)
                                                              gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                      <  2:684588/100‑1 (MQ=255)
                                                                aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                    >  2:116322/1‑100 (MQ=255)
                                                                aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                    >  1:3662221/1‑100 (MQ=255)
                                                                 cGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                   <  1:1282677/100‑1 (MQ=255)
                                                                 cGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                   >  2:2942355/1‑100 (MQ=255)
                                                                 cGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                   <  2:1883211/100‑1 (MQ=255)
                                                                    agagTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGa                                >  1:3499498/1‑100 (MQ=255)
                                                                     gagTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAt                               <  2:2434210/100‑1 (MQ=255)
                                                                      agTAAGCGCGAATATGCCTGATGGTGAAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAACCCGGCTCGGTGGTGAtt                              <  1:3333080/100‑1 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                             <  2:3396389/100‑1 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                             <  1:2067406/100‑1 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                             >  2:1659844/1‑100 (MQ=255)
                                                                         aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                           <  1:3764470/100‑1 (MQ=255)
                                                                          aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCAt                          <  2:3483966/100‑1 (MQ=255)
                                                                           gcgcGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                           <  1:3384822/98‑1 (MQ=255)
                                                                           gcgcGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                         <  2:238007/100‑1 (MQ=255)
                                                                            cgcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcc                        >  2:3058364/1‑100 (MQ=255)
                                                                             gcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcca                       >  2:409208/1‑100 (MQ=255)
                                                                              cgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccac                      <  1:3728695/100‑1 (MQ=255)
                                                                                aaTATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccacca                    <  2:1240677/100‑1 (MQ=255)
                                                                                    tGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct                 >  1:2202579/1‑99 (MQ=38)
                                                                                    tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGCGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctg                >  1:1529300/1‑100 (MQ=255)
                                                                                    tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctg                >  1:2630948/1‑100 (MQ=255)
                                                                                    tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctg                <  2:1838024/100‑1 (MQ=255)
                                                                                     gCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCAGCTCGGTGGTGATTTCATCCACCAgctgg               >  2:1352682/1‑100 (MQ=38)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGGTTTCATCCACCAgctgg               >  2:3025671/1‑100 (MQ=255)
                                                                                     gCCTGATGGGGCAACACCATCAGGCATATTAAATTAGGCTTTCAGTACATCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTGCATCCACCAgctgg               <  1:889788/100‑1 (MQ=255)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  1:2966098/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  1:359640/1‑100 (MQ=39)
                                                                                      ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              <  2:2115751/100‑1 (MQ=255)
                                                                                       cTGAGGGTGCAACACCATCAGGCATATTAAATTAGGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggct             <  1:442575/100‑1 (MQ=255)
                                                                                        tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggctg            >  1:790220/1‑100 (MQ=255)
                                                                                         gATGGTGCAACGCGTTTAGGCATATTAAATTATGCTTTGAGTACAGCCTGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCCCGAgctggc              <  1:1897526/97‑1 (MQ=255)
                                                                                         gATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGa           >  2:108351/1‑100 (MQ=255)
                                                                                          aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa          >  2:2987001/1‑100 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa        >  1:265347/1‑100 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa        <  2:192285/100‑1 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa        <  2:1525092/100‑1 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa        >  1:637531/1‑100 (MQ=255)
                                                                                             gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAAt       <  2:2368615/100‑1 (MQ=255)
                                                                                                  acacCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAt  >  2:585514/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAT  >  NC_000913/1388187‑1388384

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