Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,285 A→G 16.6% 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,2850AG16.6% 295.5 / 51.7 133intergenic (+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 (47/64);  new base G (14/8);  total (61/72)
Fisher's exact test for biased strand distribution p-value = 9.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 136 total aligned reads displayed.

CTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAT  >  NC_000913/1388187‑1388384
                                                                                                  |                                                                                                   
    gCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                >  2:1076316/1‑100 (MQ=255)
      gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                              >  1:2343923/1‑100 (MQ=255)
      gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAAGATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                              <  1:2608601/100‑1 (MQ=255)
       tttCACATACCGCGATTGCCAATAAGTTGGGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTTAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                             <  1:1742743/100‑1 (MQ=255)
        ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAg                                                                                            <  2:2298046/100‑1 (MQ=39)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:337108/100‑1 (MQ=39)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:2440297/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                           <  2:2037286/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                           >  1:1775982/1‑100 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAggn                                                                                          <  2:1746332/100‑2 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGc                                                                                          >  1:85874/1‑100 (MQ=255)
           acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCa                                                                                         <  2:929192/100‑1 (MQ=255)
           acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCa                                                                                         <  2:1687928/100‑1 (MQ=255)
            caTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                        <  1:907298/100‑1 (MQ=255)
            caTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                        >  1:1137941/1‑100 (MQ=255)
            caTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                        >  2:2829173/1‑100 (MQ=255)
              tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCatat                                                                                      >  2:1300587/1‑100 (MQ=255)
              tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCata                                                                                       >  1:829841/1‑99 (MQ=255)
                ccGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTa                                                                                    <  2:255610/100‑1 (MQ=255)
                  gcgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                             <  1:2531000/89‑1 (MQ=255)
                     aTTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGCATAttcgcgct                                                                               <  2:1854517/100‑7 (MQ=25)
                     aTTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTa                                                                               <  2:2834533/100‑1 (MQ=255)
                       tGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATg                                                                             >  1:1262103/1‑100 (MQ=255)
                        gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTACGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGc                                                                            >  1:2249464/1‑100 (MQ=255)
                         ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGCTGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                           >  1:2724408/1‑100 (MQ=255)
                         ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                           <  2:2006597/100‑1 (MQ=255)
                         ccAATAAGTTGCGGGAATATGGGCTGAGTCAGAAGAAGAACGCAGAGTAAGCGCGAATATGCCTGGTGGTGCAACAGCATTAGGCATATTAAAGTa                                                                               >  1:2392009/1‑96 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                        <  1:337520/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                        <  1:1987030/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                        <  1:2194349/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                        <  2:1697656/100‑1 (MQ=255)
                               aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                     >  1:2725660/1‑100 (MQ=255)
                               aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                     <  2:2308928/100‑1 (MQ=255)
                               aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                     >  1:806830/1‑100 (MQ=255)
                                gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                    <  2:2113739/100‑1 (MQ=255)
                                   gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                 <  2:1265537/100‑1 (MQ=255)
                                    cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                >  1:2075393/1‑100 (MQ=255)
                                     gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                               <  2:1582656/100‑1 (MQ=255)
                                     gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                               <  2:1505175/100‑1 (MQ=255)
                                     gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                               <  2:1212761/100‑1 (MQ=255)
                                      ggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCa                                                              >  2:1884669/1‑100 (MQ=255)
                                        aaTATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                            <  1:1917840/100‑1 (MQ=255)
                                          tatGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTAATGGTGCAACACCATCAGGCATATTAAATTAAGCTTTCAGTACAACCagaac                                                          >  1:1061939/1‑98 (MQ=255)
                                           atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgct                                                         <  2:546560/100‑1 (MQ=255)
                                            tGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGTCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgct                                                         <  1:1805259/99‑1 (MQ=255)
                                            tGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctg                                                        >  2:195661/1‑100 (MQ=255)
                                             ggTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgc                                                       <  2:867419/100‑1 (MQ=255)
                                              gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgct                                                      <  2:2252030/100‑1 (MQ=255)
                                                 tGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                   <  1:733381/100‑1 (MQ=255)
                                                  gAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                  <  2:1827853/100‑1 (MQ=255)
                                                   aGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTa                                                 >  1:205363/1‑100 (MQ=255)
                                                    ggcAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAGCAGGCATATTAAATTATGCTTTCAGTACAGGCAGAGCTGCTTCGt                                                  <  1:2239233/96‑1 (MQ=255)
                                                    gTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                                >  1:1061377/1‑100 (MQ=255)
                                                     tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAAt                                               >  1:200759/1‑100 (MQ=255)
                                                       agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                             >  2:890181/1‑100 (MQ=255)
                                                       agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                             >  1:644455/1‑100 (MQ=255)
                                                        gaagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCg                                            >  1:451441/1‑100 (MQ=255)
                                                         aagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTAAAGCCAGAGCTGCTTCGTAAtac                                             >  1:1920755/1‑96 (MQ=255)
                                                           gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCt                                         >  2:508971/1‑100 (MQ=255)
                                                            aagaaCGAAGAGTAAGCGCGAATATGCCTGATGTTGCAACACCATCATGCATATTAAATTACGCTTTCTGTACAGCCAGAGCCGCTTCGTACTCCGGCTc                                        <  1:354512/100‑1 (MQ=255)
                                                             agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                        >  1:2461372/1‑99 (MQ=255)
                                                             agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCg                                       <  2:2433676/100‑1 (MQ=255)
                                                              gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                      <  2:2779761/100‑1 (MQ=255)
                                                                aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                    >  2:233097/1‑100 (MQ=255)
                                                                 cGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                   >  2:13441/1‑100 (MQ=255)
                                                                  gAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCGGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                  <  1:1562847/100‑1 (MQ=255)
                                                                  gAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                  <  1:1758172/100‑1 (MQ=255)
                                                                   aaGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                  >  2:1340805/1‑99 (MQ=255)
                                                                   aaGAGTAAGCGCGAATATACCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggcg                                 <  2:1969605/100‑3 (MQ=255)
                                                                    agagTAAGCGCGAATATGCCGGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGgtggtg                                 <  1:1471532/99‑1 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATTAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTAggtgtgcggcg                             <  1:1531963/100‑8 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                             <  1:2618425/100‑1 (MQ=255)
                                                                       gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                             >  2:2663677/1‑100 (MQ=255)
                                                                        tAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTc                            <  2:1137943/100‑1 (MQ=255)
                                                                         aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTACATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                           <  1:1867553/100‑1 (MQ=255)
                                                                         aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                           <  2:1581433/100‑1 (MQ=255)
                                                                           gcgcGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                         <  2:800818/100‑1 (MQ=255)
                                                                             gcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcca                       >  2:639058/1‑100 (MQ=255)
                                                                               gAATATGCCGGATGGTGCAACACCATAAGGCATCTTAAATTATGCTTTCAGTACAGCCAGAGCTGCTGCGTAATCCGGCTCGGTGGTGATTTCATccacc                     <  1:28924/100‑1 (MQ=255)
                                                                                ttccgatcTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccacca                    <  1:802212/93‑1 (MQ=255)
                                                                                aaTATGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccacca                    >  2:1751442/1‑100 (MQ=38)
                                                                                   atGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct                 >  2:2665081/1‑100 (MQ=38)
                                                                                   atGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct                 >  2:291288/1‑100 (MQ=38)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg               >  2:2259/1‑100 (MQ=255)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg               <  2:806832/100‑1 (MQ=255)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg               >  2:1438060/1‑100 (MQ=255)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  2:930280/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  2:722492/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  1:1089882/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  2:1573562/1‑100 (MQ=39)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              >  2:24215/1‑100 (MQ=38)
                                                                                      ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc              <  1:121025/100‑1 (MQ=255)
                                                                                      ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgatggc              >  1:1929103/1‑100 (MQ=255)
                                                                                       cTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggct             <  1:594695/100‑1 (MQ=255)
                                                                                        tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggctg            >  1:2108403/1‑100 (MQ=255)
                                                                                          aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa          <  2:1156513/100‑1 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaan        <  2:2258028/100‑2 (MQ=255)
                                                                                              tGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc      <  2:1315118/100‑1 (MQ=255)
                                                                                               gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca     <  2:1176169/100‑1 (MQ=255)
                                                                                                cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcac    >  2:509082/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: