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

Only 100 of 191 total aligned reads displayed.

CTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCC  >  NC_000913/4376671‑4376966
                                                                                                                                                    |                                                                                                                                                   
cTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                                                                                                                                     <  1:1507719/149‑1 (MQ=255)
cTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                                                                                                                                     <  4:559239/149‑1 (MQ=255)
 tGGCGCAGAAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgac                                                                                                                                                    <  1:1003759/149‑1 (MQ=255)
   gcgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  >  1:1046295/1‑149 (MQ=255)
   gcgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  <  3:1454718/149‑1 (MQ=255)
   gcgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacga                                                                                                                                                  <  3:792394/149‑1 (MQ=255)
    cgcAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                                                                 <  2:8552/149‑1 (MQ=255)
      cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                                                               >  4:632919/1‑149 (MQ=255)
      cagcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                                                               <  2:813441/149‑1 (MQ=255)
       agcaACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTg                                                                                                                                              >  1:1031151/1‑149 (MQ=255)
          aaCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaa                                                                                                                                           >  1:575204/1‑149 (MQ=255)
           aCGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCaaa                                                                                                                                          >  1:510101/1‑149 (MQ=255)
             gAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                                                                                                                        >  4:1329719/1‑149 (MQ=255)
                 gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctc                                                                                                                                    <  2:822066/149‑1 (MQ=255)
                 gcgcAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCctc                                                                                                                                    <  4:371468/149‑1 (MQ=255)
                   gcggcaataagcaataaCGGTACGACAGCTGTGTCGTGCCGGTTGTTTTTTCGGCGATAGTCATAAAGGTAATAGTTGAAATTCCCCTTCCACCTGGCAAAATATCCGTTCAACCATCAGCTTGGCAGGACGACCTGCAAACGCCtttt                                                                                                                                  <  3:1043894/146‑1 (MQ=255)
                   gcagcaATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCtt                                                                                                                                  <  1:320583/149‑1 (MQ=255)
                        aataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAc                                                                                                                             <  1:1354793/149‑1 (MQ=255)
                         ataagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACACCTCTTTTCAcc                                                                                                                            >  3:1317370/1‑149 (MQ=255)
                           aagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg                                                                                                                          <  1:472535/149‑1 (MQ=255)
                           aagcaataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg                                                                                                                          >  4:911132/1‑149 (MQ=255)
                               aataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGAc                                                                                                                      >  1:243526/1‑149 (MQ=255)
                                ataaCGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACg                                                                                                                     >  4:118238/1‑149 (MQ=255)
                                          gACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTACACCGGGGACGGCCCCAATtc                                                                                                           >  4:1479655/1‑149 (MQ=255)
                                           aCAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAtt                                                                                                            >  3:602589/1‑147 (MQ=255)
                                           aCAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAtt                                                                                                            <  4:602589/147‑1 (MQ=255)
                                               cTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACg                                                                                                                     >  2:104542/1‑134 (MQ=255)
                                               cTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCgg                                                                                                      >  2:249175/1‑149 (MQ=255)
                                                tgtgTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGa                                                                                                     >  2:496454/1‑149 (MQ=255)
                                                       tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAATATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              >  1:1014630/1‑149 (MQ=255)
                                                       tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              <  1:496454/149‑1 (MQ=255)
                                                        gCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGata                                                                                             <  3:431614/149‑1 (MQ=255)
                                                         ccGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGatat                                                                                            <  2:143700/149‑1 (MQ=255)
                                                         ccGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGatat                                                                                            >  2:687071/1‑149 (MQ=255)
                                                            tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGat                                                                                              <  1:906026/144‑1 (MQ=255)
                                                            tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         >  4:530739/1‑149 (MQ=255)
                                                            tttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATACCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         >  1:1264032/1‑149 (MQ=255)
                                                              tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       <  4:1149896/149‑1 (MQ=255)
                                                              tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       <  1:237565/149‑1 (MQ=255)
                                                                ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTc                                                                                         <  3:532169/145‑1 (MQ=255)
                                                                ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       <  2:1307276/147‑1 (MQ=255)
                                                                ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCt                                                                                       >  1:1307276/1‑147 (MQ=255)
                                                                ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTgg                                                                                     >  4:447972/1‑149 (MQ=255)
                                                                     tCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTAt                                                                                <  1:1459058/149‑1 (MQ=255)
                                                                      cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                                                               <  1:687071/149‑1 (MQ=255)
                                                                      cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                                                               <  4:1207301/149‑1 (MQ=255)
                                                                       tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                                                                                                                                     >  1:1573664/1‑78 (MQ=255)
                                                                       tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGga                                                                                                                                                     <  2:1573664/78‑1 (MQ=255)
                                                                       tGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                                                              >  3:1118017/1‑149 (MQ=255)
                                                                         cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  2:138544/149‑1 (MQ=255)
                                                                          gATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAg                                                                           >  2:432984/1‑149 (MQ=255)
                                                                           aTAGTCACAAAGCTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGt                                                                          <  4:76443/149‑1 (MQ=255)
                                                                               tCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTAtt                                                                      >  2:1150401/1‑149 (MQ=255)
                                                                                cacaAAGGTAATAGTTGAAATTCCCCTGCCACCCGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                                                                 <  1:1194820/73‑1 (MQ=255)
                                                                                cacaAAGGTAATAGTTGAAATTCCCCTGCCACCCGGCAAAATATCCGTTCAACCATCAGCTTTGCAGgacgac                                                                                                                                                 >  2:1194820/1‑73 (MQ=255)
                                                                                 acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  4:748992/1‑149 (MQ=255)
                                                                                 acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  1:1496953/1‑149 (MQ=255)
                                                                                 acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  3:640433/1‑149 (MQ=255)
                                                                                 acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    <  3:1530392/149‑1 (MQ=255)
                                                                                 acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                                                                    >  3:1202611/1‑149 (MQ=255)
                                                                                  caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                                                                   <  2:243526/149‑1 (MQ=255)
                                                                                  caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCt                                                                   <  4:803257/149‑1 (MQ=255)
                                                                                   aaaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAATATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTg                                                                  <  2:1014630/149‑1 (MQ=255)
                                                                                       gTAATAGTTGAAATTCCTCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGa                                                                                                     <  2:1289828/110‑1 (MQ=255)
                                                                                       gTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTct                                                              <  2:1179988/149‑1 (MQ=255)
                                                                                         aaTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                                                                            <  3:270800/133‑1 (MQ=255)
                                                                                            aGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTGTTGCTGGTCTGCTgg                                                         <  3:1479655/149‑1 (MQ=255)
                                                                                              ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa                                                       >  3:685900/1‑149 (MQ=255)
                                                                                               tGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAg                                                      <  4:221230/149‑1 (MQ=255)
                                                                                                 aaaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGc                                                                                                   >  2:1066196/1‑102 (MQ=255)
                                                                                                 aaaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg                                                    <  2:480989/149‑1 (MQ=255)
                                                                                                  aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCGAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTgg                                                   <  2:147658/149‑1 (MQ=255)
                                                                                                  aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTgg                                                   >  2:390976/1‑149 (MQ=255)
                                                                                                  aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTgg                                                   <  3:997062/149‑1 (MQ=255)
                                                                                                   aTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGt                                                  >  1:48347/1‑149 (MQ=255)
                                                                                                   aTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGt                                                  <  4:343878/149‑1 (MQ=255)
                                                                                                    ttCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTa                                                 <  2:170215/149‑1 (MQ=255)
                                                                                                    ttCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTa                                                 <  1:214732/149‑1 (MQ=255)
                                                                                                        ccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                               <  4:283217/99‑1 (MQ=255)
                                                                                                        ccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                                                               >  3:283217/1‑99 (MQ=255)
                                                                                                             cACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGtt                                        <  2:767332/149‑1 (MQ=255)
                                                                                                                cTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGc                                     >  1:658202/1‑149 (MQ=255)
                                                                                                                  ggCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCt                                   <  3:86399/149‑1 (MQ=255)
                                                                                                                     aaaaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGaa                                >  2:66858/1‑149 (MQ=255)
                                                                                                                      aaaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGaaa                               <  1:1150401/149‑1 (MQ=255)
                                                                                                                          atCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAAtata                           >  3:566914/1‑149 (MQ=255)
                                                                                                                               ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCAc                      >  4:1465992/1‑149 (MQ=255)
                                                                                                                                tCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAAtata                           <  2:968616/143‑1 (MQ=255)
                                                                                                                                tCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAAtata                           >  1:968616/1‑143 (MQ=255)
                                                                                                                                     cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCttt                >  1:41624/1‑149 (MQ=255)
                                                                                                                                     cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCttt                >  2:521895/1‑149 (MQ=255)
                                                                                                                                         aGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTc            >  4:1029651/1‑149 (MQ=255)
                                                                                                                                         aGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTc            >  2:1415198/1‑149 (MQ=255)
                                                                                                                                           gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTg                                       <  1:84280/119‑1 (MQ=255)
                                                                                                                                           gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt          >  1:1538269/2‑149 (MQ=255)
                                                                                                                                           gtttGCAGGTCGTCCTGCAAACGCCTCGTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAGTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt          >  2:553945/2‑149 (MQ=255)
                                                                                                                                           cTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGt          <  1:142454/149‑1 (MQ=255)
                                                                                                                                            tttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGtt         >  4:82749/1‑149 (MQ=255)
                                                                                                                                             ttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGttt        <  4:1083972/149‑1 (MQ=255)
                                                                                                                                              tGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTg       >  3:482951/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGCGCAGCAACGAAAGCGCAGCAATAAGCAATAACGGTACGACAGCTGTGTCGTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCC  >  NC_000913/4376671‑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: