Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,643,428 C→G 5.3% S373C (TCC→TGC)  PP_4951 → EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,643,4280CG 5.3% 510.2 / 13.5 152S373C (TCC→TGC) PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (85/59);  new base G (5/3);  total (90/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 156 total aligned reads displayed.

ACCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTACTG  >  NC_002947/5643279‑5643577
                                                                                                                                                     |                                                                                                                                                     
aCCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTc                                                                                                                                                       >  2:1267227/1‑150 (MQ=255)
aCCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTc                                                                                                                                                       <  1:963874/150‑1 (MQ=255)
   gCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCAt                                                                                                                                                    >  1:2183470/1‑150 (MQ=255)
    cTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATg                                                                                                                                                   <  1:894336/150‑1 (MQ=255)
      gATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCa                                                                                                                                                 <  1:2613110/150‑1 (MQ=255)
         cGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGtt                                                                                                                                              >  1:1463739/1‑150 (MQ=255)
         cGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGtt                                                                                                                                              >  1:2359068/1‑150 (MQ=255)
          gAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTc                                                                                                                                             >  2:3608191/1‑150 (MQ=255)
           aaCGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCa                                                                                                                                            >  2:4578162/1‑150 (MQ=255)
            aCGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCAc                                                                                                                                           >  2:110258/1‑150 (MQ=255)
               gTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTGCATGGAGTTCACTGc                                                                                                                                        <  1:4458751/150‑1 (MQ=255)
                  cGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGAt                                                                                                                                     <  2:5083772/150‑1 (MQ=255)
                  cGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGAt                                                                                                                                     <  2:1876062/150‑1 (MQ=255)
                   ggCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATg                                                                                                                                    <  2:1593616/150‑1 (MQ=255)
                     cTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCgaagaa                                                                                                                                  <  1:2384958/150‑1 (MQ=255)
                     cTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCgaagaa                                                                                                                                  <  2:148621/150‑1 (MQ=255)
                     cTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGaa                                                                                                                                  <  2:4320698/150‑1 (MQ=255)
                        tCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAacac                                                                                                                               >  1:2855417/1‑150 (MQ=255)
                         cGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCGGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAAcaca                                                                                                                              <  1:4142488/150‑1 (MQ=255)
                               atcCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACg                                                                                                                        >  1:1949842/1‑150 (MQ=255)
                                  cTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTc                                                                                                                     >  2:1598386/1‑150 (MQ=255)
                                   tCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCa                                                                                                                    >  1:2447494/1‑150 (MQ=255)
                                   tCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCa                                                                                                                    >  2:960783/1‑150 (MQ=255)
                                         gTAATACCTTGCTGTTGGGCATTCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtc                                                                                                              <  2:3887572/150‑1 (MQ=255)
                                           aaTACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGAtctc                                                                                                            <  1:4403596/150‑1 (MQ=255)
                                             tACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcga                                                                                                          <  1:1352296/150‑1 (MQ=255)
                                              aCCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTcgac                                                                                                         <  1:4561276/150‑1 (MQ=255)
                                                 ttGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGAGGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCgacgac                                                                                                      >  2:2826549/1‑150 (MQ=255)
                                                 ttGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCgacgac                                                                                                      >  2:5027738/1‑150 (MQ=255)
                                                     tGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                  >  1:197280/1‑150 (MQ=255)
                                                     tGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgcca                                                                                                  <  1:1877555/150‑1 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccagc                                                                                                >  1:3334825/1‑150 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccagc                                                                                                >  2:215089/1‑150 (MQ=255)
                                                       ttGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccagc                                                                                                <  2:2438790/150‑1 (MQ=255)
                                                         gggCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGAcgccagcgc                                                                                              >  2:4227681/1‑150 (MQ=255)
                                                              tCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCcagca                                                                                         >  1:2772568/1‑150 (MQ=255)
                                                                 tgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                      >  2:773896/1‑150 (MQ=255)
                                                                 tgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                      >  1:1192142/1‑150 (MQ=255)
                                                                 tgctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCGTATTGGCTGCTGTTGGTGCAACTTGGCCTGCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGCTGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCg                                                                                      >  2:2989888/1‑150 (MQ=255)
                                                                  gctgGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCgg                                                                                     <  1:446774/150‑1 (MQ=255)
                                                                     ggCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTaa                                                                                  >  2:3810098/1‑150 (MQ=255)
                                                                      gCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAAc                                                                                 <  2:2817947/150‑1 (MQ=255)
                                                                       cAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACa                                                                                <  1:4659551/150‑1 (MQ=255)
                                                                          gcctgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGcc                                                                             >  1:2002777/1‑150 (MQ=255)
                                                                          gcctgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTAGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGcc                                                                             <  2:2312587/150‑1 (MQ=255)
                                                                            ctgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctg                                                                           <  2:4936097/150‑1 (MQ=255)
                                                                             tgggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctgc                                                                          >  2:2989543/1‑150 (MQ=255)
                                                                               ggcctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctgctg                                                                        >  1:4300557/1‑150 (MQ=255)
                                                                                 cctggTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTc                                                                      >  1:1645231/1‑150 (MQ=255)
                                                                                   tggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCgg                                                                    >  2:4271987/1‑150 (MQ=255)
                                                                                    ggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCggt                                                                   <  2:4803650/150‑1 (MQ=255)
                                                                                    ggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCggt                                                                   <  2:4803502/150‑1 (MQ=255)
                                                                                    ggtggACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCggt                                                                   <  2:4571937/150‑1 (MQ=255)
                                                                                      tggACGAGCAGACCCCTTCTTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCggtgg                                                                 <  2:4465246/150‑1 (MQ=255)
                                                                                        gACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTc                                                               <  2:1645231/150‑1 (MQ=255)
                                                                                             gCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCa                                                          >  2:2190932/1‑150 (MQ=255)
                                                                                                 aCCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTg                                                      <  1:2942502/150‑1 (MQ=255)
                                                                                                  ccccTTATTGCCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGt                                                     <  2:2067204/150‑1 (MQ=255)
                                                                                                   cccTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTc                                                    >  2:2045747/1‑150 (MQ=255)
                                                                                                     cTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGc                                                  >  1:2023369/1‑150 (MQ=255)
                                                                                                     cTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGc                                                  >  1:1520911/1‑150 (MQ=255)
                                                                                                      ttATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCt                                                 <  2:3373117/150‑1 (MQ=255)
                                                                                                      ttATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCt                                                 <  2:3373122/150‑1 (MQ=255)
                                                                                                      ttATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCt                                                 >  1:4310774/1‑150 (MQ=255)
                                                                                                        aTTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                               <  2:798372/150‑1 (MQ=255)
                                                                                                        aTTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                               <  2:798510/150‑1 (MQ=255)
                                                                                                        aTTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGa                                               <  1:3288494/150‑1 (MQ=255)
                                                                                                            gctgctGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCt                                           >  1:4040176/1‑150 (MQ=255)
                                                                                                              tgctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                         >  2:3215687/1‑150 (MQ=255)
                                                                                                              tgctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                         <  1:2746334/150‑1 (MQ=255)
                                                                                                              tgctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                         <  1:4648913/150‑1 (MQ=255)
                                                                                                              tgctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTgg                                         >  1:756167/1‑150 (MQ=255)
                                                                                                                ctgTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGgt                                       <  2:4464629/150‑1 (MQ=255)
                                                                                                                   ttGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGtgtg                                    >  2:2905070/1‑150 (MQ=255)
                                                                                                                   ttGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGtgtg                                    >  1:1081267/1‑150 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGcc                                 >  2:3872049/1‑150 (MQ=255)
                                                                                                                      gTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGcc                                 >  2:3872115/1‑150 (MQ=255)
                                                                                                                             ttGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCg                          >  1:2657939/1‑150 (MQ=255)
                                                                                                                             ttGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCg                          >  1:2657988/1‑150 (MQ=255)
                                                                                                                             ttGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCg                          >  1:516612/1‑150 (MQ=255)
                                                                                                                             ttGGCCTTCTGGGCGCAGTGAACTCCATGCAGATCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCg                          >  1:516919/1‑150 (MQ=255)
                                                                                                                              tGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCgg                         >  2:2182451/1‑150 (MQ=255)
                                                                                                                                 ccTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGgcgc                      >  1:1679209/1‑150 (MQ=255)
                                                                                                                                  cTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACGGTCGCTGAGCCTGGGGGTGGCcggcgccggcgcc                     >  2:2282541/1‑149 (MQ=255)
                                                                                                                                   ttCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCAc                    >  2:4667967/1‑150 (MQ=255)
                                                                                                                                   ttCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCAc                    >  2:4667968/1‑150 (MQ=255)
                                                                                                                                      tGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCCGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGt                 >  1:4362218/1‑150 (MQ=255)
                                                                                                                                       gggCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGtt                <  2:2400724/150‑1 (MQ=255)
                                                                                                                                         tcgcAGTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAAGTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTgg              >  2:4041287/2‑150 (MQ=255)
                                                                                                                                         gcgcAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTgg              >  1:778304/1‑150 (MQ=255)
                                                                                                                                         gcgcAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTgg              >  2:696586/1‑150 (MQ=255)
                                                                                                                                         gcgcAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTCTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTgg              >  2:186662/1‑150 (MQ=255)
                                                                                                                                             aGTGAACTGCATGGAGTTCACTGCGCTGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcgg          >  2:1881646/1‑150 (MQ=255)
                                                                                                                                             aGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcgg          >  2:2588122/1‑150 (MQ=255)
                                                                                                                                             aGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcgg          >  2:2588100/1‑150 (MQ=255)
                                                                                                                                             aGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGggcgg          <  2:3796236/150‑1 (MQ=255)
                                                                                                                                              gTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctg                                                                           <  2:3086911/84‑1 (MQ=255)
                                                                                                                                              gTGAACTGCATGGAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCctg                                                                           >  1:3086911/1‑84 (MQ=255)
                                                                                                                                                gAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCtt       <  1:110258/150‑1 (MQ=255)
                                                                                                                                                   cTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTAc    <  2:1856698/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCTGATCGAACGGTTCGGCTATCGCATCATCCTCACCGGTAATACCTTGCTGTTGGGCATCCTGCTGGCAAGCCTGGGCCTGGTGGACGAGCAGACCCCTTATTGGCTGCTGTTGGTGCAACTTGGCCTTCTGGGCGCAGTGAACTCCATGCAGTTCACTGCGATGAACACAGTGACGCTCATCGATCTCGACGACGCCAGCGCCAGCAGCGGTAACAGCCTGCTGTCGGTGGTCGCGCAACTGTCGCTGAGCCTGGGTGTGGCCTGCGCCGGCGCACTGTTGGGCGGCTTTACTG  >  NC_002947/5643279‑5643577

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