Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,495 C→G 6.5% P173P (CCG→CCC PP_1388 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,584,4950CG 6.5% 230.4 / 4.3 73P173P (CCG→CCCPP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base G (2/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGC  >  NC_002947/1584348‑1584640
                                                                                                                                                   |                                                                                                                                                 
ctggcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                                  <  4:407283/149‑1 (MQ=255)
     ccaccaAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAgg                                                                                                                                             >  1:65467/1‑149 (MQ=255)
      caccaAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGt                                                                                                                                            >  4:81254/1‑149 (MQ=255)
          aaGGCGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAt                                                                                                                                        >  7:65835/1‑149 (MQ=255)
                    ccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACAcc                                                                                                                              >  2:186264/1‑149 (MQ=255)
                          ggagGCTGCCAAGGCCGATAATCATCATCATGGCGCGATGGTAGTCCACGTGTGCCTCACGACTCTGCAGTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCTACCAGCCCCAGCGGCATGTTGATCCAGAACACGCAgcgc                                                                                                                        <  8:65835/147‑1 (MQ=255)
                                tGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCCGCGGCAGGTTGATCCAGAACACCCAGCGCCACGa                                                                                                                   >  2:40631/1‑148 (MQ=255)
                                  ccAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGAc                                                                                                                  <  6:13817/147‑1 (MQ=255)
                                   cAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAgg                                                                                                               >  8:323308/1‑149 (MQ=255)
                                    aaGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGt                                                                                                              <  3:81254/149‑1 (MQ=255)
                                      ggCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTAc                                                                                                            <  5:144024/149‑1 (MQ=255)
                                      ggCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTAc                                                                                                            <  8:307302/149‑1 (MQ=255)
                                         cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                                    <  8:154994/122‑1 (MQ=255)
                                         cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                                    >  7:154994/1‑122 (MQ=255)
                                         cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCg                                                                                                         <  6:33134/149‑1 (MQ=255)
                                           agaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGt                                                                                                       <  7:122003/149‑1 (MQ=255)
                                            gaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCATCCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTc                                                                                                      <  1:186264/149‑1 (MQ=255)
                                               tCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   >  3:38305/1‑149 (MQ=255)
                                                  aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAg                                                                                                >  7:48861/1‑149 (MQ=255)
                                                       aTGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                           <  4:38305/149‑1 (MQ=255)
                                                         ggCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                           >  1:23782/1‑147 (MQ=255)
                                                            gcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCa                                                                                      <  2:65467/149‑1 (MQ=255)
                                                             ccAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCAc                                                                                     <  8:48861/149‑1 (MQ=255)
                                                                      tCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGc                                                                            <  7:323308/149‑1 (MQ=255)
                                                                       ccACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGc                                                                            >  4:40083/1‑148 (MQ=255)
                                                                           cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCAc                                                                         <  5:115211/147‑1 (MQ=255)
                                                                           cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCAc                                                                         >  6:115211/1‑147 (MQ=255)
                                                                                  tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCcgc                                                                >  6:127587/1‑149 (MQ=255)
                                                                                  tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCcgc                                                                >  6:385464/1‑149 (MQ=255)
                                                                                    acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   <  4:38232/112‑1 (MQ=255)
                                                                                    acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                                   >  3:38232/1‑112 (MQ=255)
                                                                                    acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcgg                                                              >  2:303425/1‑149 (MQ=255)
                                                                                    acgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcgg                                                              >  4:367319/1‑149 (MQ=255)
                                                                                               cTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACAtgctg                                                   <  5:207803/149‑1 (MQ=255)
                                                                                                  gCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATgctgct                                                 >  3:195254/1‑148 (MQ=255)
                                                                                                  gCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATgctgct                                                 >  6:324272/1‑148 (MQ=255)
                                                                                                     tGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAg                                             <  5:127587/149‑1 (MQ=255)
                                                                                                       tCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCa                                                                                                                      >  3:182382/1‑74 (MQ=255)
                                                                                                       tCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCa                                                                                                                      <  4:182382/74‑1 (MQ=255)
                                                                                                          gcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGcc                                        >  2:354653/1‑149 (MQ=255)
                                                                                                             agcgCACGGCTAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTg                                     <  3:147602/149‑1 (MQ=255)
                                                                                                               cgcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGt                                   <  5:324272/149‑1 (MQ=255)
                                                                                                                  ccgtctcaTGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGACCCCCCAGCGACACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGTCGTACATGCTGATGAAGTAGCCCTg                                     <  4:77243/137‑1 (MQ=255)
                                                                                                                  aCGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGc                                <  5:385464/149‑1 (MQ=255)
                                                                                                                  aCGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGc                                >  3:144805/1‑149 (MQ=255)
                                                                                                                  aCGGCGAATGGCCCACAAAGCAACCAGCCACAGCGGCAGGTTGATCCATAACACCCAGCGCCACGACAGGTACTCTGTCAACCAGCAGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGccatc                                     >  3:77243/1‑141 (MQ=255)
                                                                                                                    ggCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCgg                              >  3:322895/1‑149 (MQ=255)
                                                                                                                      cGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGcc                            >  7:89087/1‑149 (MQ=255)
                                                                                                                       gAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcg                           >  6:245877/1‑149 (MQ=255)
                                                                                                                        aaTGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcgc                          <  5:245877/149‑1 (MQ=255)
                                                                                                                            gCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGtt                       <  4:195254/148‑1 (MQ=255)
                                                                                                                             cccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcg                     >  1:170479/1‑149 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACAtgc                                                     <  8:182547/116‑1 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACAtgc                                                     >  7:182547/1‑116 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcgc                          >  8:296959/1‑143 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCcgc                          <  7:296959/143‑1 (MQ=255)
                                                                                                                              ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCACCGCTGCCCGCGGCGTACAAGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgc                    >  2:20574/1‑149 (MQ=255)
                                                                                                                               cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcg                   >  6:244941/1‑149 (MQ=255)
                                                                                                                               cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcg                   >  4:243193/1‑149 (MQ=255)
                                                                                                                                  ggatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggc                                                             >  7:84127/5‑104 (MQ=38)
                                                                                                                                  ggatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggc                                                             <  8:84127/100‑1 (MQ=38)
                                                                                                                                  aaaGCAACCAGCCCCAGGGGGAGGTTGATCCAGAACACCGAGCGCCTCGACAGGTACTCGGTCAACGAGCCGGCCAGCAGTGGCCCGGCCACGCTGGCAGCGGCGTACATGGTGCTGAAGTAGCCCTGGTAGCGGCCGCTTTCGCTCgg                <  2:325144/149‑1 (MQ=255)
                                                                                                                                  aaaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACt                                                                                                           >  5:79865/1‑58 (MQ=255)
                                                                                                                                  aaaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACt                                                                                                           <  6:79865/58‑1 (MQ=255)
                                                                                                                                  aaaGCAACCAGCACCAGCGGCAGGTTGCTCCAGCACACCCAGCGCCACGACAGGTACACCGGCAACCAGACGCCCAGCACTGGCACGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgg                >  1:325144/1‑149 (MQ=255)
                                                                                                                                       aaCCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGcc                                        <  7:318227/120‑1 (MQ=255)
                                                                                                                                       aaCCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGcc                                        >  8:318227/1‑120 (MQ=255)
                                                                                                                                       aaCCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCa           >  2:276325/1‑149 (MQ=255)
                                                                                                                                        accaccCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAc          <  6:425715/149‑1 (MQ=255)
                                                                                                                                           aGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGaa       <  5:244941/149‑1 (MQ=255)
                                                                                                                                            gCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAg      <  4:322895/149‑1 (MQ=255)
                                                                                                                                               acAGCGGCAGGTTGATCAAGAACACCCAGAGCCACGACATGTACTCGGTCAACCAGCAGCCCAGCACTGGCCCGGCCACTCTGCAAGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCg   >  2:93072/2‑149 (MQ=255)
                                                                                                                                                cAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGc  <  2:23782/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGTCGC  >  NC_002947/1584348‑1584640

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