Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,285 A→G 11.5% 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,2850AG11.5% 256.1 / 23.3 104intergenic (+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 (46/46);  new base G (8/4);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 107 total aligned reads displayed.

AACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGT  >  NC_000913/1388138‑1388434
                                                                                                                                                   |                                                                                                                                                     
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)
                    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)
                             cgcgCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGAAGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAg                                                                                                                        >  2:1363821/1‑150 (MQ=255)
                              gcgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGt                                                                                                                       <  2:802528/150‑1 (MQ=255)
                               cgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                                                                      >  1:130595/1‑150 (MQ=255)
                               cgcAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                                                                                                                      <  1:926172/150‑1 (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:1214459/150‑1 (MQ=255)
                                    aCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGcc                                                                                                                 <  2:705135/150‑1 (MQ=255)
                                       ggCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCaga                                                                                                              >  2:1479477/1‑150 (MQ=255)
                                                tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                     <  1:957729/150‑1 (MQ=255)
                                                tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                     >  1:1056981/1‑150 (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)
                                                                  cgcgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGgtggtg                                                                                   >  2:424034/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)
                                                                            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:776042/1‑150 (MQ=255)
                                                                                                       cAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAt                                              >  1:31140/1‑150 (MQ=255)
                                                                                                         gaagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAttt                                            <  2:1639240/150‑1 (MQ=255)
                                                                                                           agaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTc                                          >  1:71970/1‑150 (MQ=255)
                                                                                                              agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                       <  2:1204731/150‑1 (MQ=255)
                                                                                                              agaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTc                                       <  2:1056981/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                          <  2:1020043/150‑1 (MQ=255)
                                                                                                                           aGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACg                          >  1:972136/1‑150 (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                                     <  2:842753/129‑1 (MQ=255)
                                                                                                                                     tGCCTGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCaa                                     >  1:842753/1‑129 (MQ=255)
                                                                                                                                     tGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTg                >  1:557436/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)
                                                                                                                                        cTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa             >  1:1485985/1‑150 (MQ=255)
                                                                                                                                         tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAg            >  1:1335730/1‑150 (MQ=255)
                                                                                                                                            tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAcc         <  1:40092/150‑1 (MQ=255)
                                                                                                                                            tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAcc         <  1:1479477/150‑1 (MQ=255)
                                                                                                                                              gTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCtt       <  1:1710153/150‑1 (MQ=255)
                                                                                                                                                gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTc     <  2:693843/150‑1 (MQ=255)
                                                                                                                                                gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTc     <  1:1681517/150‑1 (MQ=255)
                                                                                                                                                 cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTTAATAACCACAACGGCACGCGCTGCCAGACCTTTCa    >  2:1074248/1‑150 (MQ=255)
                                                                                                                                                  aaCACCATCAGGCATATTAAATTATGCTATCAGTACACCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCa  <  1:655100/150‑2 (MQ=255)
                                                                                                                                                   acacCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGCGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGt  >  1:902338/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGT  >  NC_000913/1388138‑1388434

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