Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,720 G→C 5.1% intergenic (+21/‑155) 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,7200GC 5.1% 523.4 / 12.6 157intergenic (+21/‑155)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (69/80);  new base C (3/5);  total (72/85)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 162 total aligned reads displayed.

GCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAG  >  NC_000913/4376571‑4376868
                                                                                                                                                     |                                                                                                                                                    
     gATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCg                                                                                                                                                 <  1:225540/150‑1 (MQ=255)
        cTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGc                                                                                                                                              >  1:977695/1‑150 (MQ=255)
              tgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTTGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGCACGACACAGCTGTCGTACCgttat                                                                                                                                        <  2:1506759/150‑3 (MQ=255)
              tgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGAAAGCTGTGTCGTGCCgtttg                                                                                                                                        >  1:1201349/1‑150 (MQ=255)
              tgttctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACACTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCgtttg                                                                                                                                        >  1:1050532/1‑150 (MQ=255)
                 tctgGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATTGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCgtttgttt                                                                                                                                     <  2:246199/150‑1 (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)
                     gTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTc                                                                                                                                 <  2:185382/150‑1 (MQ=255)
                       gCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTg                                                                                                                               >  1:714574/1‑150 (MQ=255)
                       gCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTg                                                                                                                               <  2:374590/150‑1 (MQ=255)
                        cTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGc                                                                                                                              <  2:441677/150‑1 (MQ=255)
                         tttCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCg                                                                                                                             >  1:1369335/1‑150 (MQ=255)
                            cAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATa                                                                                                                          >  2:140492/1‑150 (MQ=255)
                            cAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGACTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATa                                                                                                                          <  1:1598088/150‑1 (MQ=255)
                             aaCAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAg                                                                                                                         <  1:452186/150‑1 (MQ=255)
                               cAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTc                                                                                                                       <  2:1809195/150‑1 (MQ=255)
                                 gTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTcac                                                                                                                     >  1:1528986/1‑150 (MQ=255)
                                 gTATTAACTGCCTGCAAAACCACGCGTGGCGTTGGTGATGACATTTCTGATGGCGGTAACGCGATATCTGGGGAAGCAACGAACGCGCAGCAATAAGCAATAACCGTACGACAGCTGTGACGTGCCGTTTGTTTTTTCTGCGATAGTcac                                                                                                                     >  2:733599/1‑150 (MQ=255)
                                      aaCTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAgg                                                                                                                <  1:208219/150‑1 (MQ=255)
                                        ctgcctgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTa                                                                                                              >  1:1254556/1‑150 (MQ=255)
                                         cttcaccAACGCCACGCGTGGTGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTaa                                                                                                             >  1:1426159/7‑150 (MQ=255)
                                            ctgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAg                                                                                                          <  2:941815/150‑1 (MQ=255)
                                             tgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGt                                                                                                         >  2:522934/1‑150 (MQ=255)
                                             tgcAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGt                                                                                                         <  1:1241538/150‑1 (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)
                                                       cgcgGGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGGTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTccc                                                                                               <  2:251635/150‑1 (MQ=255)
                                                          ggcgCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATGTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTg                                                                                            <  1:1158732/147‑1 (MQ=255)
                                                              cGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCAc                                                                                        >  2:1671344/1‑150 (MQ=255)
                                                                  ggTGAAGACATTTCTTATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                                    >  2:1356074/1‑150 (MQ=255)
                                                                  ggTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                                    >  1:374009/1‑150 (MQ=255)
                                                                  ggTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                                    >  1:765668/1‑150 (MQ=255)
                                                                  ggTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGACGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTgg                                                                                    <  2:1405056/150‑1 (MQ=255)
                                                                   gTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGc                                                                                   <  2:809003/150‑1 (MQ=255)
                                                                      aaGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCaaa                                                                                >  1:614005/1‑150 (MQ=255)
                                                                        gACATTTCTGATGGCGGTACCGCGATTTCTGGCTCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGGCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAat                                                                              <  2:1153911/150‑1 (MQ=255)
                                                                             ttCTGATGGCTGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCg                                                                         <  2:1473807/150‑1 (MQ=255)
                                                                               cTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGtt                                                                       >  2:561937/1‑150 (MQ=255)
                                                                                tGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTc                                                                      <  1:841192/150‑1 (MQ=255)
                                                                                 gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTc                                                                                                                                 <  1:1792429/90‑1 (MQ=255)
                                                                                 gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCa                                                                     <  1:1324628/150‑1 (MQ=255)
                                                                                 gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCa                                                                     <  2:1589389/150‑1 (MQ=255)
                                                                                 gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACATGGCAAAATATCCGTACa                                                                     >  2:736980/1‑150 (MQ=255)
                                                                                 gATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCAACCTGGCAAATTATCCGTTCa                                                                     >  1:335913/1‑150 (MQ=255)
                                                                                   tGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAAc                                                                   <  1:1003581/150‑1 (MQ=255)
                                                                                     gCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCa                                                                 <  2:761820/150‑1 (MQ=255)
                                                                                            cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                          <  2:1675094/150‑1 (MQ=255)
                                                                                            cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                          >  2:308197/1‑150 (MQ=255)
                                                                                            cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                          >  1:1535874/1‑150 (MQ=255)
                                                                                            cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCtt                                                          >  2:611829/1‑150 (MQ=255)
                                                                                            cgcgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGAtt                                                          >  1:1127817/1‑150 (MQ=255)
                                                                                              cgATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTg                                                        >  2:719838/1‑150 (MQ=255)
                                                                                                agttCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTGCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                                      <  1:1689866/148‑1 (MQ=255)
                                                                                                aTTTCTTGCGCAGCAACGAAAGCGCAGCAATAAGAAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTTCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                                      <  1:450052/150‑1 (MQ=255)
                                                                                                aTTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCa                                                      <  1:1870668/150‑1 (MQ=255)
                                                                                                   tCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                                   <  2:384385/150‑1 (MQ=255)
                                                                                                    cTGGCGCAGCAACGAAAGCGCAGCAATAATCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgac                                                  <  1:736067/150‑1 (MQ=255)
                                                                                                     tGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacg                                                 >  2:1239888/1‑150 (MQ=255)
                                                                                                      ggCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                <  2:411564/150‑1 (MQ=255)
                                                                                                       gcgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                               <  2:227047/150‑1 (MQ=255)
                                                                                                        cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCGCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTcc                                              <  1:766525/150‑1 (MQ=255)
                                                                                                        cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGAcc                                              <  2:487592/150‑1 (MQ=255)
                                                                                                            gcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCACCCATCAGCTTTGCAGGACGACCTGCa                                          >  1:1816335/1‑150 (MQ=255)
                                                                                                            gcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                          <  1:1480064/150‑1 (MQ=255)
                                                                                                            gcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                          <  1:379219/150‑1 (MQ=255)
                                                                                                                agAAATCGCAGCAATTAGCAATAACTGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                      <  2:1633583/149‑1 (MQ=255)
                                                                                                                  aaaGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                    <  1:1123276/150‑1 (MQ=255)
                                                                                                                  aaaGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                    <  2:679671/150‑1 (MQ=255)
                                                                                                                     gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctct                                 <  2:1341342/150‑1 (MQ=255)
                                                                                                                     gcgcAGCAATAAGCAATAAAGGTACGACAGCTGTGTCGGGCCGTTTGTTTTTTCTGCGATAGTGACAATGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctct                                 <  2:1153407/150‑1 (MQ=255)
                                                                                                                       gcagcaATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGATTTTTCTGCGATAGAAACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAAAATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCCAACGccacgtt                               >  2:131328/1‑145 (MQ=255)
                                                                                                                          gcaataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                            >  2:1256687/1‑150 (MQ=255)
                                                                                                                              taagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg                        >  1:1526207/1‑150 (MQ=255)
                                                                                                                               aagcaataaCGGTACGACAGCTGAGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCggg                       <  1:266197/150‑1 (MQ=255)
                                                                                                                                agcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                      >  1:1264994/1‑150 (MQ=255)
                                                                                                                                agcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                      >  2:39759/1‑150 (MQ=255)
                                                                                                                                 gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa                     <  2:771756/150‑1 (MQ=255)
                                                                                                                                 gcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa                     <  1:134387/150‑1 (MQ=255)
                                                                                                                                 aaacaaaCGGCACGACACAGCTGTCGTACCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGACCTGCAAACTCCTATTTTCACCGGGGa                     >  2:1302703/6‑150 (MQ=255)
                                                                                                                                 aaacaaaCGGCACGACACAGCTGTCGTACCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGa                     >  1:1220916/6‑150 (MQ=255)
                                                                                                                                     taaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGc                 <  1:1596856/150‑1 (MQ=255)
                                                                                                                                       aCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGccc               >  1:1022443/1‑150 (MQ=255)
                                                                                                                                        cGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGcccc              <  2:1128626/150‑1 (MQ=255)
                                                                                                                                            aCGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAtt          <  2:271568/150‑1 (MQ=255)
                                                                                                                                            aCGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAtt          <  2:1528986/150‑1 (MQ=255)
                                                                                                                                             cGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc         <  2:104163/150‑1 (MQ=255)
                                                                                                                                             cGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGAAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc         <  1:1654577/150‑1 (MQ=255)
                                                                                                                                             cGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGACACCTGGCAAAATATCCGTTCAGCCATCCGCTTTGCAGGTCGACCTGGAAACGCCTCTTTTCACCGGGGACGGCCGCAATtc         >  2:284070/1‑150 (MQ=255)
                                                                                                                                              gACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtct        >  2:1809812/1‑150 (MQ=255)
                                                                                                                                               aCAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc       >  2:1112330/1‑150 (MQ=255)
                                                                                                                                                cAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTcc      >  1:874304/1‑150 (MQ=255)
                                                                                                                                                  tcTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCgg    <  1:141671/149‑1 (MQ=255)
                                                                                                                                                  gCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCgg    <  2:1738365/150‑1 (MQ=255)
                                                                                                                                                    tgtgTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAg  >  2:744930/1‑150 (MQ=255)
                                                                                                                                                    tgtgGCGTGCCGTTTGTTTTTTCTGCGATAGTCACAACGGAAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAg  <  1:891155/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAG  >  NC_000913/4376571‑4376868

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: