Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 120,158:1 +A 7.1% intergenic (+23/+20) ampE → / ← aroP ampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913120,1581.A 7.1% 480.8 / 30.9 127intergenic (+23/+20)ampE/aroPampicillin resistance inner membrane protein; putative signaling protein in beta‑lactamase regulation/aromatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (50/68);  new base A (4/5);  total (54/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 133 total aligned reads displayed.

TCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACC  >  NC_000913/120009‑120307
                                                                                                                                                      |                                                                                                                                                     
tctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTAc                                                                                                                                                        <  1:308935/150‑1 (MQ=255)
 ctGGCGCGTGAACCGCATGTCGATACGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                                                                                      <  2:456353/150‑1 (MQ=255)
 ctGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCTGTACg                                                                                                                                                      >  2:733072/1‑150 (MQ=255)
   ggCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACggg                                                                                                                                                    >  2:1538152/1‑150 (MQ=255)
       cgTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCAACCTGTACGGGTGAg                                                                                                                                                >  1:458560/1‑150 (MQ=255)
        gTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAgg                                                                                                                                               <  1:1538152/150‑1 (MQ=255)
         tGAACCGCATGTCGATATGGGGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGGGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                                                                              <  1:964105/150‑1 (MQ=255)
         tGAACCGCATGTCGATAAGGTGGAGACGGCGAAGGCAGCGGTTTCTATGGCGAAGAAAACCTCGTTCGTGGGCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAggg                                                                                                                                              <  2:1625315/150‑1 (MQ=255)
         tGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAggg                                                                                                                                              >  1:1276043/1‑150 (MQ=255)
          gAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGc                                                                                                                                             <  1:776235/150‑1 (MQ=255)
           aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCg                                                                                                                                            >  1:1001986/1‑150 (MQ=255)
           aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCGTACAGGGTGAGGGca                                                                                                                                            <  1:39901/150‑2 (MQ=255)
           aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGt                                                                                                                                           >  1:1168545/1‑150 (MQ=255)
           aaCCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTATCAATGCCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCg                                                                                                                                            <  1:1526910/150‑1 (MQ=255)
           aaCCGCATGTCGATAAGGTGGAAACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGt                                                                                                                                           >  1:605812/1‑150 (MQ=255)
           aaCCGCATGGCGATAAGGTGGAGAGGCCGATGGCAGCGGTGTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTTGGGGGCGTTGGTGTACAGATTATTGCCCTCACCCGTTCAGGGTGGGGGca                                                                                                                                            <  1:503766/150‑2 (MQ=255)
              cGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGATGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                         <  2:138973/150‑1 (MQ=255)
              cGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCTCTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                         >  2:1090906/1‑150 (MQ=255)
              cGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTag                                                                                                                                         >  2:348345/1‑150 (MQ=255)
                cATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGGTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTagag                                                                                                                                       <  2:1343902/150‑1 (MQ=255)
                   gTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAgagaga                                                                                                                                    <  2:1418072/150‑1 (MQ=255)
                    tCGATATGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGATAACCTCGTGCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAt                                                                                                                                   <  1:665321/150‑1 (MQ=255)
                      gATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTa                                                                                                                                 <  2:1277603/150‑1 (MQ=255)
                         aaGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                              <  2:743532/150‑1 (MQ=255)
                         aaGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATg                                                                                                                              <  2:1615109/150‑1 (MQ=255)
                           ggTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCAACCTGTACGGGTGAGGGCGTATACAGATAAATgcg                                                                                                                            >  2:762268/1‑150 (MQ=255)
                               gagaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCttt                                                                                                                        >  2:1075210/1‑150 (MQ=255)
                                agaCGCCGAAGGCAGCGGTTTCAATGGCGAATAAAACCTCGTTCGTGGTCGTGGTGGTGACTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                       <  1:1371751/150‑1 (MQ=255)
                                agaCGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCtttt                                                                                                                       >  2:347300/1‑150 (MQ=255)
                                    gCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACgg                                                                                                                   >  2:1469664/1‑150 (MQ=255)
                                           gCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGATGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                            >  2:1553727/1‑150 (MQ=255)
                                           gCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGc                                                                                                            >  1:1301841/1‑150 (MQ=255)
                                             aGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCgg                                                                                                          <  2:1476237/150‑1 (MQ=255)
                                              gCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGt                                                                                                         <  2:1594545/150‑1 (MQ=255)
                                               cGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGtt                                                                                                        >  2:1111445/1‑150 (MQ=255)
                                                  tttCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTTAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGttttg                                                                                                     >  2:167129/1‑149 (MQ=255)
                                                    tCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTtctc                                                                                                   <  2:140900/150‑1 (MQ=255)
                                                     cAATGGCGAAGAAAACCTCGTTCGTGGGCGTGGTGGTGAATGCACTACTGACGATTTACGGGGCGGTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTctct                                                                                                  <  1:849616/150‑1 (MQ=255)
                                                       aTGGCGTAGAAAACGTCTTTCGGGGTCGTGGTGGTGATTTCACTACTGACGTATTTCTGGGCGTTGGTGTTAAGATTATTTCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCttt                                                                                                <  2:1279975/150‑1 (MQ=255)
                                                         ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                                                                              >  1:928537/1‑150 (MQ=255)
                                                         ggCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTaa                                                                                              <  1:5785/150‑1 (MQ=255)
                                                           cGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAc                                                                                            >  2:1229766/1‑150 (MQ=255)
                                                            gaagaaAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACa                                                                                           >  2:421810/1‑150 (MQ=255)
                                                                   aCCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                                    <  1:393333/150‑1 (MQ=255)
                                                                   aCCTCGTTAGTGGTCGTTGTTGTGAATGGACTATTGAAGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGTTTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGcc                                                                                    <  1:643457/150‑1 (MQ=255)
                                                                    ccTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCg                                                                                   >  1:1014201/1‑150 (MQ=255)
                                                                      tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                                                                                 <  1:1466725/150‑1 (MQ=255)
                                                                      tcgttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGAt                                                                                 <  2:484656/150‑1 (MQ=255)
                                                                        gttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAc                                                                               <  2:1693784/150‑1 (MQ=255)
                                                                         ttcgtGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATAcc                                                                              <  2:1178548/150‑1 (MQ=255)
                                                                          tcgtggtcgtggTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCt                                                                             >  2:348975/1‑150 (MQ=255)
                                                                                 cgtggtGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGa                                                                      <  2:1244566/150‑1 (MQ=255)
                                                                                       ggtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGcc                                                                >  2:1031865/1‑150 (MQ=255)
                                                                                        gtgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCa                                                               <  2:1458157/150‑1 (MQ=255)
                                                                                         tgATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATNNCCCTCANCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGAGGTTTTCTACTTAAACAGATAGCCAATACATAACACGATCAGACAt                                                              >  1:1565017/1‑150 (MQ=255)
                                                                                         tgATTGCACTACTGACGATTAACGGGGCGTTGGTGTAAAGATTATTGCACTCACACTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCCTTAAACAGATAGCCGATACATAACAAGATAAGACAt                                                              >  1:1571975/1‑150 (MQ=255)
                                                                                          gATTGCACTACTGACGATTTACGGGGCGGGGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATa                                                             <  2:1479199/150‑1 (MQ=255)
                                                                                               cACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCggg                                                        >  1:532543/1‑150 (MQ=255)
                                                                                                   actGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGAATAATGCGCTTTTACGGCTTTGGCGGTATTCTCTTTAAACAGATAGCCGATACCTAACACGAACAGCCATACCGTGATCa                                                    >  1:867121/1‑150 (MQ=255)
                                                                                                    ctGACGATTTACGGGGCGTTGGGGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAg                                                   <  2:1001415/150‑1 (MQ=255)
                                                                                                     tGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAgg                                                  >  1:1220146/1‑150 (MQ=255)
                                                                                                      gACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACACTGTACGGGGGAGGGCCTAGAGAGATTAATGCGCTATAACGGCTTTGGCGGTTTTCACTTTAAACAGATATCCTATACCAAACAAGCTCAGCCATACAGGGATCAGGt                                                 >  2:853999/1‑150 (MQ=255)
                                                                                                        cGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGtat                                               <  2:1301841/150‑1 (MQ=255)
                                                                                                          aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATAc                                             <  2:928537/150‑1 (MQ=255)
                                                                                                          aTTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACACTCTAAGGGTGAGGGCGTAGAGAGATTAATTCTCTTTTACGGCTTTGGCGGATTTATCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGCTGATGtataa                                             >  1:1429058/1‑149 (MQ=255)
                                                                                                           tttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATAcc                                            <  2:430804/150‑1 (MQ=255)
                                                                                                            ttACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATActg                                           <  1:1467042/150‑3 (MQ=255)
                                                                                                            ttACGCGGCGTTGGTGTAAAGATTATTGCGCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCg                                           <  2:925742/150‑1 (MQ=255)
                                                                                                            taaCGGGGCGTTGGTGAAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATCCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCg                                           <  2:218939/148‑1 (MQ=255)
                                                                                                               cGGGGCGTTTGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGaaa                                        <  1:1075416/150‑1 (MQ=255)
                                                                                                                ggggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAAt                                       <  1:787953/150‑1 (MQ=255)
                                                                                                                 gggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                      <  1:1111445/150‑1 (MQ=255)
                                                                                                                 gggCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATc                                      <  1:348975/150‑1 (MQ=255)
                                                                                                                       tGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCAtt                                <  1:527969/150‑1 (MQ=255)
                                                                                                                          tgtAAAGATTATTGCCCTCACCCAGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCt                             >  2:889903/1‑150 (MQ=255)
                                                                                                                           gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGTTCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTg                            >  1:610598/1‑150 (MQ=255)
                                                                                                                           gtAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGGAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATTGCCATTCCTg                            >  1:309420/1‑150 (MQ=255)
                                                                                                                            tAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATAAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTgg                           <  2:888153/150‑1 (MQ=255)
                                                                                                                             taaGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTggg                          >  1:193411/2‑150 (MQ=255)
                                                                                                                                gATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtc                       <  2:675271/150‑1 (MQ=255)
                                                                                                                                    attGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGtcatca                   <  1:140083/150‑1 (MQ=255)
                                                                                                                                     ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTTCGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAg                  <  1:860672/150‑1 (MQ=255)
                                                                                                                                     ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAAGAGATAGCCGATAGCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCGATTCCTGGGGTCATCAg                  <  1:762268/150‑1 (MQ=255)
                                                                                                                                     ttGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAg                  <  1:1621741/150‑1 (MQ=255)
                                                                                                                                       gCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa                <  1:1090906/150‑1 (MQ=255)
                                                                                                                                       ccccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATCCGCTTTTACGGCTTTGGCGGGTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCa                <  1:1531492/149‑1 (MQ=255)
                                                                                                                                        cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt               >  1:190564/1‑150 (MQ=255)
                                                                                                                                        cccTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGAGCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt               >  2:846643/1‑150 (MQ=255)
                                                                                                                                        cccTCACCCGTACAGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCAt               >  1:419318/1‑150 (MQ=255)
                                                                                                                                         ccTCACCCTGGATGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATa              <  2:1334263/150‑1 (MQ=255)
                                                                                                                                          cTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATaa             >  2:1414342/1‑150 (MQ=255)
                                                                                                                                               ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAcc        <  1:1673950/150‑1 (MQ=255)
                                                                                                                                               ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAcc        <  1:1075210/150‑1 (MQ=255)
                                                                                                                                               ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAcc        <  2:702322/150‑1 (MQ=255)
                                                                                                                                               ccTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCAcc        <  1:347300/150‑1 (MQ=255)
                                                                                                                                                 tGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAg      <  1:471963/150‑1 (MQ=255)
                                                                                                                                                 tGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAg      <  1:888880/150‑1 (MQ=255)
                                                                                                                                                  gTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATATTCACCAGt     >  1:1402200/1‑150 (MQ=255)
                                                                                                                                                    aCGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTAc   <  1:582066/150‑1 (MQ=255)
                                                                                                                                                     cGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTAcc  >  1:898553/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCTGGCGCGTGAACCGCATGTCGATAAGGTGGAGACGCCGAAGGCAGCGGTTTCAATGGCGAAGAAAACCTCGTTCGTGGTCGTGGTGGTGATTGCACTACTGACGATTTACGGGGCGTTGGTGTAAAGATTATTGCCCTCACCCTGTACGGGTGAGGGCGTAGAGAGATTAATGCGCTTTTACGGCTTTGGCGGTTTTCTCTTTAAACAGATAGCCGATACCTAACACGATCAGCCATACCGGGATCAGGTATACCGAAATCGCCATTCCTGGGGTCATCAGCATAATCACCAGTACC  >  NC_000913/120009‑120307

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: