Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,822 A→T 15.7% intergenic (+123/‑53) 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,8220AT15.7% 272.3 / 42.9 115intergenic (+123/‑53)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (45/51);  new base T (9/9);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

Only 100 of 122 total aligned reads displayed.

CGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATG  >  NC_000913/4376724‑4376921
                                                                                                  |                                                                                                    
cGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                     >  1:2592007/1‑100 (MQ=255)
  tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                   >  2:2336591/1‑100 (MQ=255)
    ccGTTTGTTTTTTCTGCGATAGTCACAAAGGCAATAGGTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGc                                                                                                 <  1:1747281/100‑1 (MQ=255)
     cGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                                                                                <  2:83819/100‑1 (MQ=255)
     cGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                                                                                >  2:644670/1‑100 (MQ=255)
     cGTTTGTTTTTTCTGCGATAGTCACAAAGCTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                                                                                <  1:1562061/100‑1 (MQ=255)
        ttgtttTTTCTGCGAAAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaag                                                                                             <  2:2549874/100‑2 (MQ=255)
         tgtttTTTCTGCGATAGTCACAAAGTTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagc                                                                                            <  2:2376721/100‑3 (MQ=255)
         tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagc                                                                                            <  2:194425/100‑3 (MQ=255)
         tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                                                                            >  2:758499/1‑100 (MQ=255)
         tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCGGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAAc                                                                                             <  1:1797053/99‑1 (MQ=255)
          gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagct                                                                                           <  2:2507197/100‑4 (MQ=255)
          gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagct                                                                                           <  1:543466/100‑4 (MQ=255)
          gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGc                                                                                           >  2:2076531/1‑100 (MQ=255)
           ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                          >  1:1099465/1‑100 (MQ=255)
           ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                          >  1:2059908/1‑100 (MQ=255)
           ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                          <  1:176919/100‑1 (MQ=255)
           ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATACGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagt                                                                                            <  1:1783780/98‑3 (MQ=255)
            tttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCct                                                                                         >  2:2581291/1‑100 (MQ=255)
               ttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgatgg                                                                                      <  2:1921716/100‑9 (MQ=38)
               ttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCtt                                                                                      >  1:2076994/1‑100 (MQ=255)
                tCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCttt                                                                                     >  2:141047/1‑100 (MQ=255)
                tCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCttt                                                                                     <  2:1811798/100‑1 (MQ=255)
                 cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCtttt                                                                                    >  2:87804/1‑100 (MQ=255)
                 cTGCGATAGTCACAAAGGTAATAGATGAAATTCCCCCGCCCCCTGGCAAAATACCCCTGCACCCATCAGCCTTGCAGGCCGACCTGCAAACGCCTCtttt                                                                                    <  1:215211/100‑1 (MQ=255)
                   gCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                                                                                  <  1:1342692/100‑1 (MQ=255)
                    cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAc                                                                                 >  2:385903/1‑100 (MQ=255)
                    cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTAGCAGGACGACCTGCAAACGCCTCTTTTCAc                                                                                 >  2:909493/1‑100 (MQ=255)
                      aTAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑Cg                                                                              >  1:2095022/1‑100 (MQ=255)
                       tAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑Cgg                                                                             >  1:1973673/1‑100 (MQ=255)
                         gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑Cgggg                                                                           >  2:1320536/1‑100 (MQ=255)
                         gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑Cgggg                                                                           >  1:1973705/1‑100 (MQ=255)
                           cacaAATGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGAc                                                                         >  2:823845/1‑100 (MQ=255)
                           cacaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGAc                                                                         >  2:2534192/1‑100 (MQ=255)
                           cacaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACATGCAAACGCCTCTTTTCAC‑CGGGGAc                                                                         >  1:246551/1‑100 (MQ=255)
                             caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACgg                                                                       >  2:1962605/1‑100 (MQ=255)
                             caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACg                                                                        <  2:246553/99‑1 (MQ=255)
                              aaaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACg                                                                        <  2:2344759/98‑1 (MQ=255)
                              aaaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGc                                                                      <  2:1653840/100‑1 (MQ=255)
                               aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGcc                                                                     <  2:614130/100‑1 (MQ=255)
                               aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGcc                                                                     <  2:851435/100‑1 (MQ=255)
                                 gtgaaTAGTTGAAGTTCCCCTGAAAAATGACAAAATATACGTTCAACCATAAGCTGTGCAGGACGACCTGCAAACGCCTC‑TTTCACACGGGAACGGcc                                                                     >  1:98178/4‑98 (MQ=255)
                                 ggTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGcccc                                                                   >  2:1112366/1‑100 (MQ=255)
                                 ggTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGcccc                                                                   >  2:1112067/1‑100 (MQ=255)
                                 ggTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGcccc                                                                   <  2:154847/100‑1 (MQ=255)
                                       aGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATtct                                                             >  1:2334669/1‑100 (MQ=255)
                                         ttGAAATTCCTCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTcc                                                           <  1:1703920/100‑1 (MQ=255)
                                         ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTcc                                                           <  1:1645776/100‑1 (MQ=255)
                                         ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTcc                                                           <  2:410963/100‑1 (MQ=255)
                                         ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTcc                                                           <  2:1405521/100‑1 (MQ=255)
                                          tGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGCCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCg                                                          >  2:1298889/1‑100 (MQ=255)
                                          tGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCg                                                          <  1:1105953/100‑1 (MQ=255)
                                           gAAATTCCCCTGCCTCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCgg                                                         <  2:1542671/100‑1 (MQ=255)
                                           gAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAA‑AGCGGCCGGACAAAATTCTCCg                                                          >  2:476742/1‑99 (MQ=255)
                                            aaaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGa                                                        <  1:940138/100‑1 (MQ=255)
                                             aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAg                                                       >  2:1406477/1‑100 (MQ=255)
                                             aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAg                                                       <  2:1105078/100‑1 (MQ=255)
                                                tCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCt                                                    <  1:660791/100‑1 (MQ=255)
                                                  cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAT‑CGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                  >  1:1328412/1‑100 (MQ=255)
                                                  cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                  >  1:1760086/1‑100 (MQ=255)
                                                     cgatctcTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGatat                                               <  2:939930/94‑1 (MQ=255)
                                                      gCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATg                                              >  1:465496/1‑100 (MQ=255)
                                                        cACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGACTGATATGTc                                            <  1:148404/100‑1 (MQ=255)
                                                           cTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGatatagatcg                                         >  1:939927/1‑94 (MQ=255)
                                                             ggCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGa                                       <  2:1037681/100‑1 (MQ=255)
                                                              gCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAt                                      >  1:1192960/1‑100 (MQ=255)
                                                               cAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAtt                                     >  1:1325328/1‑100 (MQ=255)
                                                                aaaaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTa                                    <  1:2403084/100‑1 (MQ=255)
                                                                   atatCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCTCCGATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                 <  1:572921/100‑1 (MQ=255)
                                                                   atatCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                 <  2:1812002/100‑1 (MQ=255)
                                                                    tatCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACTGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCtt                                <  2:962334/100‑1 (MQ=255)
                                                                       ccGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCTTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAg                              <  1:2512108/99‑1 (MQ=255)
                                                                       ccGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGt                             >  2:1475560/1‑100 (MQ=255)
                                                                          ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTAt                          <  2:1952549/100‑1 (MQ=255)
                                                                           tCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTAtt                         <  2:170265/100‑1 (MQ=255)
                                                                            cAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTg                        >  2:2404182/1‑100 (MQ=255)
                                                                              aCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                      >  2:411158/1‑100 (MQ=255)
                                                                               ccATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTg                     >  1:156667/1‑100 (MQ=255)
                                                                                cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                    <  2:1965174/100‑1 (MQ=255)
                                                                                cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                    <  1:1123303/100‑1 (MQ=255)
                                                                                cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                    >  2:1632827/1‑100 (MQ=255)
                                                                                   ggcgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTct                 >  2:1995391/5‑100 (MQ=255)
                                                                                    gcgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctg                >  2:86915/4‑100 (MQ=255)
                                                                                     gCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               <  2:1069691/100‑1 (MQ=255)
                                                                                     gCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               <  1:411156/100‑1 (MQ=255)
                                                                                     cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  1:2135735/3‑100 (MQ=255)
                                                                                     cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  1:702434/3‑100 (MQ=255)
                                                                                     cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  2:2352307/3‑100 (MQ=255)
                                                                                     cgtttGCAGGTCATCCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  2:1706007/3‑100 (MQ=255)
                                                                                      cTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgct              <  1:2348977/100‑1 (MQ=255)
                                                                                       tttGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgctg             >  2:1555685/1‑100 (MQ=255)
                                                                                       tttGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATAGGTCCTGGATTATCTTAGTTATTGCTGGTctgctg             >  1:1362664/1‑100 (MQ=255)
                                                                                        ttGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTgg            <  1:1214787/100‑1 (MQ=255)
                                                                                            aGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGt        <  1:1897181/100‑1 (MQ=255)
                                                                                            aGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGt        <  1:444133/100‑1 (MQ=255)
                                                                                             ggACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       <  2:1731456/100‑1 (MQ=255)
                                                                                              gacgacCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       <  1:1592777/99‑1 (MQ=255)
                                                                                                cgacCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTa    <  2:423241/100‑1 (MQ=255)
                                                                                                 gacCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTAt   <  2:1172784/100‑1 (MQ=255)
                                                                                                  acCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATg  <  2:603407/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
CGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAC‑CGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATG  >  NC_000913/4376724‑4376921

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑