Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,883 C→T 9.6% P921L (CCA→CTA) ‡ PP_3413 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,863,8830CT 9.6% 269.8 / 17.3 104P921L (CCA→CTA) ‡PP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (49/45);  new base T (7/3);  total (56/48)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

Only 100 of 113 total aligned reads displayed.

TGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCGT  >  NC_002947/3863736‑3864032
                                                                                                                                                     |                                                                                                                                                     
tgACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCcag                                                                                                                                                     >  2:1212090/1‑150 (MQ=255)
   cacTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCcagcag                                                                                                                                                  <  2:1564344/150‑1 (MQ=255)
    acTGGCGGAGAATGGCGAGATCGCCCTTGCACGGTGGCAGGCAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCcagctcg                                                                                                                                                 >  1:197461/1‑147 (MQ=255)
      tGGCGGAGACTGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGAGTTGGTGATCACCGATTGCAGCATGCCGGTCATTGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGcc                                                                                                                                               >  1:1185021/1‑150 (MQ=255)
        gCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGt                                                                                                                                             >  1:1128029/1‑150 (MQ=255)
            agaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCTGGCAGGCCgcgctc                                                                                                                                         <  2:1436130/150‑6 (MQ=255)
             gaATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCt                                                                                                                                        >  1:1085352/1‑150 (MQ=255)
                tGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCCGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCCCCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTgg                                                                                                                                     <  2:1741253/150‑1 (MQ=255)
                    gagaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                    <  1:584414/147‑1 (MQ=255)
                    gagaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                    >  2:584414/1‑147 (MQ=255)
                           cccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCACG‑GCCTGCTGGCAGGCCgcgc                                                                                                                                           >  2:1824695/1‑130 (MQ=255)
                            ccTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCaa                                                                                                                         >  2:160060/1‑150 (MQ=255)
                              ttGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGCTTGCAGCATGCCGGTCATGGACGGGCATGAACTCCCCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATg                                                                                                                       <  2:2047643/150‑1 (MQ=255)
                              ttGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGTAGGCTGTATTCTTGGCGTGACCGCCAATg                                                                                                                       <  1:571158/150‑1 (MQ=255)
                                  aCGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGa                                                                                                                                >  2:1517229/1‑137 (MQ=255)
                                    gCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAgg                                                                                                                 >  1:124794/1‑150 (MQ=255)
                                      tGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGAGTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCACTTGCAAGGCCCCTGTAGGAGCG‑G‑GCTTGCCAGCAGGCCGTATTCTTGGCATGACCGCCACTGCCTAGGct                                                                                                               >  1:1488251/1‑149 (MQ=255)
                                        gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCga                                                                                                             <  2:124794/150‑1 (MQ=255)
                                         cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCgag                                                                                                            <  1:1934821/150‑1 (MQ=255)
                                         cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCgag                                                                                                            <  1:1644473/150‑1 (MQ=255)
                                         cAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCgag                                                                                                            <  1:160060/150‑1 (MQ=255)
                                          aGGAAGCGCGCTTCGACCTTGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGTCAAGCGCG‑GCATGCCAGCAGGCCGCATTCTTGGCGTGACCGCCAATGCCCAGGCCgagg                                                                                                           <  1:54439/150‑1 (MQ=255)
                                             aaGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGc                                                                                                        >  2:2276597/1‑150 (MQ=255)
                                               gcgcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGTGTGACCGCCAATGCCCAGGCCGAGGAGCGt                                                                                                      >  2:183780/1‑150 (MQ=255)
                                                   gcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCCCCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTc                                                                                                  <  1:107837/150‑1 (MQ=255)
                                                   gcTTCGACCTGGTGATCACCGATTGAGTCCAGTCGGTCATGGACGGGCATGAACTCACCCGGCGCATTGGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTc                                                                                                  >  1:1269101/1‑150 (MQ=255)
                                                        gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                                                                                             <  1:2276597/150‑1 (MQ=255)
                                                         atcttgTGAGCAGCGTTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGcc                                                                                            >  2:1610335/6‑150 (MQ=255)
                                                          ccTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCt                                                                                           <  1:1255364/150‑1 (MQ=255)
                                                           cTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTg                                                                                          >  2:1880855/1‑150 (MQ=255)
                                                            tGGTGATCACCGATTGCAGCATGCCGGTCCTGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTgg                                                                                         >  2:1849489/1‑150 (MQ=255)
                                                            tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                    >  1:91640/1‑107 (MQ=255)
                                                            tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGc                                                                                                                                    <  2:91640/107‑1 (MQ=255)
                                                            tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTgg                                                                                         >  2:848794/1‑150 (MQ=255)
                                                            tGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTgg                                                                                         >  1:757302/1‑150 (MQ=255)
                                                             ggTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGc                                                                                        <  2:1572697/150‑1 (MQ=255)
                                                                gATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAg                                                                                     >  1:1201724/1‑150 (MQ=255)
                                                                 aTCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGc                                                                                    <  1:1501825/150‑1 (MQ=255)
                                                                 aTCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGc                                                                                    <  2:1128029/150‑1 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTAtt                                                                                                                                          <  1:450242/93‑1 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTAtt                                                                                                                                          >  2:450242/1‑93 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATg                                                                                                                       <  1:1133161/112‑1 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGc                                                                                 >  1:2159359/1‑150 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCACG‑GCCTGCTGGCAGGCCgcgc                                                                                                                                           >  2:524885/1‑89 (MQ=255)
                                                                    aCCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCACG‑GCCTGCTGGCAGGCCgcgc                                                                                                                                           <  1:524885/92‑4 (MQ=255)
                                                                      cGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCAt                                                                               >  2:792016/1‑150 (MQ=255)
                                                                        aTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATgg                                                                             <  1:1212090/150‑1 (MQ=255)
                                                                        aTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATgg                                                                             <  1:1610335/150‑1 (MQ=255)
                                                                            cagcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGT‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGa                                                                         <  1:105278/150‑1 (MQ=255)
                                                                             agcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                      <  1:324739/88‑1 (MQ=255)
                                                                             agcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                      >  2:324739/1‑88 (MQ=255)
                                                                              gcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCg                                                                                   >  1:1881087/1‑138 (MQ=255)
                                                                              gcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAAt                                                                       <  1:848794/150‑1 (MQ=255)
                                                                                  gCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGc                                                                                                        <  2:1201859/113‑1 (MQ=255)
                                                                                     ggTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGCGCTC‑ATCTTGTCCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGtt                                                                <  2:1269101/150‑1 (MQ=255)
                                                                                     ggTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGtt                                                                <  1:536538/150‑1 (MQ=255)
                                                                                      gTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGGCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTAGTCGCCGCCTGGCCAGGGGCATGGACGAATGCCTGTTc                                                               >  1:812544/1‑150 (MQ=255)
                                                                                        cATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCaa                                                             <  2:2159359/150‑1 (MQ=255)
                                                                                         aTGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAg                                                            >  1:207918/1‑150 (MQ=255)
                                                                                          tGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGc                                                           >  2:1212443/1‑150 (MQ=255)
                                                                                           ggACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCGGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGcc                                                          >  2:474680/1‑150 (MQ=255)
                                                                                              ggggCATGAACTCACCCGGGGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCTTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAt                                                       <  2:566499/149‑1 (MQ=255)
                                                                                              cGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAt                                                       <  1:1634414/150‑1 (MQ=255)
                                                                                               gggCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAtt                                                      >  2:437382/1‑150 (MQ=255)
                                                                                                 gCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTgg                                                    >  1:1159254/1‑150 (MQ=255)
                                                                                                   aTGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCcagca                                                                                                                                                   <  2:2006265/53‑1 (MQ=255)
                                                                                                   aTGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCcagca                                                                                                                                                   >  1:2006265/1‑53 (MQ=255)
                                                                                                    tGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCt                                                 >  2:5683/1‑150 (MQ=255)
                                                                                                        cTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAAt                                                                       >  1:219543/1‑124 (MQ=255)
                                                                                                        cTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACGCGCTGCCTGGCCAGCGGCATGGACGAAt                                                                       <  2:219543/124‑1 (MQ=255)
                                                                                                           aCCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                      <  2:1717374/58‑1 (MQ=255)
                                                                                                           aCCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTg                                                                                                                                      >  1:1717374/1‑58 (MQ=255)
                                                                                                             ccGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGACGTCTTATAGGCGTGAGCACCAATGGCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCt                                        >  2:33808/1‑150 (MQ=255)
                                                                                                                  gcATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGa                                   <  2:1085352/150‑1 (MQ=255)
                                                                                                                  gcATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGa                                   >  2:451962/1‑150 (MQ=255)
                                                                                                                    aTTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGAcc                                 >  2:1973829/1‑150 (MQ=255)
                                                                                                                      tCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCa                               <  2:1159254/150‑1 (MQ=255)
                                                                                                                           aGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTgc                                               <  1:2209657/130‑1 (MQ=255)
                                                                                                                           aGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTgc                                               >  2:2209657/1‑130 (MQ=255)
                                                                                                                                 gAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAg                     <  1:1278590/150‑1 (MQ=255)
                                                                                                                                   gggCGACCGCG‑GCCGGCCAGTTGTCCGTATTCTTGGCGTGACCGCCAATGCCCAGTCCGAGTAGCGTACTCGCTGCTTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGCCCCACCTGCCCCAGGt                   <  2:45020/150‑1 (MQ=255)
                                                                                                                                     gCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGcccc                       >  2:2035919/1‑144 (MQ=255)
                                                                                                                                      agaatacg‑GCCTGCTGGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTgcg                >  1:2001556/7‑150 (MQ=255)
                                                                                                                                        aGCGCG‑GCCGTCCAACAGGCTGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATCGACGTATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGcc              <  2:1185021/150‑1 (MQ=255)
                                                                                                                                         gcgcG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGccc             >  1:2400515/1‑150 (MQ=255)
                                                                                                                                         atacg‑GCCTGCTGGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGccc             >  1:2032258/4‑150 (MQ=255)
                                                                                                                                          tacg‑GCCTGCTGGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCt            >  1:1806815/3‑150 (MQ=255)
                                                                                                                                          tacg‑GCCTGCTGGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCt            >  2:260873/3‑150 (MQ=255)
                                                                                                                                          tacg‑GCCTGCTGGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCt            >  1:732948/3‑150 (MQ=255)
                                                                                                                                          cgcg‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCt            >  2:1566767/1‑150 (MQ=255)
                                                                                                                                           gcg‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGa                                   <  2:1670112/126‑1 (MQ=255)
                                                                                                                                           gcg‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGa                                   >  1:1670112/1‑126 (MQ=255)
                                                                                                                                            cg‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTca          <  1:1698089/150‑1 (MQ=255)
                                                                                                                                                 cTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTcagcag      >  1:59548/1‑150 (MQ=255)
                                                                                                                                                  tGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGc                                                                     <  1:5077/86‑1 (MQ=255)
                                                                                                                                                  tGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCAc                              <  1:1208525/125‑1 (MQ=255)
                                                                                                                                                  tGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCAc                              >  2:1208525/1‑125 (MQ=255)
                                                                                                                                                    ccAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGTTCGGATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCg   >  1:2098795/1‑150 (MQ=255)
                                                                                                                                                    ccAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCg   >  2:182181/1‑150 (MQ=255)
                                                                                                                                                     cagcagGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCGt  <  1:1880855/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGC‑AGCCTGGAGGGCGAGCGCG‑GCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGATTGGTCTGCGCACGCTCAAGACCCACCTGCCCCAGGTGCGCCCTCAGCAGCCGT  >  NC_002947/3863736‑3864032

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