Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,285 A→G 17.7% 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,2850AG17.7% 284.1 / 56.6 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 (50/62);  new base G (13/11);  total (64/73)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 147 total aligned reads displayed.

TCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACA  >  NC_000913/1388186‑1388383
                                                                                                   |                                                                                                  
tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                    >  2:1979476/1‑100 (MQ=255)
 ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                   >  1:3355074/1‑100 (MQ=255)
 ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAac                                                                                                   >  2:3988077/1‑100 (MQ=255)
   cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                 >  2:1255478/1‑100 (MQ=255)
   cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                 >  1:1680578/1‑100 (MQ=255)
   cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAAAAGAAAGAAGAGTAAGCGCGAATATGCCTGATGGTGAAacac                                                                                                 >  1:740626/1‑100 (MQ=255)
   cGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGGCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                 <  1:782145/100‑1 (MQ=255)
    ggCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACAcc                                                                                                <  2:1765410/100‑1 (MQ=255)
       gTTTCACATACCGCGATTGCCACTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                             <  1:266414/100‑1 (MQ=255)
       gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                                                             <  1:524585/100‑1 (MQ=255)
         ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAg                                                                                           <  1:266672/100‑1 (MQ=39)
         ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAACAg                                                                                           <  1:3055415/100‑1 (MQ=255)
          tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                          <  1:1641461/100‑1 (MQ=39)
          tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                          <  1:3761741/100‑1 (MQ=38)
          tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                          <  2:703269/100‑1 (MQ=255)
          tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                          <  1:1369050/100‑1 (MQ=39)
          tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                          <  2:1049050/100‑1 (MQ=255)
           cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                         <  2:2372907/100‑1 (MQ=38)
           cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAgg                                                                                          <  2:1266118/99‑1 (MQ=255)
           cacaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGc                                                                                         <  2:2834031/100‑1 (MQ=255)
            acaTACCGCGCTTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                         <  1:3692169/99‑1 (MQ=38)
            acaTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTACAACACCATCAGGCa                                                                                        <  1:4096626/100‑1 (MQ=255)
            acaTACCGCGATTGCCAATAAGGTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                          <  1:2172810/98‑1 (MQ=38)
               tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCatat                                                                                     >  1:1341084/1‑100 (MQ=255)
                   gcgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTaaa                                                                                 >  2:2496550/1‑100 (MQ=255)
                    cgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAAt                                                                                <  2:1056028/100‑1 (MQ=255)
                    cgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAAt                                                                                <  2:429126/100‑1 (MQ=255)
                     gATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAAAGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAAtt                                                                               <  1:568536/100‑1 (MQ=255)
                       ttGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTAt                                                                             <  1:1261161/100‑1 (MQ=255)
                        tGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTAt                                                                             <  1:806491/99‑1 (MQ=255)
                        tGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATg                                                                            >  2:847488/1‑100 (MQ=255)
                            aaTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCt                                                                          <  1:2858936/98‑1 (MQ=255)
                             agaaGGTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                       <  1:3098796/98‑1 (MQ=255)
                             aTAAGTTGCTGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGGGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                       <  2:366275/100‑1 (MQ=255)
                             aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGCTGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                       <  1:4338208/100‑1 (MQ=255)
                             aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                       <  2:3376196/100‑1 (MQ=255)
                             aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTc                                                                       <  1:1787563/100‑1 (MQ=255)
                                aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTTATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                    <  1:361474/100‑1 (MQ=255)
                                aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                    <  1:953086/100‑1 (MQ=255)
                                aGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                    >  2:3939184/1‑100 (MQ=255)
                                 gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                   <  1:4014965/100‑1 (MQ=255)
                                   tGCTGGAATATGGTCTGAGTCAGAAGAATAATTAAGAGTAATCGCCAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                   <  1:3009518/98‑1 (MQ=38)
                                   tGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACa                                                                 >  2:1652078/1‑100 (MQ=255)
                                    gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                <  1:3506628/100‑1 (MQ=255)
                                    gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                >  1:4396957/1‑100 (MQ=255)
                                    gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                >  1:4081306/1‑100 (MQ=255)
                                    gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                >  1:396085/1‑100 (MQ=255)
                                     cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACTATCAGGCATATTAAATTATGCTTTGAGTACAg                                                                >  1:1534396/1‑99 (MQ=255)
                                      gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                              >  1:1879699/1‑100 (MQ=255)
                                      gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                              <  2:3041973/100‑1 (MQ=255)
                                       ggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTACAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCa                                                             >  2:3877819/1‑100 (MQ=255)
                                        gAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCag                                                            >  2:990587/1‑100 (MQ=255)
                                         aaTATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                           >  1:3810682/1‑100 (MQ=255)
                                           tatGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgc                                                         <  2:1827653/100‑1 (MQ=255)
                                            atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgct                                                        <  2:4287865/100‑1 (MQ=255)
                                              ggTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgc                                                      >  2:101099/1‑100 (MQ=255)
                                               gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgct                                                     <  1:439029/100‑1 (MQ=255)
                                                   gAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                 <  2:1486934/100‑1 (MQ=255)
                                                    aGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTa                                                <  1:620984/100‑1 (MQ=255)
                                                        agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTTCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                            <  1:4274552/100‑1 (MQ=255)
                                                        agaagaagaACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                            <  1:1801763/100‑1 (MQ=255)
                                                           agaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGc                                         >  1:654743/1‑100 (MQ=255)
                                                                 aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                   >  2:1840608/1‑100 (MQ=255)
                                                                 aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                   >  2:1523721/1‑100 (MQ=255)
                                                                 aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATCATGCTTTCACTACAGCCAGAGCTGCTTCCTAATCCGTCTCggtg                                   >  1:1273649/1‑100 (MQ=39)
                                                                  cGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                  <  2:1303807/100‑1 (MQ=255)
                                                                  cGAAGAGTAAGAGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                  >  2:648909/1‑100 (MQ=255)
                                                                        gTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAttt                            >  2:1609889/1‑100 (MQ=255)
                                                                           aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCan                         <  1:1423321/100‑2 (MQ=255)
                                                                            gcgcGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                        <  2:684494/100‑1 (MQ=255)
                                                                             cgcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcc                       <  1:1840602/100‑1 (MQ=255)
                                                                                 aaTATGCTTGATGGTGCAACAGCCTCAAGCATATTAAATTTTGCTTTCTGTACAGCCAGAGCTGCTTCGTAATCCTGCTCGGTGGTGATTTCATcc                       <  1:3696701/96‑1 (MQ=255)
                                                                                 aaTATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATccacca                   <  1:891752/100‑1 (MQ=255)
                                                                                   tatGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgc                 <  2:3986095/100‑1 (MQ=255)
                                                                                    atGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct                >  2:2767728/1‑100 (MQ=38)
                                                                                    atGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgct                >  2:136577/1‑100 (MQ=255)
                                                                                     tGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctg               >  1:2178198/1‑100 (MQ=38)
                                                                                      gCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg              >  2:3214706/1‑100 (MQ=38)
                                                                                      gCCTGATGGTGTTACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg              >  2:2964983/1‑100 (MQ=255)
                                                                                      gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg              >  2:2958508/1‑100 (MQ=255)
                                                                                      gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg              >  2:343541/1‑100 (MQ=255)
                                                                                      gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCAGCCACCAgctgg              >  1:1460157/1‑100 (MQ=255)
                                                                                      gCCTGATGGGGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg              >  2:3267006/1‑100 (MQ=39)
                                                                                       ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCTGGCTCGGTGGAGATTTCATCCACCAgat                >  1:3966046/1‑95 (MQ=37)
                                                                                       ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc             >  1:3248011/1‑100 (MQ=38)
                                                                                       ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc             >  2:4388785/1‑100 (MQ=38)
                                                                                       ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc             >  2:3378798/1‑100 (MQ=38)
                                                                                       ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc             >  2:3020577/1‑100 (MQ=38)
                                                                                       ccTGATGGTGCATCACCACCAGGCATATTAAATTATGCTTGCAGTTCAGCCAGAGCTGCTTCGTAATCCAGCTCGGAAGTGAATTCATCCACCAactggc             >  1:1949073/1‑100 (MQ=255)
                                                                                       ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc             <  1:782579/100‑1 (MQ=255)
                                                                                              gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAAt      >  2:1999816/1‑100 (MQ=255)
                                                                                               tGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc     >  2:2343869/1‑100 (MQ=255)
                                                                                               tGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc     >  2:3672660/1‑100 (MQ=255)
                                                                                               tGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATc     >  1:2450935/1‑100 (MQ=255)
                                                                                                gCAACACCATCATGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCCTAATCCGGCTCGGTGGTGATTTGATTCTCCCTCTGGCTGAAACTc     >  1:4017342/1‑99 (MQ=255)
                                                                                                gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca    >  1:1980758/1‑100 (MQ=255)
                                                                                                gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca    <  1:2373720/100‑1 (MQ=255)
                                                                                                gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca    >  1:2476676/1‑100 (MQ=255)
                                                                                                 cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcac   >  1:3840894/1‑100 (MQ=255)
                                                                                                  aaCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcaca  <  2:2582728/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACA  >  NC_000913/1388186‑1388383

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: