Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,874 G→C 14.8% A216A (GCG→GCC yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,8740GC14.8% 1051.4 / 145.2 359A216A (GCG→GCCyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (141/164);  new base C (27/26);  total (169/190)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 361 total aligned reads displayed.

CCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGG  >  NC_002947/1486729‑1487023
                                                                                                                                                 |                                                                                                                                                     
ccAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAAcc                                                                                                                                                   >  2:6781678/1‑150 (MQ=255)
  aaGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCgg                                                                                                                                                 <  2:3329238/150‑1 (MQ=255)
   aGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCggg                                                                                                                                                <  1:688173/150‑1 (MQ=255)
    gATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGt                                                                                                                                               >  2:6603751/1‑150 (MQ=255)
        aTGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATc                                                                                                                                           <  2:6316842/150‑1 (MQ=255)
           cTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTc                                                                                                                                        <  1:7804623/150‑1 (MQ=9)
            tgtgCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCa                                                                                                                                       >  1:1464122/1‑150 (MQ=255)
            tGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCa                                                                                                                                       >  1:7124385/1‑150 (MQ=12)
            tGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGACCCGGGTTATCTTCa                                                                                                                                       >  2:5023956/1‑150 (MQ=255)
            tGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCa                                                                                                                                       <  2:3030469/150‑1 (MQ=255)
             gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                                                                                                                                      <  1:2801957/150‑1 (MQ=12)
             gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                                                                                                                                      >  1:6649297/1‑150 (MQ=12)
             gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                                                                                                                                      <  1:1259590/150‑1 (MQ=12)
             gTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCAt                                                                                                                                      <  2:1322233/150‑1 (MQ=12)
              ttCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCAt                                                                                                                                      >  2:496942/1‑149 (MQ=255)
              ttCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATg                                                                                                                                     <  1:6700837/150‑1 (MQ=255)
                 gacgaGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGAc                                                                                                                                  >  2:3114338/1‑150 (MQ=255)
                     aGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGatcagg                                                                                                                              <  1:5571488/150‑6 (MQ=255)
                      gCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAggg                                                                                                                             >  2:6011744/1‑150 (MQ=255)
                      gCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAggg                                                                                                                             >  1:2938278/1‑150 (MQ=255)
                       ccGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGc                                                                                                                            >  2:4017292/1‑150 (MQ=255)
                       ccGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGc                                                                                                                            <  1:398615/150‑1 (MQ=255)
                        cGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGcc                                                                                                                           <  2:3236883/150‑2 (MQ=255)
                        cGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGAACTTCATGGATCAGGGcc                                                                                                                           <  1:4564118/150‑2 (MQ=255)
                         gACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGcca                                                                                                                          <  1:5764292/150‑3 (MQ=255)
                            cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                       >  2:8285526/1‑150 (MQ=255)
                            cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGGTATCTTCATGGACAAGGGCGAGAt                                                                                                                       >  1:7341508/1‑150 (MQ=255)
                             tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                      <  1:8700521/150‑1 (MQ=255)
                             tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                      >  1:2067581/1‑150 (MQ=255)
                             tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                      <  2:7091924/150‑1 (MQ=255)
                              cGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGAtcg                                                                                                                     <  1:6950756/150‑3 (MQ=255)
                               ggCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgt                                                                                                                    >  2:2526735/1‑150 (MQ=255)
                                  gctcgatccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGa                                                                                                                 >  2:1112654/8‑150 (MQ=255)
                                  gctcgatccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTCGa                                                                                                                 >  1:1007624/8‑150 (MQ=255)
                                   ctcgACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCATCACAATCATCTGCGTAACGCATGAAATGGGTTTTGCGCGCAAGATCGCCGATCGCGTCATCTTTATGGATCAGGGCAAGATCGTCGAg                                                                                                                <  1:8066303/147‑1 (MQ=255)
                                    tcgACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGACGGCATCACAATCATCTGCGTAACGCATGAAATGGGTTTTGCGCGCAAGATCGCCGATCGCGTCATCTTTATGGATCAGGGCAAGATCGTCGAgg                                                                                                               <  2:8117738/148‑2 (MQ=255)
                                    tGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGc                                                                                                               >  1:6242539/1‑150 (MQ=255)
                                     ggACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCa                                                                                                              <  2:1003370/150‑1 (MQ=255)
                                     ggACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCa                                                                                                              >  1:548508/1‑150 (MQ=255)
                                     ggACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCa                                                                                                              <  1:3226059/150‑1 (MQ=255)
                                        tccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAgaa                                                                                                           >  1:7496934/2‑148 (MQ=255)
                                         ccGGAAGTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                          >  1:7297190/1‑150 (MQ=255)
                                         ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                          >  2:1914822/1‑150 (MQ=255)
                                           ggAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCcgc                                                                                                        >  2:5973795/1‑150 (MQ=255)
                                           ggAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCcgc                                                                                                        >  1:5330545/1‑150 (MQ=255)
                                             aaaTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGTCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCgcgc                                                                                                      <  1:898658/150‑1 (MQ=255)
                                             aaaTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCgcgc                                                                                                      >  2:5288195/1‑150 (MQ=255)
                                               aTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCg                                                                                                    >  1:1637500/1‑150 (MQ=255)
                                                 ggTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGa                                                                                                  <  2:164898/150‑1 (MQ=255)
                                                  gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                 <  1:3343451/150‑1 (MQ=255)
                                                     aaGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGAc                                                                                              <  2:5658208/150‑1 (MQ=255)
                                                     aaGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGAc                                                                                              >  2:4517400/1‑150 (MQ=255)
                                                          aGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActtctt                                                                                         <  2:2689559/150‑1 (MQ=255)
                                                           gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                        <  1:7709062/150‑1 (MQ=255)
                                                                 gATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCgg                                                                                  >  1:6545367/1‑150 (MQ=255)
                                                                     ccATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCcgc                                                                              <  1:881419/150‑1 (MQ=255)
                                                                      cATGGTGGGCCTGGCCGATGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCcgcg                                                                             >  1:4959914/1‑150 (MQ=255)
                                                                      cATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCcgcg                                                                             >  1:1214025/1‑150 (MQ=255)
                                                                      cATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCcgcg                                                                             >  2:4589725/1‑150 (MQ=255)
                                                                      cATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCcgcg                                                                             >  1:2280442/1‑150 (MQ=255)
                                                                       aTGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCgcgc                                                                            >  2:9070188/1‑150 (MQ=255)
                                                                            gggCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGa                                                                       <  1:419730/150‑1 (MQ=255)
                                                                             ggCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGAc                                                                      <  2:4783465/150‑1 (MQ=255)
                                                                              gCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGTCTTCTTCGACCGGCCGCGCAGTGAcc                                                                     >  1:1910922/1‑150 (MQ=255)
                                                                               ccTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCg                                                                    <  1:7042985/150‑1 (MQ=255)
                                                                                     cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCaa                                                              >  2:4718389/1‑150 (MQ=255)
                                                                                     cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCaa                                                              <  2:6184369/150‑1 (MQ=255)
                                                                                     cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCaa                                                              >  2:650138/1‑150 (MQ=255)
                                                                                     cGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCaa                                                              <  1:8770042/150‑1 (MQ=255)
                                                                                        ggATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgtt                                                           <  2:3783357/150‑1 (MQ=255)
                                                                                           tGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGt                                                                         >  2:7370335/1‑133 (MQ=255)
                                                                                           tGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAgttgtt                                                        <  2:2346913/150‑1 (MQ=255)
                                                                                              cATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCt                                                     >  2:2805076/1‑150 (MQ=255)
                                                                                              cATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCt                                                     >  1:8385700/1‑150 (MQ=255)
                                                                                               aTGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTc                                                    >  1:7823962/1‑150 (MQ=255)
                                                                                               aTGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTc                                                    >  1:4979335/1‑150 (MQ=255)
                                                                                                tGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCa                                                   <  2:8471207/150‑1 (MQ=255)
                                                                                                 gACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAg                                                  >  1:3697155/1‑150 (MQ=255)
                                                                                                  aCCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGTGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGc                                                 >  2:2248148/1‑150 (MQ=255)
                                                                                                  aCCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGc                                                 <  2:8592204/150‑1 (MQ=255)
                                                                                                   ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                        <  2:8701857/110‑1 (MQ=255)
                                                                                                   ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGcc                                                >  2:8333444/1‑150 (MQ=255)
                                                                                                   ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGcc                                                <  2:4917509/150‑1 (MQ=255)
                                                                                                    cATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCa                                               <  1:5700622/150‑1 (MQ=255)
                                                                                                    cATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCa                                               <  1:632221/150‑1 (MQ=255)
                                                                                                     aTGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGCACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAg                                              >  1:4441280/1‑150 (MQ=255)
                                                                                                          ctGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCt                                         <  2:2458112/150‑1 (MQ=255)
                                                                                                           tGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTg                                        >  1:2481919/1‑150 (MQ=255)
                                                                                                             cGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCa                                      <  2:2090989/150‑1 (MQ=255)
                                                                                                              gTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAt                                     >  2:7583186/1‑150 (MQ=255)
                                                                                                                   ccACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAtt                                    <  1:4676076/146‑1 (MQ=255)
                                                                                                                   ccACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAtt                                    >  2:4676076/1‑146 (MQ=255)
                                                                                                                          aTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                         <  1:6781678/150‑1 (MQ=255)
                                                                                                                                tttGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTgc                   <  2:1543660/150‑1 (MQ=255)
                                                                                                                                    cccGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCgg               >  2:4104719/1‑150 (MQ=255)
                                                                                                                                       gTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAG            <  1:5906133‑M1/150‑1 (MQ=255)
                                                                                                                                          ccGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCC         >  1:7568129‑M1/1‑150 (MQ=255)
                                                                                                                                              ggCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCCGCGCAGCCCaaac     <  1:3151367‑M1/150‑5 (MQ=255)
                                                                                                                                               gCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGC    <  2:5411122‑M1/150‑1 (MQ=255)
                                                                                                                                                cGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCG   <  2:7823955‑M1/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGG  >  NC_002947/1486729‑1487023

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: