Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,046,233 C→G 5.2% S716S (TCG→TCC PP_0906 ← RND family multidrug transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,046,2330CG 5.2% 312.4 / 6.8 97S716S (TCG→TCCPP_0906RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base C (51/41);  new base G (2/3);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

TTCACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCGGGTA  >  NC_002947/1046084‑1046382
                                                                                                                                                       |                                                                                                                                                     
ttcACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATc                                                                                                                                                       <  2:2423211/150‑1 (MQ=255)
   aCGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGac                                                                                                                                                    <  1:2129768/150‑1 (MQ=255)
   aCGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGac                                                                                                                                                    >  2:60458/1‑150 (MQ=255)
        gtagtTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACAcccc                                                                                                                                               >  1:2009065/1‑150 (MQ=255)
         tagtTGTTCAGCCAGCCGGGGTTGTCGCGATATGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGGCAccctc                                                                                                                                              <  1:1704340/150‑3 (MQ=255)
          agtTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCa                                                                                                                                             <  2:1653784/150‑1 (MQ=255)
            ttgttCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATc                                                                                                                                           >  1:2436417/1‑150 (MQ=255)
             tgttCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGACCCCCTTAGCACCGGCGCCTTCGAG‑‑TGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCt                                                                                                                                          <  1:2113962/150‑1 (MQ=255)
               ttCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCCCGCCCAGCAGGGTCGCCAAGGTGGCGCCCAGGGTGTCCATCGACACCCCCAGCTgc                                                                                                                                        >  1:275645/1‑150 (MQ=255)
                   gccagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGGGCCGCCGAGGGTGTCCATCGACACCCCCATCTGcgccg                                                                                                                                    >  1:947672/1‑150 (MQ=255)
                    ccagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCGTCGCCACCCCCCTCTGcgccgc                                                                                                                                   >  2:1959529/1‑150 (MQ=255)
                    ccagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGcgccgc                                                                                                                                   >  1:978556/1‑150 (MQ=255)
                     cagccGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCgccgcc                                                                                                                                  >  1:1961011/1‑150 (MQ=255)
                        ccGGGGTTGTCGCGGTAATGTCGTTCAACCTGGGCGATCACCTTGTAACTGCCGCGTTCCAGAAGGTGACCGGGTATTTCAGCCTCGCCGAGCAGGGTCGTCATGATGCATACGTGTGTTGCCATCGACACCCCCATCTGCGCCGCCttt                                                                                                                               <  1:1125278/150‑1 (MQ=255)
                           gggTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGcc                                                                                                                            >  1:518602/1‑150 (MQ=255)
                                gTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTc                                                                                                                       <  2:1961011/150‑1 (MQ=255)
                                gTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTc                                                                                                                       <  1:1959529/150‑1 (MQ=255)
                                 tCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCg                                                                                                                      >  1:2172394/1‑150 (MQ=255)
                                   gcgGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGAt                                                                                                                    <  2:2009065/150‑1 (MQ=255)
                                          ggCCGTTCAACCTGGGCGATCACAAAGACCTTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccacc                                                                                                             >  1:1242414/1‑150 (MQ=255)
                                           gCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccacca                                                                                                            >  2:1122301/1‑150 (MQ=255)
                                            ccGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccaccac                                                                                                           >  2:552558/1‑150 (MQ=255)
                                              gTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTccaccacca                                                                                                         >  1:501917/1‑150 (MQ=255)
                                                  aaCCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTc                                                                                                     >  1:571324/1‑150 (MQ=255)
                                                    ccTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCAAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCgg                                                                                                   >  2:1146968/1‑150 (MQ=255)
                                                     cTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCggg                                                                                                  <  2:615776/150‑1 (MQ=255)
                                                        ggCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGttt                                                                                               <  1:1792247/150‑1 (MQ=255)
                                                           gATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTc                                                                                            <  2:2436417/150‑1 (MQ=255)
                                                           gATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTc                                                                                            >  1:1099789/1‑150 (MQ=255)
                                                              cacaTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGaa                                                                                         >  2:1017929/1‑150 (MQ=255)
                                                                 cTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctc                                                                                                                                              >  1:2384900/1‑92 (MQ=255)
                                                                 cTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCGATGGACAccctc                                                                                                                                              <  2:2384900/94‑3 (MQ=255)
                                                                 cTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCCGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGc                                                                                      >  1:1562382/1‑150 (MQ=255)
                                                                  ttGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCa                                                                                     >  1:103587/1‑150 (MQ=255)
                                                                         tGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGa                                                                              >  1:103008/1‑150 (MQ=255)
                                                                            ggCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCgg                                                                                                   <  2:708510/126‑1 (MQ=255)
                                                                            ggCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCgg                                                                                                   >  1:708510/1‑126 (MQ=255)
                                                                            ggCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGatat                                                                           <  2:1099789/150‑1 (MQ=255)
                                                                            ggCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTGGGCGAAGGCAAGATCGatat                                                                           >  2:1762392/1‑150 (MQ=255)
                                                                               cTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCTGGGCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGa                                                                        <  1:60458/150‑1 (MQ=255)
                                                                               cTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGa                                                                        <  2:501917/150‑1 (MQ=255)
                                                                                ttCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAg                                                                       >  2:1528796/1‑150 (MQ=255)
                                                                                ttCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAg                                                                       <  2:103008/150‑1 (MQ=255)
                                                                                ttCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAg                                                                       >  2:2259570/1‑150 (MQ=255)
                                                                                ttCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCACCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAg                                                                       >  2:1884260/1‑150 (MQ=255)
                                                                                  cAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGa                                                                     <  1:873640/150‑1 (MQ=255)
                                                                                    aG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAg                                                                   >  1:858827/1‑150 (MQ=255)
                                                                                     g‑‑ggTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAgg                                                                  <  2:879318/150‑1 (MQ=255)
                                                                                        ggTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGc                                                                 >  1:2046997/1‑150 (MQ=255)
                                                                                          tGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGa                                                               >  2:720145/1‑150 (MQ=255)
                                                                                            aaaCGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAAc                                                             <  1:1122301/150‑1 (MQ=255)
                                                                                               cGGTTGATTTACGCCTCGCCAAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTg                                                          >  1:1066276/1‑150 (MQ=255)
                                                                                                 gTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGata                                                                            <  1:1300022/130‑1 (MQ=255)
                                                                                                 gTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGata                                                                            >  2:1300022/1‑130 (MQ=255)
                                                                                                 gTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGcc                                                        >  1:2452013/1‑150 (MQ=255)
                                                                                                  ttGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGccc                                                       <  2:571324/150‑1 (MQ=255)
                                                                                                   tGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGTCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCg                                                      >  2:2419993/1‑150 (MQ=255)
                                                                                                   tGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGaa                                                              <  2:910835/142‑1 (MQ=255)
                                                                                                   tGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGaa                                                              <  2:2149366/142‑1 (MQ=255)
                                                                                                   tGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGaa                                                              >  1:910835/1‑142 (MQ=255)
                                                                                                   tGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGaa                                                              >  1:2149366/1‑142 (MQ=255)
                                                                                                          cGCCTCGCCCAGCAGGGTCGCCAAGGTGTCGCCGATGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTGGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTg                                               >  2:299083/1‑150 (MQ=255)
                                                                                                          cGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTg                                               >  2:1899314/1‑150 (MQ=255)
                                                                                                           gCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTccgta                                              <  1:91250/150‑4 (MQ=255)
                                                                                                           gCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGa                                              <  1:1132386/150‑1 (MQ=255)
                                                                                                            ccTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAg                                             >  2:775769/1‑150 (MQ=255)
                                                                                                               cGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAgcgc                                          <  2:1066276/150‑1 (MQ=255)
                                                                                                                 cccAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAgcgcgc                                        <  2:2046997/150‑1 (MQ=255)
                                                                                                                 cccAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCCGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAgcgcgc                                        <  2:1562382/150‑1 (MQ=255)
                                                                                                                     gcagGGTCGCCAAGGTGTCGCCGATGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGCAAGCAAGATCGATATCGAGGAAGGCGAACTTTCCCGATTCCTGAGCGCGCGTCt                                    <  1:299083/150‑1 (MQ=255)
                                                                                                                      cagGGTCGCCAAGGTGCCGCCTAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCtt                                   <  2:589483/150‑1 (MQ=255)
                                                                                                                         ggTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGcc                                                        >  1:591457/1‑126 (MQ=255)
                                                                                                                         ggTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGcc                                                        <  2:591457/126‑1 (MQ=255)
                                                                                                                         ggTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGAt                                <  2:518602/150‑1 (MQ=255)
                                                                                                                              ccAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGCTGCGtt                           >  1:2249938/1‑150 (MQ=255)
                                                                                                                               cAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTg                          >  2:1936694/1‑150 (MQ=255)
                                                                                                                                aaGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCTCGTCTTGATGCGTTgg                         <  1:1146968/150‑1 (MQ=255)
                                                                                                                                aaGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         <  2:1551862/150‑1 (MQ=255)
                                                                                                                                aaGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         <  1:1899314/150‑1 (MQ=255)
                                                                                                                                        gccgAGGGTGTCCATCGACACCCCCACCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTc                 <  1:1884260/150‑1 (MQ=255)
                                                                                                                                         ccgAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTcc                >  2:1606196/1‑150 (MQ=255)
                                                                                                                                             gggTGTCGATGGACACCCTCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGccc                                                       <  1:1894665/107‑1 (MQ=255)
                                                                                                                                             gggTGTCGATGGACACCCTCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGccc                                                       >  2:1894665/1‑107 (MQ=255)
                                                                                                                                             gggTGTCAATCGACACCCCCTTCTGCGCCGCCTGTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCcagca            <  2:275645/150‑1 (MQ=255)
                                                                                                                                              ggTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCcagcag           >  1:1931200/1‑150 (MQ=255)
                                                                                                                                              ggTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGAAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCcagcag           >  2:652194/1‑150 (MQ=255)
                                                                                                                                               gtgtCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGt          >  1:2404116/1‑150 (MQ=255)
                                                                                                                                               gtgtCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGt          <  1:1501594/150‑1 (MQ=255)
                                                                                                                                                   ccATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCg      >  2:1724175/1‑150 (MQ=255)
                                                                                                                                                   ccATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCg      >  1:219070/1‑150 (MQ=255)
                                                                                                                                                    cATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGTTTCCGGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCgg     <  1:552558/150‑1 (MQ=255)
                                                                                                                                                    cATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         >  2:552878/1‑130 (MQ=255)
                                                                                                                                                    cATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         <  1:552878/130‑1 (MQ=255)
                                                                                                                                                    cATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCgg     <  1:775769/150‑1 (MQ=255)
                                                                                                                                                     aTCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         >  1:154806/1‑129 (MQ=255)
                                                                                                                                                     aTCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTgg                         <  2:154806/129‑1 (MQ=255)
                                                                                                                                                     aTCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCggg    >  2:1759600/1‑150 (MQ=255)
                                                                                                                                                     aTCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCggg    >  2:1484201/1‑150 (MQ=255)
                                                                                                                                                       cGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCGGGTa  <  1:794020/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TTCACGTAGTAGTTGTTCAGCCAGCCGGGGTTGTCGCGGTAAGGCCGTTCAACCTGGGCGATCACCTTGTAACTGCGGCCTTCAAG‑‑GGTGAAACGGTTGATTTCCGCCTCGCCCAGCAGGGTCGCCAAGGTGCCGCCGAGGGTGTCCATCGACACCCCCATCTGCGCCGCCTTTGCCCGGTCGATGTCCACCACCACTTCGGGTTTGTCGAAGGCAAGATCGATATCGAGGAAGGCGAACTTGCCCGATTCCTGAGCGCGCGTCTTGATGCGTTGGGCCACGTCCAGCAGTGCCGGGTA  >  NC_002947/1046084‑1046382

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