Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,597,893 T→C 5.7% L465P (CTG→CCG)  PP_4923 → outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,597,8930TC 5.7% 266.5 / 3.9 88L465P (CTG→CCG) PP_4923outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base T (44/39);  new base C (3/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATAT  >  NC_002947/5597745‑5598039
                                                                                                                                                         |                                                                                                                                                  
gACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCt                                                                                                                                                    >  5:337528/1‑149 (MQ=255)
gACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCt                                                                                                                                                    >  1:38682/1‑149 (MQ=255)
  cTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTgg                                                                                                                                                  <  1:210338/149‑1 (MQ=255)
  cTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTgg                                                                                                                                                  <  6:135871/149‑1 (MQ=255)
      aacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAACCCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTggcgg                                                                                                                                               >  6:303968/1‑148 (MQ=255)
        caaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGgcag                                                                                                                                            <  1:345803/149‑1 (MQ=255)
           tAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGcagcag                                                                                                                                         >  8:57060/1‑149 (MQ=255)
                 ttACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGcc                                                                                                                                    >  8:309921/1‑148 (MQ=255)
                 ttACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTCACT‑AC‑AACCCGGACAAG‑‑‑GACTTCATGCCGCCGGACCTGGCGGCAGCAGCGGcc                                                                                                                                    >  3:121124/1‑148 (MQ=255)
                   aCGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCaa                                                                                                                                  <  2:36679/148‑1 (MQ=255)
                     gACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGa                                                                                                                                >  7:117945/1‑148 (MQ=255)
                       cTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCACAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACt                                                                                                                             >  1:274705/1‑149 (MQ=255)
                        tACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCAGGACAAG‑‑‑GACTTCCTGCCGCAGGCCCTGGCGGCCGCAGCGGCCATGAAcgt                                                                                                                            >  8:185970/1‑147 (MQ=255)
                          cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           <  2:38682/148‑1 (MQ=255)
                          cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           <  7:57060/148‑1 (MQ=255)
                               tcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCAGGTCCGGCGgcaggaa                                                                                                                                         >  7:147062/1‑126 (MQ=255)
                               tcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCAGGTCCGGCGgcaggaa                                                                                                                                         <  8:147062/129‑4 (MQ=255)
                               tcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                            >  2:348406/1‑126 (MQ=255)
                               tcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                            <  1:348406/126‑1 (MQ=255)
                                 gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           <  1:46817/141‑1 (MQ=255)
                                 gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                           >  2:46817/1‑141 (MQ=255)
                                 gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAg                                                                                                                   <  3:265289/149‑1 (MQ=255)
                                      ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                              >  5:27890/1‑149 (MQ=255)
                                       gtgtgCCTGAAGCAGGCAGCCTGCACCCTGAGCCCGCAAGACCTGCAGGAACTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCTGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGGAGCGGCCAAGTACTTCGAGCGCAGGCGtt                                                                                                              <  4:121124/144‑1 (MQ=255)
                                         gAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                              >  3:327122/1‑146 (MQ=255)
                                          aGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTgaga                                                                                                          <  5:283781/149‑1 (MQ=255)
                                          aGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTga                                                                                                            <  4:262979/147‑1 (MQ=255)
                                              tGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                         >  2:24704/1‑130 (MQ=255)
                                              tGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                         <  1:24704/130‑1 (MQ=255)
                                                aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAggg                                                                                                     <  4:176026/148‑1 (MQ=255)
                                                aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                      <  6:337528/147‑1 (MQ=255)
                                                   caggcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTg                                                                                                 >  1:224302/1‑149 (MQ=255)
                                                      gcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGAc                                                                                              >  3:313389/1‑149 (MQ=255)
                                                      gcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGAc                                                                                              >  1:319281/1‑149 (MQ=255)
                                                          ccGGCCCCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTCGCGAACCCGGACACGCCCTACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGa                                                                                               <  1:121722/149‑1 (MQ=255)
                                                           cGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTc                                                                                         <  8:116982/149‑1 (MQ=255)
                                                            ggCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTcg                                                                                        >  7:290196/1‑149 (MQ=255)
                                                               aCCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCgg                                                                                     <  2:274705/149‑1 (MQ=255)
                                                                      gCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAACCCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCTGACCTGGCGGCAGAAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGcc                                                                              <  8:280361/149‑1 (MQ=255)
                                                                        ccGCAAGACCTGCAGGACCTGAAGCGCTACCAGAAACCTGACA‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACATCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCg                                                                            >  3:348077/1‑149 (MQ=255)
                                                                            aaGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTcc                                                                        >  3:75597/1‑149 (MQ=255)
                                                                            aaGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCAAGCCCGcgaa                                                                        <  6:41220/149‑4 (MQ=255)
                                                                                ccTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCcaca                                                                    >  2:372167/1‑149 (MQ=255)
                                                                                 cTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCacac                                                                   >  7:113357/1‑149 (MQ=255)
                                                                                  tGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGTGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACg                                                                  >  3:224306/1‑149 (MQ=255)
                                                                                    cAGGACCTGAAGCGCTACCTGAAACATGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGt                                                                >  6:365833/1‑149 (MQ=255)
                                                                                          cTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAgg                                                          <  2:319281/149‑1 (MQ=255)
                                                                                          cTGAAGCGCTACC‑GAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAgg                                                          >  5:180159/1‑148 (MQ=255)
                                                                                            gAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGcc                                                        >  7:71333/1‑149 (MQ=255)
                                                                                            gAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGcc                                                        >  3:225938/1‑149 (MQ=255)
                                                                                             aaGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGATTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCt                                                       >  4:195189/1‑149 (MQ=255)
                                                                                             aaGCGCTACC‑GAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCt                                                       <  6:180159/148‑1 (MQ=255)
                                                                                               gcgcTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCGTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGa                                                     <  5:303968/149‑1 (MQ=255)
                                                                                                  cTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTcgcg                                                                                      <  4:257152/113‑1 (MQ=255)
                                                                                                  cTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTcgcg                                                                                      >  3:257152/1‑113 (MQ=255)
                                                                                                   tACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAgg                                                  >  5:164845/1‑148 (MQ=255)
                                                                                                     ccTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGATTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGc                                                 <  3:195189/147‑1 (MQ=255)
                                                                                                      cTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGc                                                         <  7:73841/138‑1 (MQ=255)
                                                                                                      cTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGc                                                         >  8:73841/1‑138 (MQ=255)
                                                                                                      cTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCat                                               >  3:319405/1‑148 (MQ=255)
                                                                                                        gAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCatgac                                            >  2:82910/1‑148 (MQ=255)
                                                                                                           aCCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACATCGAGCGCAGCCCATGAGATAGCGGCTGGACTCTTCGCGGGCTTGCACGCTCCCACACCGAGCCAGGCCTCAAGGCATGATGCa                                         >  2:313842/1‑149 (MQ=255)
                                                                                                            ccTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                        <  1:372167/149‑1 (MQ=255)
                                                                                                                aCT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGACGCACCtg                                     <  1:82910/148‑1 (MQ=255)
                                                                                                                 cT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATCGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtgt                                    <  7:185970/148‑1 (MQ=255)
                                                                                                                  t‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtgt                                    >  2:12077/1‑147 (MQ=255)
                                                                                                                     c‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTga                                  <  6:27890/147‑1 (MQ=255)
                                                                                                                     c‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTga                                  <  4:327122/147‑1 (MQ=255)
                                                                                                                     c‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACATGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCACACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTgaga                                >  6:30698/1‑149 (MQ=255)
                                                                                                                       aaCCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGt                                                                >  6:307527/1‑116 (MQ=255)
                                                                                                                            ggACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGc                          <  4:225938/149‑1 (MQ=255)
                                                                                                                             gACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCt                         >  8:60299/1‑149 (MQ=255)
                                                                                                                                  g‑‑‑gACTTCCTGCCGCCGGACCTGTCGGCAGCAGCGGCCTAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTg                    <  4:348077/149‑1 (MQ=255)
                                                                                                                                      gACTTCCTGCCTCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTgg                   >  2:125794/1‑149 (MQ=255)
                                                                                                                                      gACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTgg                   <  1:12077/149‑1 (MQ=255)
                                                                                                                                      gACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTgg                   <  7:309921/149‑1 (MQ=255)
                                                                                                                                        cTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGc                  >  8:26714/1‑148 (MQ=255)
                                                                                                                                          tgcTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGa                >  8:176360/3‑148 (MQ=255)
                                                                                                                                          tCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGa                >  6:348093/1‑148 (MQ=255)
                                                                                                                                            cTGCCGCCGGACCTGGTGGCAGCAGCGGCCAAGAACTTCGAGCGCAGTCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg             <  4:224306/149‑1 (MQ=255)
                                                                                                                                             tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAgg            >  7:185332/1‑149 (MQ=255)
                                                                                                                                                 gccgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCtt        >  5:101267/1‑149 (MQ=255)
                                                                                                                                                  ccgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTc       <  4:187172/149‑1 (MQ=255)
                                                                                                                                                   cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCa      <  8:113357/149‑1 (MQ=255)
                                                                                                                                                     gACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCaa     >  6:26475/1‑148 (MQ=255)
                                                                                                                                                      aCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAata   <  8:117945/149‑1 (MQ=255)
                                                                                                                                                      aCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAata   <  5:365833/149‑1 (MQ=255)
                                                                                                                                                      aCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAata   >  6:123846/1‑149 (MQ=255)
                                                                                                                                                        cTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAatat  >  6:268927/1‑148 (MQ=255)
                                                                                                                                                        cTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAACGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGTCTTCAatat  >  5:3171/1‑148 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
GACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACT‑AC‑AACCCGGACAAG‑‑‑GACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATAT  >  NC_002947/5597745‑5598039

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: ‑