Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,283 C→T 9.6% 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,2830CT 9.6% 305.9 / 14.4 104intergenic (+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 (46/48);  new base T (7/3);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 106 total aligned reads displayed.

TATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCA  >  NC_000913/1388134‑1388432
                                                                                                                                                     |                                                                                                                                                     
tATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAATCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGc                                                                                                                                                       <  2:1205300/150‑1 (MQ=255)
    aaCCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAaca                                                                                                                                                   <  2:284562/150‑1 (MQ=255)
       ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGCACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                >  1:625381/1‑150 (MQ=255)
       ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                >  1:528418/1‑150 (MQ=255)
       ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCa                                                                                                                                                <  2:109044/150‑1 (MQ=255)
        cAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                                                                                                               >  2:163381/1‑150 (MQ=255)
           cTTTATCGCAATAATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATTGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGAGGGTGCAACACCATCAg                                                                                                                                            <  1:1059100/150‑1 (MQ=255)
            tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAACCGTCTCGGCGTTTCACATACCGCGATTGCCAATACGTTGCGGGAATATGGTCGGCGTCAGAAGCAGAACGAAGAGGACGCGCGAATATGCCTGATGGTGTTGCACCATCAgg                                                                                                                                           <  1:1158856/150‑1 (MQ=255)
             ttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCAGGc                                                                                                                                          <  1:894964/150‑1 (MQ=255)
              tATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGACTCGGAGTTTCACATACAGCGATTGCCAATAAGTAGAGGGAATAAGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATAAGCCTGATGGTGCAACACCATCAGGCa                                                                                                                                         >  1:238495/1‑150 (MQ=255)
               aTCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCat                                                                                                                                        >  2:154767/1‑150 (MQ=255)
                 cGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTTCAACACCATCAGGCatat                                                                                                                                      >  2:441343/1‑150 (MQ=255)
                       taacccAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAAtt                                                                                                                                <  2:1598479/147‑1 (MQ=255)
                        aTCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTa                                                                                                                               <  2:528418/150‑1 (MQ=255)
                         tCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTAt                                                                                                                              >  2:434084/1‑150 (MQ=255)
                          cccAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATg                                                                                                                             >  2:1654353/1‑150 (MQ=255)
                           ccAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGc                                                                                                                            <  1:355365/150‑1 (MQ=255)
                               cACGCGCAACCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAAGAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGCGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATAGTAAATTATGCTTTc                                                                                                                        >  1:1163360/1‑150 (MQ=255)
                                aCGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCa                                                                                                                       <  2:1045917/150‑1 (MQ=255)
                                aCGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCa                                                                                                                       <  1:935598/150‑1 (MQ=255)
                                  gcgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                                                                     <  2:802528/150‑1 (MQ=255)
                                   cgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                                                                    >  1:130595/1‑150 (MQ=255)
                                     cAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACa                                                                                                                  <  1:982771/150‑1 (MQ=255)
                                      aaaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAg                                                                                                                 <  2:322009/150‑1 (MQ=255)
                                       aaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                                                                >  2:934003/1‑150 (MQ=255)
                                       aaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                                                                <  1:1428036/150‑1 (MQ=255)
                                       aaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                                                                                                                <  1:1214459/150‑1 (MQ=255)
                                        aCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                                                                               <  2:705135/150‑1 (MQ=255)
                                           ggCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                                                                            >  2:1479477/1‑150 (MQ=255)
                                           ggCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCagag                                                                                                           >  2:439602/1‑150 (MQ=255)
                                                    tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                   >  1:1056981/1‑150 (MQ=255)
                                                    tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                   <  1:957729/150‑1 (MQ=255)
                                                    tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                   >  1:1204731/1‑150 (MQ=255)
                                                    gctacgCGTTTCACATACCTCGATTGCAAATAAGATGAGTTAATAGGTGCTGTGTCAGACGAAGAACGAATATTAATCTCGAATATGCCTGATGGTGCAACACCAGCAGGCATATTAAATTCTGGTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                   <  2:264466/145‑1 (MQ=255)
                                                     ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                                                                  <  1:434084/150‑1 (MQ=255)
                                                      tcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCCTCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCCTa                                                                                                 >  1:70710/1‑150 (MQ=255)
                                                             ttCACATACCGCGATTGAGAATAAGTTGCTGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACCGCCAGAGCTGCTTCGTAATCCGGc                                                                                          <  1:1137378/150‑1 (MQ=255)
                                                              tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGCGCTGCTTCGTAATCCGGCt                                                                                         >  1:739497/1‑150 (MQ=255)
                                                              tCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCATCGTAATCCGGCt                                                                                         >  1:693843/1‑150 (MQ=255)
                                                                 caTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                                                                      >  1:770609/1‑150 (MQ=255)
                                                                   tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtg                                                                                    >  2:40092/1‑150 (MQ=255)
                                                                    aCCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggtgg                                                                                   >  1:49914/1‑150 (MQ=255)
                                                                       gcgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGa                                                                                >  1:1020043/1‑150 (MQ=255)
                                                                         gATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGAtt                                                                              <  1:934003/150‑1 (MQ=255)
                                                                           ttGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTc                                                                            <  2:1411755/150‑1 (MQ=255)
                                                                           ttGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCAGGCTCGGTGGTGAttta                                                                            >  2:915154/1‑149 (MQ=255)
                                                                              ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTTCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGATGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                         >  1:105757/1‑150 (MQ=255)
                                                                              ccAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                         <  2:130595/150‑1 (MQ=255)
                                                                               cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcc                                                                        <  2:1567488/150‑1 (MQ=255)
                                                                               cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACAACATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcc                                                                        <  2:101640/150‑1 (MQ=255)
                                                                                aaTATGTTGCGTGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcca                                                                       <  1:1637281/150‑1 (MQ=255)
                                                                                         cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                              >  2:562631/1‑150 (MQ=255)
                                                                                         cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                              <  1:154767/150‑1 (MQ=255)
                                                                                         cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGCCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCTTAATCCGTCTCGGTGGTCATTTCATCCACCAgctggc                                                              >  1:920109/1‑150 (MQ=255)
                                                                                           ggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggctg                                                            <  1:1363821/150‑1 (MQ=255)
                                                                                                atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATCTGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAAt                                                       <  1:642077/150‑1 (MQ=255)
                                                                                                atGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAAt                                                       >  2:1681517/1‑150 (MQ=255)
                                                                                                  ggTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACGCCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATca                                                     >  2:733605/1‑150 (MQ=255)
                                                                                                   gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATcac                                                    >  2:570735/1‑150 (MQ=255)
                                                                                                     cTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAt                                                  >  2:165935/1‑150 (MQ=255)
                                                                                                        aGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGt                                               >  1:180260/1‑150 (MQ=255)
                                                                                                           cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAt                                            >  1:31140/1‑150 (MQ=255)
                                                                                                           cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAt                                            >  1:776042/1‑150 (MQ=255)
                                                                                                           cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAt                                            >  1:769600/1‑150 (MQ=255)
                                                                                                             gaagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAttt                                          <  2:1639240/150‑1 (MQ=255)
                                                                                                               agaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTc                                        >  1:71970/1‑150 (MQ=255)
                                                                                                                  agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                     <  2:1056981/150‑1 (MQ=255)
                                                                                                                  agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                     <  2:1204731/150‑1 (MQ=255)
                                                                                                                  agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGACAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                     <  2:685359/150‑1 (MQ=255)
                                                                                                                   gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCa                                    <  1:439602/150‑1 (MQ=255)
                                                                                                                             tAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACaa                          >  2:1669794/1‑150 (MQ=255)
                                                                                                                              aaGCGCGAATATGACTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAAc                         <  2:196240/150‑1 (MQ=255)
                                                                                                                               aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACg                        >  1:972136/1‑150 (MQ=255)
                                                                                                                               aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACg                        <  2:1020043/150‑1 (MQ=255)
                                                                                                                                 cgcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCCATAACCACACCGGc                      >  1:1196866/1‑150 (MQ=255)
                                                                                                                                 cgcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGc                      <  1:570735/150‑1 (MQ=255)
                                                                                                                                  tagAATATGCCTGATTGTGCAACACCATAAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTAGTAATCCGGATAGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCa                     <  1:747930/148‑1 (MQ=255)
                                                                                                                                  gcgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCa                     <  2:1258142/150‑1 (MQ=255)
                                                                                                                                   cgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAc                    <  1:1512500/150‑1 (MQ=255)
                                                                                                                                   cgAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAc                    <  1:424034/150‑1 (MQ=255)
                                                                                                                                     aaTATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTCATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAcgc                  >  1:1362895/1‑150 (MQ=255)
                                                                                                                                         tGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCaa                                   >  1:842753/1‑129 (MQ=255)
                                                                                                                                         tGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCaa                                   <  2:842753/129‑1 (MQ=255)
                                                                                                                                         tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTg              >  1:431112/1‑150 (MQ=255)
                                                                                                                                          gCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGc             >  1:966375/1‑150 (MQ=255)
                                                                                                                                          gCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGc             <  1:924277/150‑1 (MQ=255)
                                                                                                                                           ccTGATGGTTCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGcc            <  2:49914/150‑1 (MQ=255)
                                                                                                                                           ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGcc            >  2:93283/1‑150 (MQ=255)
                                                                                                                                           ccTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGAGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGcc            >  1:90941/1‑150 (MQ=255)
                                                                                                                                           ccTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGcc            <  1:1654353/150‑1 (MQ=255)
                                                                                                                                            cTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa           >  2:1024845/1‑150 (MQ=255)
                                                                                                                                            cTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATATTCGTCAATAACCACACCGGCACGCGCTGccc           >  1:1062126/1‑149 (MQ=255)
                                                                                                                                            cTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCGGGCTGAAAATCACATTGTCATTTTCGTCAATAACCCCAACGGCACGCGCTGCCa           >  2:400213/1‑150 (MQ=255)
                                                                                                                                             tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAg          >  1:1335730/1‑150 (MQ=255)
                                                                                                                                                tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAcc       <  1:1479477/150‑1 (MQ=255)
                                                                                                                                                tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAcc       <  1:40092/150‑1 (MQ=255)
                                                                                                                                                  gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCtt     <  1:1710153/150‑1 (MQ=255)
                                                                                                                                                    gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTc   <  1:1681517/150‑1 (MQ=255)
                                                                                                                                                    gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTc   <  2:693843/150‑1 (MQ=255)
                                                                                                                                                     cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTTAATAACCACAACGGCACGCGCTGCCAGACCTTTCa  >  2:1074248/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCA  >  NC_000913/1388134‑1388432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: