Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,143,406 C→A 6.5% R173S (CGC→AGC)  yadG → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,143,4060CA 6.5% 1035.4 / 43.0 308R173S (CGC→AGC) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (135/153);  new base A (12/8);  total (147/161)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

AGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATGCGCCA  >  NC_002947/2143257‑2143554
                                                                                                                                                     |                                                                                                                                                    
aGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGcgcc                                                                                                                                                      >  2:3101261/1‑150 (MQ=255)
 gCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGcgccg                                                                                                                                                     >  2:6337136/1‑150 (MQ=255)
  cTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGcgccgc                                                                                                                                                    >  1:6721695/1‑150 (MQ=255)
  cTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGcgccgc                                                                                                                                                    >  2:2367751/1‑150 (MQ=255)
   tGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCt                                                                                                                                                   >  2:1201618/1‑150 (MQ=255)
   tGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCt                                                                                                                                                   >  1:7240573/1‑150 (MQ=255)
    gggCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTc                                                                                                                                                  <  2:7090202/150‑1 (MQ=255)
    gggCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTc                                                                                                                                                  <  2:5339063/150‑1 (MQ=255)
      gCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGa                                                                                                                                                <  2:839169/150‑1 (MQ=255)
       ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                               <  1:4825584/150‑1 (MQ=255)
       ccTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAt                                                                                                                                               <  1:3750020/150‑1 (MQ=255)
        cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtg                                                                                                                                              >  1:3891406/1‑150 (MQ=255)
        cTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCGGCGCAGCTCGAta                                                                                                                                              <  1:6584406/150‑2 (MQ=255)
          gtgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGAtgtg                                                                                                                                            <  1:697100/150‑1 (MQ=255)
           tgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTgg                                                                                                                                           >  1:2331431/1‑150 (MQ=255)
           tgGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTgg                                                                                                                                           <  2:8171027/150‑1 (MQ=255)
            gggACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGa                                                                                                                                          <  1:5376030/150‑1 (MQ=255)
                 aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTc                                                                                                                                     >  1:8073317/1‑150 (MQ=255)
                 aaaCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTc                                                                                                                                     <  2:6406464/150‑1 (MQ=255)
                   aCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCt                                                                                                                                   >  2:513345/1‑150 (MQ=255)
                   aCGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCt                                                                                                                                   >  1:4688285/1‑150 (MQ=255)
                     gTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCa                                                                                                                                 <  1:120980/150‑1 (MQ=255)
                         tgtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGCTGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGa                                                                                                                             <  2:8560803/150‑1 (MQ=255)
                         tgtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGa                                                                                                                             <  1:916109/150‑1 (MQ=255)
                         tgtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGa                                                                                                                             <  1:6793282/150‑1 (MQ=255)
                          gtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAg                                                                                                                            <  2:4468424/150‑1 (MQ=255)
                          gtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAg                                                                                                                            <  1:5816037/150‑1 (MQ=255)
                          gtgCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAg                                                                                                                            >  1:410274/1‑150 (MQ=255)
                              aGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGa                                                                                                                        >  1:7822441/1‑150 (MQ=255)
                               gTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGaa                                                                                                                       <  2:1078397/150‑1 (MQ=255)
                               gTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGaa                                                                                                                       <  1:5561489/150‑1 (MQ=255)
                                  gcgTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCa                                                                                                                    >  1:6279917/1‑150 (MQ=255)
                                      tcgttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAg                                                                                                                <  2:7170501/150‑1 (MQ=255)
                                       cgttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAgg                                                                                                               <  2:453410/150‑1 (MQ=255)
                                        gttGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAggg                                                                                                              >  2:6226398/1‑150 (MQ=255)
                                           gTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCAt                                                                                                           >  2:8766412/1‑150 (MQ=255)
                                            tCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATc                                                                                                          <  2:3753065/150‑1 (MQ=255)
                                             ccGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCa                                                                                                         <  1:4290740/150‑1 (MQ=255)
                                              cggcggCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCAc                                                                                                        >  2:3544683/1‑150 (MQ=255)
                                                gcggcATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACca                                                                                                      >  1:12835/1‑150 (MQ=255)
                                                  ggcATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACcatc                                                                                                    <  1:2985980/150‑1 (MQ=255)
                                                     aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                                 <  2:3272396/150‑1 (MQ=255)
                                                     aTGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCatcatc                                                                                                 >  2:5724494/1‑150 (MQ=255)
                                                          gcggcgCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTcac                                                                                            >  1:1547743/1‑150 (MQ=255)
                                                             gcgcCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTcaccgc                                                                                         >  2:8072213/1‑148 (MQ=255)
                                                                ccTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAGCCAGAAGGGCATCACCATCATCCTCAccaccca                                                                                      >  2:1547240/1‑150 (MQ=255)
                                                                   gatgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTa                                                                                   >  1:8725561/1‑150 (MQ=255)
                                                                   gatgatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTa                                                                                   >  2:2656177/1‑150 (MQ=255)
                                                                      gatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCt                                                                                <  2:3702381/150‑1 (MQ=255)
                                                                      gatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCt                                                                                <  1:4553053/150‑1 (MQ=255)
                                                                      gatTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCt                                                                                >  2:1902139/1‑150 (MQ=255)
                                                                       atTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                               >  1:2362089/1‑150 (MQ=255)
                                                                       atTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTg                                                                               <  1:6140171/150‑1 (MQ=255)
                                                                        ttGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTgg                                                                              <  2:8623226/150‑1 (MQ=255)
                                                                         tGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGa                                                                             <  1:3303771/150‑1 (MQ=255)
                                                                          gCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGaa                                                                            <  2:3340704/150‑1 (MQ=255)
                                                                          gCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGaa                                                                            <  1:1174037/150‑1 (MQ=255)
                                                                            ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                          >  2:4992606/1‑150 (MQ=255)
                                                                            ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                          <  1:5635647/150‑1 (MQ=255)
                                                                            ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                          <  1:7318598/150‑1 (MQ=255)
                                                                            ccGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAaga                                                                          >  1:1044101/1‑150 (MQ=255)
                                                                                gcATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCt                                                                      <  2:1890917/150‑1 (MQ=255)
                                                                                  aTTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGa                                                                    >  1:646699/1‑150 (MQ=255)
                                                                                        ccACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCt                                                              >  1:8597974/1‑150 (MQ=255)
                                                                                        ccACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCt                                                              <  1:4365080/150‑1 (MQ=255)
                                                                                         cACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtg                                                             >  1:8287811/1‑150 (MQ=255)
                                                                                         cACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtg                                                             >  2:6690592/1‑150 (MQ=255)
                                                                                          aCGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCtgt                                                            >  2:2748785/1‑150 (MQ=255)
                                                                                            gAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGc                                                          <  2:3616444/150‑1 (MQ=255)
                                                                                              aCCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCg                                                        >  1:7459660/1‑150 (MQ=255)
                                                                                                     ctgctgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATc                                                 <  1:8746050/150‑1 (MQ=255)
                                                                                                      tgctgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCg                                                >  2:6893396/1‑150 (MQ=255)
                                                                                                       gctgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCt                                                                                >  1:3674714/1‑117 (MQ=255)
                                                                                                         tgATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGca                                             >  2:6837189/1‑150 (MQ=255)
                                                                                                             ccTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGcatcat                                         <  1:4647311/150‑1 (MQ=255)
                                                                                                             ccTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGcatcat                                         <  2:5522850/150‑1 (MQ=255)
                                                                                                              cTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCatcatc                                        <  2:2124819/150‑1 (MQ=255)
                                                                                                               tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                                       >  2:869548/1‑150 (MQ=255)
                                                                                                               tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                                       >  1:8184073/1‑150 (MQ=255)
                                                                                                               tGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCg                                       >  2:4438864/1‑150 (MQ=255)
                                                                                                                       tccACCCCGGCGGTGGGTATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGc                               >  1:4897851/2‑150 (MQ=255)
                                                                                                                        ccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCa                              >  2:4286400/1‑150 (MQ=255)
                                                                                                                        ccACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCa                              >  2:6018095/1‑150 (MQ=255)
                                                                                                                         cACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCAc                             >  1:2478421/1‑150 (MQ=255)
                                                                                                                         cACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCAc                             >  1:6558246/1‑150 (MQ=255)
                                                                                                                         cACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCAc                             <  1:1958102/150‑1 (MQ=255)
                                                                                                                         cACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCAc                             <  1:4828481/150‑1 (MQ=255)
                                                                                                                                  ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTCTCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga                    <  2:4002130/150‑1 (MQ=255)
                                                                                                                                  ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCgg                                               <  2:290335/123‑1 (MQ=255)
                                                                                                                                  ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga                    <  2:4800054/150‑1 (MQ=255)
                                                                                                                                  ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga                    <  1:7278100/150‑1 (MQ=255)
                                                                                                                                  ggTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCga                    <  2:8229173/150‑1 (MQ=255)
                                                                                                                                   gTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCgag                   <  2:6279917/150‑1 (MQ=255)
                                                                                                                                        tatCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAacac              >  1:8346385/1‑150 (MQ=255)
                                                                                                                                          tCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCa            >  1:287985/1‑150 (MQ=255)
                                                                                                                                             aGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCa         >  1:5079419/1‑150 (MQ=255)
                                                                                                                                              gCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCAt        >  2:7809027/1‑150 (MQ=255)
                                                                                                                                                tGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATgc      <  2:4876748/150‑1 (MQ=255)
                                                                                                                                                tGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATgc      <  2:7545456/150‑1 (MQ=255)
                                                                                                                                                tGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATgc      <  2:410274/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTGGGCCTGTGGGACAAACGTGATGTGCAGTCGCGTTCGTTGTCCGGCGGCATGAAGCGGCGCCTGATGATTGCCCGCGCATTGATCCACGAACCGCGCCTGCTGATCCTGGACGAACCCACCGCCGGGGTGGATATCGAGCTGCGCCGCTCGATGTGGAGCTTCCTCACTGAGCTGAACCAGAAGGGCATCACCATCATCCTCACCACCCACTATCTGGAAGAGGCTGAGCAGCTGTGCCGTAACATCGGCATCATCGACCATGGCACCATCGTCGAGAACACCAGCATGCGCCA  >  NC_002947/2143257‑2143554

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: