Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,710 T→C 5.3% intergenic (+11/‑165) ecnB → / → sugE entericidin B membrane lipoprotein/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,7100TC 5.3% 458.8 / 9.9 150intergenic (+11/‑165)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (68/74);  new base C (3/5);  total (71/79)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 151 total aligned reads displayed.

AGACAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAT  >  NC_000913/4376563‑4376859
                                                                                                                                                   |                                                                                                                                                     
agaCAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTAc                                                                                                                                                     >  2:1132903/1‑150 (MQ=255)
agaCAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTAc                                                                                                                                                     >  2:1324628/1‑150 (MQ=255)
   cAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGAc                                                                                                                                                  <  2:1784712/150‑1 (MQ=255)
      ttGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACTCGATATCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGc                                                                                                                                               >  2:1381252/1‑150 (MQ=255)
        gCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCtg                                                                                                                                             >  2:1602536/1‑150 (MQ=255)
            cGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTc                                                                                                                                         <  1:1123507/150‑1 (MQ=255)
                  tttCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCg                                                                                                                                   <  2:746440/150‑1 (MQ=255)
                  tttCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCg                                                                                                                                   <  2:1644102/150‑1 (MQ=255)
                   ttctgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgt                                                                                                                                  <  1:32751/150‑1 (MQ=255)
                    tctgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgtt                                                                                                                                 <  2:1001544/150‑1 (MQ=255)
                      tgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTTGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgttat                                                                                                                               <  2:1506759/150‑3 (MQ=255)
                      tgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACACTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCgtttg                                                                                                                               >  1:1050532/1‑150 (MQ=255)
                         tctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCgtttgttt                                                                                                                            <  1:1905338/150‑1 (MQ=255)
                          gtgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGGGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgttattgct                                                                                                                           <  2:1593486/149‑7 (MQ=255)
                          ctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgttattgct                                                                                                                           <  1:1820579/150‑7 (MQ=255)
                            ggTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGtttttt                                                                                                                         >  1:1099650/1‑150 (MQ=255)
                             gTGCTTTCAACAGTATTAACTGCCTGCAACACGACGCGTGGGGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTc                                                                                                                        <  2:1699472/150‑1 (MQ=255)
                             gTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTc                                                                                                                        <  2:1131879/150‑1 (MQ=255)
                               gCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTg                                                                                                                      >  1:373625/1‑150 (MQ=255)
                               gCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTg                                                                                                                      >  1:714574/1‑150 (MQ=255)
                                 tttCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCg                                                                                                                    >  1:1369335/1‑150 (MQ=255)
                                    cAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATa                                                                                                                 >  2:140492/1‑150 (MQ=255)
                                    cAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGACTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATa                                                                                                                 <  1:1598088/150‑1 (MQ=255)
                                       cAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTc                                                                                                              <  2:1809195/150‑1 (MQ=255)
                                       cAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTc                                                                                                              <  1:479721/150‑1 (MQ=255)
                                         gTATTAACTGCCTGCAAAACCACGCGTGGCGTTGGTGATGACATTTCTGATGGCGGTAACGCGATATCTGGGGAAGCAACGAACGCGCAGCAATAAGCAATAACCGTACGACAGCTGTGACGTGCCGTTTGTTTTTTCTGCGATAGTcac                                                                                                            >  2:733599/1‑150 (MQ=255)
                                              aaCTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAgg                                                                                                       <  1:208219/150‑1 (MQ=255)
                                               aCTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGt                                                                                                      >  1:771756/1‑150 (MQ=255)
                                                ctgcctgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTa                                                                                                     >  1:1254556/1‑150 (MQ=255)
                                                 tgcctgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTaa                                                                                                    >  1:885882/1‑150 (MQ=255)
                                                 cttcaccAACGCCACGCGTGGTGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTaa                                                                                                    >  1:1426159/7‑150 (MQ=255)
                                                    ctgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAg                                                                                                 <  2:941815/150‑1 (MQ=255)
                                                     tgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGt                                                                                                >  2:268931/1‑150 (MQ=255)
                                                     accAACGCCACGCGTGGTGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGt                                                                                                >  2:599480/3‑150 (MQ=255)
                                                      gcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTATTTCTGCGATAGTCACAAAGGTAATAGtt                                                                                               >  1:661101/1‑150 (MQ=255)
                                                         acacCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGaa                                                                                            <  1:990136/150‑1 (MQ=255)
                                                           accacGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAAt                                                                                          <  2:1825860/150‑1 (MQ=255)
                                                               cgcgTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTccc                                                                                      >  1:39060/1‑150 (MQ=255)
                                                               cgcgGGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGGTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTccc                                                                                      <  2:251635/150‑1 (MQ=255)
                                                                  ggcgCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATGTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTg                                                                                   <  1:1158732/147‑1 (MQ=255)
                                                                      cGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCAc                                                                               >  2:1671344/1‑150 (MQ=255)
                                                                      cGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTACCCTGCCAc                                                                               >  2:1416846/1‑150 (MQ=255)
                                                                          ggTGAAGACATTTCTTATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                           >  2:1356074/1‑150 (MQ=255)
                                                                          ggTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                           >  1:374009/1‑150 (MQ=255)
                                                                           gTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGc                                                                          <  2:809003/150‑1 (MQ=255)
                                                                            tGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCa                                                                         >  2:1578995/1‑150 (MQ=255)
                                                                              aaGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCaaa                                                                       >  1:614005/1‑150 (MQ=255)
                                                                               aGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCaaaa                                                                      >  2:459418/1‑150 (MQ=255)
                                                                                gACATTTCTGATGGCGGTACCGCGATTTCTGGCTCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGGCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAat                                                                     <  2:1153911/150‑1 (MQ=255)
                                                                                     ttCTGATGGCTGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCg                                                                <  2:1473807/150‑1 (MQ=255)
                                                                                     ttCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCg                                                                <  1:1180418/150‑1 (MQ=255)
                                                                                       cTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGtt                                                              >  2:561937/1‑150 (MQ=255)
                                                                                         gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTc                                                                                                                        <  1:1792429/90‑1 (MQ=255)
                                                                                         gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCa                                                            <  1:1324628/150‑1 (MQ=255)
                                                                                          aTGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCaa                                                           >  2:73466/1‑150 (MQ=255)
                                                                                           tGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAAc                                                          <  1:1003581/150‑1 (MQ=255)
                                                                                            ggCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAAcc                                                         >  2:1349254/1‑150 (MQ=255)
                                                                                             gCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCa                                                        <  2:761820/150‑1 (MQ=255)
                                                                                                    cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCCCAAAGGTAATAGTTGAAATTCCCCTGCCCCCTGGCAAAATATCCGTTCCACCATCAGCtt                                                 >  1:926270/1‑150 (MQ=255)
                                                                                                    cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                 >  1:1535874/1‑150 (MQ=255)
                                                                                                    cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                 >  2:611829/1‑150 (MQ=255)
                                                                                                    cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                 >  2:308197/1‑150 (MQ=255)
                                                                                                    cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                 >  1:1460605/1‑150 (MQ=255)
                                                                                                      cgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTg                                               <  1:218588/150‑1 (MQ=255)
                                                                                                        agttCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTGCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                             <  1:1689866/148‑1 (MQ=255)
                                                                                                        aTTTCTTGCGCAGCAACGAAAGCGCAGCAATAAGAAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTTCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                             <  1:450052/150‑1 (MQ=255)
                                                                                                        aTTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                             <  1:1870668/150‑1 (MQ=255)
                                                                                                        aTTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                             <  1:1311972/150‑1 (MQ=255)
                                                                                                           tCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                          <  2:384385/150‑1 (MQ=255)
                                                                                                            cTGGCGCAGCAACGAAAGCGCAGCAATAATCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgac                                         <  1:736067/150‑1 (MQ=255)
                                                                                                             tGGCGCATCAACGAAAGCGCAGCAATAAGCAATAACGGTACGTCAGCTGTGTTGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacg                                        <  1:1512547/150‑1 (MQ=255)
                                                                                                             tGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacg                                        >  2:1239888/1‑150 (MQ=255)
                                                                                                              ggCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                       <  2:411564/150‑1 (MQ=255)
                                                                                                              ggCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTGTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACATGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                       >  2:384808/1‑150 (MQ=255)
                                                                                                               gcgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                      <  2:227047/150‑1 (MQ=255)
                                                                                                                cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCGCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTcc                                     <  1:766525/150‑1 (MQ=255)
                                                                                                                cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGAcc                                     <  2:487592/150‑1 (MQ=255)
                                                                                                                  cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTg                                   <  2:689712/150‑1 (MQ=255)
                                                                                                                  cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTg                                   <  2:233542/150‑1 (MQ=255)
                                                                                                                    gcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCACCCATCAGCTTTGCAGGACGACCTGCa                                 >  1:1816335/1‑150 (MQ=255)
                                                                                                                          aaaGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                           <  1:1123276/150‑1 (MQ=255)
                                                                                                                          aaaGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                           <  2:679671/150‑1 (MQ=255)
                                                                                                                              cgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCtt                       >  2:399801/1‑150 (MQ=255)
                                                                                                                               gcagcaATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGATTTTTCTGCGATAGAAACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAAAATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCCAACGccacgtt                      >  2:131328/1‑145 (MQ=255)
                                                                                                                                  gcaataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                   >  2:1256687/1‑150 (MQ=255)
                                                                                                                                   caataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAc                  <  2:374009/150‑1 (MQ=255)
                                                                                                                                      taagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg               >  1:1533771/1‑150 (MQ=255)
                                                                                                                                      taagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg               >  1:1526207/1‑150 (MQ=255)
                                                                                                                                      aaagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg               <  1:12219/149‑1 (MQ=255)
                                                                                                                                       aagcaataaCGGTACGACAGCTGAGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCggg              <  1:266197/150‑1 (MQ=255)
                                                                                                                                        aaaacaaaCGGCACGACACAGCTGTCGTACCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg             >  1:1245251/7‑150 (MQ=255)
                                                                                                                                         gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa            <  1:134387/150‑1 (MQ=255)
                                                                                                                                         gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa            <  2:771756/150‑1 (MQ=255)
                                                                                                                                         aaacaaaCGGCACGACACAGCTGTCGTACCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGACCTGCAAACTCCTATTTTCACCGGGGa            >  2:1302703/6‑150 (MQ=255)
                                                                                                                                         aaacaaaCGGCACGACACAGCTGTCGTACCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa            >  1:1220916/6‑150 (MQ=255)
                                                                                                                                          caataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGAc           >  2:531890/1‑150 (MQ=255)
                                                                                                                                             taaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGc        <  1:1596856/150‑1 (MQ=255)
                                                                                                                                               aCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGccc      <  2:1067392/150‑1 (MQ=255)
                                                                                                                                               aCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGccc      >  1:1022443/1‑150 (MQ=255)
                                                                                                                                                cGGTACGTCAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCAAAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGcccc     <  1:736980/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGACAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAT  >  NC_000913/4376563‑4376859

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