Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,227,441 A→G pseudogene (1318/1436 nt) mdtQ ← putative multidrug resistance outer membrane protein MdtQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,227,4410AG98.9% 320.5 / ‑4.7 90pseudogene (1318/1436 nt)mdtQputative multidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (62/27);  total (63/27)
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.58e-01

GAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCCGCGCCA  >  NC_000913/2227309‑2227576
                                                                                                                                    |                                                                                                                                       
gAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGa                                                                                                                                  >  1:500247/1‑140 (MQ=255)
gAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGa                                                                                                                                  <  2:397897/140‑1 (MQ=255)
 aaGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGaa                                                                                                                                 >  2:321814/1‑140 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAGTGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCa                                                                                                                              >  1:1194842/1‑140 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCTTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCa                                                                                                                              <  1:1472677/140‑1 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCa                                                                                                                              >  1:791671/1‑140 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCa                                                                                                                              >  2:200141/1‑140 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCa                                                                                                                              >  1:480052/1‑140 (MQ=255)
    gcgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAg                                                                                                                             <  1:174440/140‑1 (MQ=255)
     cgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAg                                                                                                                             >  1:1313373/1‑140 (MQ=255)
     cgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAg                                                                                                                             >  1:424194/1‑140 (MQ=255)
     cgcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAg                                                                                                                             >  2:1091693/1‑140 (MQ=255)
      gcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGc                                                                                                                            >  2:907059/1‑140 (MQ=255)
      gcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGc                                                                                                                            >  1:1324626/1‑140 (MQ=255)
      gcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGc                                                                                                                            >  1:940405/1‑140 (MQ=255)
      gcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGc                                                                                                                            >  2:462148/1‑140 (MQ=255)
      gcATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGc                                                                                                                            >  2:1115941/1‑140 (MQ=255)
       cATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCa                                                                                                                           >  1:1529882/1‑140 (MQ=255)
        atatCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCaa                                                                                                                          >  1:457596/1‑140 (MQ=255)
          atcatcAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAatg                                                                                                                        <  2:1220869/140‑1 (MQ=255)
                    tttGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGt                                                                                                              >  2:461358/1‑140 (MQ=255)
                    tttGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGt                                                                                                              >  2:1168231/1‑140 (MQ=255)
                      tGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTa                                                                                                            >  2:1523543/1‑140 (MQ=255)
                      tGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTa                                                                                                            >  1:1103680/1‑140 (MQ=255)
                         aCCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAg                                                                                                         <  1:1382492/140‑1 (MQ=255)
                            ccccGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgc                                                                                                      >  2:11462/1‑140 (MQ=255)
                             cccGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcg                                                                                                     >  2:163378/1‑140 (MQ=255)
                              ccGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcgc                                                                                                    <  1:854388/140‑1 (MQ=255)
                              ccGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcgc                                                                                                    <  2:400711/140‑1 (MQ=255)
                              ccGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAgcgcgc                                                                                                    <  2:480052/140‑1 (MQ=255)
                               cGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGcc                                                                                                   >  1:1446685/1‑140 (MQ=255)
                               cGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGcc                                                                                                   >  1:95710/1‑140 (MQ=255)
                               cGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGcc                                                                                                   >  2:1541807/1‑140 (MQ=255)
                               cGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGcc                                                                                                   >  1:347044/1‑140 (MQ=255)
                                    aaCGCACCAGTGAGTTGAATGGAGGCATCCAGCCCCTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTgcgc                                                                                              <  1:1549175/140‑1 (MQ=255)
                                      cGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCaa                                                                                            <  1:1052812/140‑1 (MQ=255)
                                        cACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGa                                                                                          >  1:1129468/1‑140 (MQ=255)
                                              tGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTAcc                                                                                    <  2:940405/140‑1 (MQ=255)
                                                aGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACcac                                                                                  <  2:1274615/140‑1 (MQ=255)
                                                 gTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACcaca                                                                                 >  2:337136/1‑140 (MQ=255)
                                                   tGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACg                                                                               >  1:335672/1‑140 (MQ=255)
                                                     aaTGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCa                                                                             >  2:1265377/1‑140 (MQ=255)
                                                     aaTGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCa                                                                             >  2:217214/1‑140 (MQ=255)
                                                     aaTGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCa                                                                             >  2:624466/1‑140 (MQ=255)
                                                       tGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAg                                                                           >  2:656427/1‑140 (MQ=255)
                                                       tGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAg                                                                           <  2:1324626/140‑1 (MQ=255)
                                                            gCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGTGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATc                                                                      <  2:1423692/140‑1 (MQ=255)
                                                            gCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATc                                                                      >  1:1361628/1‑140 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   >  2:1280881/1‑140 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   >  2:583629/1‑140 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   <  2:580417/140‑1 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   >  2:491280/1‑140 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   >  1:1358651/1‑140 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   <  1:907059/140‑1 (MQ=255)
                                                               tCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATcgcg                                                                   >  2:676611/1‑140 (MQ=255)
                                                                         gCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTgc                                                         >  1:761973/1‑140 (MQ=255)
                                                                          cccTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTgcg                                                        >  1:677920/1‑140 (MQ=255)
                                                                          cccTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTgcg                                                        >  2:1202187/1‑140 (MQ=255)
                                                                               gCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGa                                                   >  2:20862/1‑140 (MQ=255)
                                                                                cAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGAt                                                  <  2:424194/140‑1 (MQ=255)
                                                                                    aaCAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTg                                              >  1:1250535/1‑140 (MQ=255)
                                                                                     aCAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTgg                                             >  1:638810/1‑140 (MQ=255)
                                                                                        ggCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGttt                                          <  1:163378/140‑1 (MQ=255)
                                                                                        ggCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGttt                                          >  1:209939/1‑140 (MQ=255)
                                                                                        ggCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGttt                                          <  1:1523543/140‑1 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  2:108455/1‑140 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  1:797813/1‑140 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  2:850180/1‑140 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  1:301675/1‑140 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  1:1490847/1‑140 (MQ=255)
                                                                                         gCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGtttt                                         >  2:750722/1‑140 (MQ=255)
                                                                                          ccATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGttttt                                        >  1:1417070/1‑140 (MQ=255)
                                                                                                 cccGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGcc                                 <  1:1091693/140‑1 (MQ=255)
                                                                                                  ccGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCACTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCa                                <  1:321814/140‑1 (MQ=255)
                                                                                                   cGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAg                               >  2:947591/1‑140 (MQ=255)
                                                                                                     cTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAgtg                             >  2:1461676/1‑140 (MQ=255)
                                                                                                     cTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAgtg                             >  2:1260228/1‑140 (MQ=255)
                                                                                                     cTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAgtg                             >  2:904693/1‑140 (MQ=255)
                                                                                                     cTCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAggg                             >  2:1531048/1‑138 (MQ=255)
                                                                                                      tCACGCAGCGCGGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAgtgt                            <  1:656427/140‑1 (MQ=255)
                                                                                                              cgcgGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCt                    >  2:558028/1‑140 (MQ=255)
                                                                                                              cgcgGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCt                    >  1:439522/1‑140 (MQ=255)
                                                                                                              cgcgGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCt                    <  1:1541807/140‑1 (MQ=255)
                                                                                                                cgGGGATTCTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGa                  >  2:1188236/1‑140 (MQ=255)
                                                                                                                        cTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGc          >  2:738033/1‑140 (MQ=255)
                                                                                                                        cTGGCTTCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGc          <  1:1265377/140‑1 (MQ=255)
                                                                                                                        cTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGc          >  1:1080641/1‑140 (MQ=255)
                                                                                                                             ttCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCcgcg     <  2:600084/140‑1 (MQ=255)
                                                                                                                             ttCGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCcgcg     <  1:1115941/140‑1 (MQ=255)
                                                                                                                               cGCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCCGCGcc   <  1:1280881/140‑1 (MQ=255)
                                                                                                                                gCTGGCGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCCGCGCCa  <  2:761973/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GAAGGCGCATATCATCAGCGTTTGTACCCCCCGCCCAACGCACCAGTGAGTTGAATGGAGGCATCCAGCCACTGCCCTTGCAATAACAGGCCATTGGCCCGCTCACGCAGCGCGGGGATTCTGGCTTCGCTGACGCGGGAACCAGCAATGATGCCCGCGTTAAAGCGCGCCTGCGCAAGACCTACCACACGCAAGGCATCGCGCTCAATTTGCGCCTGATGCTGGTTTTTCTCCGCCAGTGTCTGAACCTGACTGGCTGCCCGCGCCA  >  NC_000913/2227309‑2227576

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