Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,283 C→T 16.4% intergenic (+22/+22) 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,2830CT16.4% 304.5 / 49.5 134intergenic (+22/+22)tyrR/tpxaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base C (48/64);  new base T (14/8);  total (62/72)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-02

Only 100 of 136 total aligned reads displayed.

GTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCAC  >  NC_000913/1388185‑1388382
                                                                                                  |                                                                                                   
      gCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                              >  2:1076316/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                                                                                         <  1:2688170/100‑1 (MQ=38)
           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:1137941/1‑100 (MQ=255)
              caTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                      <  1:907298/100‑1 (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)
                       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                                                                      <  2:1697656/100‑1 (MQ=255)
                              aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  1:2194349/100‑1 (MQ=255)
                              aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  1:1987030/100‑1 (MQ=255)
                              aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  1:337520/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:1212761/100‑1 (MQ=255)
                                       gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                             <  2:1505175/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:2472723/1‑100 (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)
                                                                                                  |                                                                                                   
GTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCAC  >  NC_000913/1388185‑1388382

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: