Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,819 A→T 7.3% intergenic (+120/‑56) 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,8190AT 7.3% 367.3 / 11.6 106intergenic (+120/‑56)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (49/49);  new base T (4/4);  total (53/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCC  >  NC_000913/4376675‑4376966
                                                                                                                                                |                                                                                                                                                   
cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  <  1:834427/148‑1 (MQ=255)
cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  <  1:30171/148‑1 (MQ=255)
cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  >  2:512341/1‑148 (MQ=255)
    gcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTg                                                                                                                                              <  2:440882/148‑1 (MQ=255)
     caACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGc                                                                                                                                             >  2:570233/1‑148 (MQ=255)
       aCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaaa                                                                                                                                          <  4:144037/149‑1 (MQ=255)
       aCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaa                                                                                                                                           >  2:637788/1‑148 (MQ=255)
        cGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaaa                                                                                                                                          >  2:35832/1‑148 (MQ=255)
         gAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaag                                                                                                                                         <  1:876706/148‑2 (MQ=255)
              cgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctct                                                                                                                                   <  4:409528/149‑1 (MQ=255)
               gcagcaATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctct                                                                                                                                   >  2:642113/1‑148 (MQ=255)
                      taagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAcc                                                                                                                            <  1:457503/148‑1 (MQ=255)
                                  gTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGcccc                                                                                                                >  2:492246/1‑148 (MQ=255)
                                   tACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCaa                                                                                                              >  4:277714/1‑149 (MQ=255)
                                      gACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc                                                                                                           >  4:741374/1‑149 (MQ=255)
                                        cAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc                                                                                                         <  3:82805/149‑1 (MQ=255)
                                          gCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCg                                                                                                       <  3:126615/149‑1 (MQ=255)
                                            tgtgTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCgg                                                                                                      >  1:315559/1‑148 (MQ=255)
                                            tgtgTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGa                                                                                                     >  4:139774/1‑149 (MQ=255)
                                                 cGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCt                                                                                                 >  2:559378/1‑148 (MQ=255)
                                                   tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  <  3:484337/97‑1 (MQ=255)
                                                   tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                               >  2:807998/1‑148 (MQ=255)
                                                   tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                               >  1:665923/1‑148 (MQ=255)
                                                    gCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              <  1:345995/148‑1 (MQ=255)
                                                    gCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              >  2:345995/1‑148 (MQ=255)
                                                     ccGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGCTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGata                                                                                             <  1:17922/148‑1 (MQ=255)
                                                        tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         >  4:227099/1‑149 (MQ=255)
                                                        tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         >  3:362295/1‑149 (MQ=255)
                                                        tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         <  4:226668/149‑1 (MQ=255)
                                                        tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         >  3:486396/1‑149 (MQ=255)
                                                         ttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTcc                                                                                        >  4:682221/1‑149 (MQ=255)
                                                           gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTg                                                                                      <  3:139774/149‑1 (MQ=255)
                                                            ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTgg                                                                                     >  3:512375/1‑149 (MQ=255)
                                                              ttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGa                                                                                    <  1:603896/148‑1 (MQ=255)
                                                                  cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGTCAAAATATACGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                                                               <  3:580958/149‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCtt                                                                             <  1:35832/148‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCtt                                                                             <  1:492246/148‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  4:486396/149‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  3:376650/149‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  3:73890/149‑1 (MQ=255)
                                                                     cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  4:389261/149‑1 (MQ=255)
                                                                      gATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            >  1:406083/1‑148 (MQ=255)
                                                                          gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAGCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTAt                                                                       >  3:780320/1‑149 (MQ=255)
                                                                           tCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              <  4:631506/125‑1 (MQ=255)
                                                                           tCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              >  3:631506/1‑125 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc                                                                                                         <  4:660293/112‑1 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc                                                                                                         >  3:660293/1‑112 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  4:52876/1‑149 (MQ=255)
                                                                             acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  4:483182/1‑149 (MQ=255)
                                                                              caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                                                                   >  4:734921/1‑149 (MQ=255)
                                                                                aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTg                                                                  <  1:570233/148‑1 (MQ=255)
                                                                                  ggTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTc                                                               <  3:112825/149‑1 (MQ=255)
                                                                                    tAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTct                                                              >  2:512966/1‑148 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCtt                                                                             <  2:19135/127‑1 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCtt                                                                             >  1:19135/1‑127 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                       >  3:899679/1‑149 (MQ=255)
                                                                                          ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                       <  3:682221/149‑1 (MQ=255)
                                                                                           cgAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                       <  1:807998/147‑1 (MQ=255)
                                                                                             aaaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg                                                    <  3:235904/149‑1 (MQ=255)
                                                                                               aTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCTAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGt                                                  >  4:322637/1‑149 (MQ=255)
                                                                                                 tCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTa                                                 <  1:809135/148‑1 (MQ=255)
                                                                                                   cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       >  3:505727/1‑108 (MQ=255)
                                                                                                   cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       <  4:505727/108‑1 (MQ=255)
                                                                                                   cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATgg                                              >  3:525784/1‑149 (MQ=255)
                                                                                                     cTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTACCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATggg                                             >  1:12646/1‑148 (MQ=255)
                                                                                                         cACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGtt                                        <  3:52876/149‑1 (MQ=255)
                                                                                                          aCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGtt                                        >  1:239323/1‑148 (MQ=255)
                                                                                                          aCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGtt                                        <  1:637788/148‑1 (MQ=255)
                                                                                                          aCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGtt                                        <  2:221296/148‑1 (MQ=255)
                                                                                                          aCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTg                                       >  4:741506/1‑149 (MQ=255)
                                                                                                              ggCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGcc                                    >  1:783876/1‑148 (MQ=255)
                                                                                                                cAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTGACCGAGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTAGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTg                                  >  1:268819/1‑148 (MQ=255)
                                                                                                                   aaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                                                               <  1:69743/100‑1 (MQ=255)
                                                                                                                   aaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAat                              <  4:525784/149‑1 (MQ=255)
                                                                                                                      atCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAAtata                           <  3:741506/149‑1 (MQ=255)
                                                                                                                         cGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATAcc                         <  1:512966/148‑1 (MQ=255)
                                                                                                                          gTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTTGAAGTGGTATg                                               <  3:582959/125‑1 (MQ=255)
                                                                                                                          gTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATg                                               >  4:582959/1‑125 (MQ=255)
                                                                                                                          gTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCa                       <  3:277714/149‑1 (MQ=255)
                                                                                                                           ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATGCCCAc                      >  4:582138/1‑149 (MQ=255)
                                                                                                                           ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCa                       <  2:239323/148‑1 (MQ=255)
                                                                                                                              aaCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCGCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGcc                                    <  4:651705/132‑1 (MQ=255)
                                                                                                                              aaCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGcc                                    >  3:651705/1‑132 (MQ=255)
                                                                                                                               aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                                                                   <  2:404158/100‑1 (MQ=255)
                                                                                                                               aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGc                   >  1:676094/1‑148 (MQ=255)
                                                                                                                                 cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCttt                >  3:177144/1‑149 (MQ=255)
                                                                                                                                  aTCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTa               >  4:493782/1‑149 (MQ=255)
                                                                                                                                      gCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCg           >  4:660012/1‑149 (MQ=255)
                                                                                                                                       gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt          <  4:355550/148‑1 (MQ=255)
                                                                                                                                       gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt          >  3:355550/2‑149 (MQ=255)
                                                                                                                                        tttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATgg                                              <  3:143420/112‑1 (MQ=255)
                                                                                                                                        tttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATgg                                              >  4:143420/1‑112 (MQ=255)
                                                                                                                                        tttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGtt         <  4:512375/149‑1 (MQ=255)
                                                                                                                                        tttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGtt         >  3:721025/1‑149 (MQ=255)
                                                                                                                                         ttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGtt         >  2:285187/1‑148 (MQ=255)
                                                                                                                                         ttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGttt        <  3:660012/149‑1 (MQ=255)
                                                                                                                                         ttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGtt         >  2:507376/1‑148 (MQ=255)
                                                                                                                                            cAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGa      <  1:559378/148‑1 (MQ=255)
                                                                                                                                              ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACg    >  2:617127/1‑148 (MQ=255)
                                                                                                                                               gacgacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGcc  >  4:235758/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCC  >  NC_000913/4376675‑4376966

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