Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,284 A→T 16.3% intergenic (+23/+21) 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,2840AT16.3% 345.4 / 58.7 135intergenic (+23/+21)tyrR/tpxaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base A (49/64);  new base T (14/8);  total (63/72)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

Only 100 of 135 total aligned reads displayed.

CTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCAC  >  NC_000913/1388187‑1388382
                                                                                                 |                                                                                                  
ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                  >  2:1156531/1‑100 (MQ=255)
  cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                >  2:1582396/1‑100 (MQ=255)
    gCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                              >  2:1076316/1‑100 (MQ=255)
      gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                            >  1:2343923/1‑100 (MQ=255)
       tttCACATACCGCGATTGCCAATAAGTTGGGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTTAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                           <  1:1742743/100‑1 (MQ=255)
        ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAg                                                                                          <  2:2298046/100‑1 (MQ=39)
        ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAg                                                                                          <  2:2668371/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                         <  2:337108/100‑1 (MQ=39)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                         <  1:2688170/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                         <  2:2037286/100‑1 (MQ=38)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                         <  2:1224642/100‑1 (MQ=255)
         tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                         >  1:1775982/1‑100 (MQ=255)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                        <  1:950024/100‑1 (MQ=38)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                        <  1:1930284/100‑1 (MQ=38)
          cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAggn                                                                                        <  2:1746332/100‑2 (MQ=255)
           acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCa                                                                                       <  2:929192/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)
                       tGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATg                                                                           >  1:1262103/1‑100 (MQ=255)
                        gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTACGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGc                                                                          >  1:2249464/1‑100 (MQ=255)
                         ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGCTGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                         >  1:2724408/1‑100 (MQ=255)
                         ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                         >  1:719383/1‑100 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  1:2194349/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  1:337520/100‑1 (MQ=255)
                            aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                      <  2:1697656/100‑1 (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:1505175/100‑1 (MQ=255)
                                     gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                             <  2:1582656/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)
                                                 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)
                                                    gTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                              >  2:1758178/1‑100 (MQ=255)
                                                     tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAAt                                             >  1:200759/1‑100 (MQ=255)
                                                       agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                           >  2:890181/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)
                                                            aagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                      <  1:1884661/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                                <  2:284132/100‑1 (MQ=255)
                                                                  gAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                <  1:811033/100‑1 (MQ=255)
                                                                  gAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                <  1:1758172/100‑1 (MQ=255)
                                                                   aaGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtggt                                >  2:1340805/1‑99 (MQ=255)
                                                                   aaGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGgtggtg                               >  1:511991/1‑100 (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)
                                                                         aaGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTACATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                         <  1:1867553/100‑1 (MQ=255)
                                                                          aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCAt                        >  2:1331340/1‑100 (MQ=255)
                                                                           gcgcGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                       <  2:800818/100‑1 (MQ=255)
                                                                             tttaaTATGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcca                     >  2:2523367/4‑100 (MQ=37)
                                                                             gcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcca                     >  2:639058/1‑100 (MQ=255)
                                                                               taaTATGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccacc                   >  2:1468996/2‑100 (MQ=37)
                                                                               gAATATGCCGGATGGTGCAACACCATAAGGCATCTTAAATTATGCTTTCAGTACAGCCAGAGCTGCTGCGTAATCCGGCTCGGTGGTGATTTCATccacc                   <  1:28924/100‑1 (MQ=255)
                                                                                   atGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct               >  2:291288/1‑100 (MQ=38)
                                                                                     gCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg             >  2:2006707/1‑100 (MQ=38)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg             <  2:806832/100‑1 (MQ=255)
                                                                                     gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg             >  2:2259/1‑100 (MQ=255)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  1:1089882/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  2:722492/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  2:930280/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  2:24215/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  2:1701222/1‑100 (MQ=38)
                                                                                      ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc            >  2:1573562/1‑100 (MQ=39)
                                                                                      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)
                                                                                        tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggctg          >  1:2472723/1‑100 (MQ=255)
                                                                                         gATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGa         >  2:2257512/1‑100 (MQ=255)
                                                                                          aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa        <  2:192012/100‑1 (MQ=255)
                                                                                          aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa        <  2:1156513/100‑1 (MQ=255)
                                                                                            ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaan      <  2:2258028/100‑2 (MQ=255)
                                                                                              tGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAACCAGAGCTGCTTCGTAATCCGACTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc    <  1:306330/100‑1 (MQ=255)
                                                                                               gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca   <  2:1176169/100‑1 (MQ=255)
                                                                                                cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcac  >  2:509082/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCAC  >  NC_000913/1388187‑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: