Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,597,893 T→C 7.6% 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 7.6% 269.0 / 8.8 92L465P (CTG→CCG) PP_4923outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base T (47/38);  new base C (5/2);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATATCA  >  NC_002947/5597746‑5598041
                                                                                                                                                    |                                                                                                                                                     
aCTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTg                                                                                                                                                      <  6:223713/149‑1 (MQ=255)
  tACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggc                                                                                                                                                    >  3:178084/1‑149 (MQ=255)
    caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                                  >  2:86621/1‑149 (MQ=255)
    caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                                  >  2:162728/1‑149 (MQ=255)
    caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcg                                                                                                                                                   <  4:178084/148‑1 (MQ=255)
             ccGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCgg                                                                                                                                         <  2:229133/149‑1 (MQ=255)
               gTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCGGCCGGACCTGGCGGGAGCAGCGGcc                                                                                                                                       <  8:227771/149‑1 (MQ=255)
               gTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGcc                                                                                                                                       <  3:351508/149‑1 (MQ=255)
                  aCGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCaa                                                                                                                                     <  3:56773/148‑1 (MQ=255)
                      cTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACt                                                                                                                                >  3:292456/1‑149 (MQ=255)
                         cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                              >  1:374218/1‑148 (MQ=255)
                         cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                              <  5:216092/148‑1 (MQ=255)
                         aaTTATCGACAACCTGAGCCTGAAGCAGGGAGCCGTAACCCTGAGCCCGCTAGACCTTGAGGACGTGAAGCGCTACCTGAAACTTGACTACAACCATGACAATGACTTACTGCGGCATGACCTGGCGGAAGCAGCTGACAATAACTTc                                                                                                                              <  5:171238/147‑1 (MQ=255)
                          aTTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGTAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                            <  6:242048/149‑1 (MQ=255)
                          aTTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                            <  4:63744/149‑1 (MQ=255)
                           ttCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcagg                                                                                                                                              <  4:110755/130‑2 (MQ=255)
                           ttCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcagg                                                                                                                                              >  3:110755/1‑129 (MQ=255)
                               cGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAgcg                                                                                                                         >  3:367249/1‑147 (MQ=255)
                                gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAgcgc                                                                                                                        <  2:256480/147‑1 (MQ=255)
                                   aaCCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGcc                                                                                                                   >  1:95170/1‑149 (MQ=255)
                                   aaCCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCCGACCTGGCGGCAGCCGCGGCCAAGAACTTCGAGCGCAGGcc                                                                                                                   >  5:319934/1‑149 (MQ=255)
                                     ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                                 >  7:296219/1‑149 (MQ=255)
                                     ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                                 <  3:243404/149‑1 (MQ=255)
                                      cTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                                 >  3:247744/1‑148 (MQ=255)
                                       tGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                                 >  7:303097/1‑147 (MQ=255)
                                       tGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                                 >  8:210922/1‑147 (MQ=255)
                                         aGCCTGAAGCATGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTga                                                                                                               <  5:293538/147‑1 (MQ=255)
                                         aGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTga                                                                                                               <  1:86621/148‑1 (MQ=255)
                                           ccTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTgaga                                                                                                             >  2:280815/1‑147 (MQ=255)
                                            cTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAg                                                                                                          >  1:9113/1‑149 (MQ=255)
                                               aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                         <  6:66910/147‑1 (MQ=255)
                                                     gcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTggtc                                                                                                 >  8:150744/1‑147 (MQ=255)
                                                                ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGc                                                                                      >  1:23801/1‑149 (MQ=255)
                                                                ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGc                                                                                      >  6:253407/1‑149 (MQ=255)
                                                                  tGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCtt                                                                                    >  2:184608/1‑149 (MQ=255)
                                                                       ccGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGCCCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCg                                                                               >  3:229095/1‑149 (MQ=255)
                                                                        cGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                         <  5:381637/122‑1 (MQ=255)
                                                                        cGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                         >  6:381637/1‑122 (MQ=255)
                                                                        cGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCAGc                                                                              >  5:308755/1‑149 (MQ=255)
                                                                         gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAg                                                                                                                           >  3:288640/1‑103 (MQ=255)
                                                                         gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAg                                                                                                                           <  4:288640/103‑1 (MQ=255)
                                                                         gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgggg                                                                                                       >  5:75417/1‑123 (MQ=255)
                                                                         gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgggg                                                                                                       <  6:75417/123‑1 (MQ=255)
                                                                             gACCTGCAGGACCTGACGAGCTACCTGAAACCTGCCTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCca                                                                         <  8:296219/149‑1 (MQ=255)
                                                                             gACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCca                                                                         >  5:339891/1‑149 (MQ=255)
                                                                              aCCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCca                                                                         >  1:231458/1‑148 (MQ=255)
                                                                              aCCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCAAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCca                                                                         >  1:336720/1‑148 (MQ=255)
                                                                                cTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCacac                                                                     >  2:343747/1‑149 (MQ=255)
                                                                                   cAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGt                                                                  <  4:367249/149‑1 (MQ=255)
                                                                                    aGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTg                                                                 <  1:184608/149‑1 (MQ=255)
                                                                                     ggACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGc                                                                <  1:162728/149‑1 (MQ=255)
                                                                                         cTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAgg                                                            >  6:203669/1‑149 (MQ=255)
                                                                                         cTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAgg                                                            <  2:56481/149‑1 (MQ=255)
                                                                                          tGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGc                                                           >  1:331132/1‑149 (MQ=255)
                                                                                           gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGcc                                                          >  3:83492/1‑149 (MQ=255)
                                                                                            aaGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGc                                                           <  1:343747/147‑1 (MQ=255)
                                                                                            aaGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCt                                                         >  3:330235/1‑149 (MQ=255)
                                                                                              gcgcTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGa                                                       <  4:292456/149‑1 (MQ=255)
                                                                                                  tACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAgg                                                    >  3:253332/1‑148 (MQ=255)
                                                                                                    ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGc                                                   <  2:23801/147‑1 (MQ=255)
                                                                                                    ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCat                                                 <  2:374218/149‑1 (MQ=255)
                                                                                                    ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCat                                                 >  3:58472/1‑149 (MQ=255)
                                                                                                     cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCat                                                 >  2:367074/1‑148 (MQ=255)
                                                                                                          aCCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCa                                           >  3:292163/1‑149 (MQ=255)
                                                                                                           ccTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                          >  1:288487/1‑149 (MQ=255)
                                                                                                           ccTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCAAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                          <  2:336720/149‑1 (MQ=255)
                                                                                                            cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                         >  8:342527/1‑149 (MQ=255)
                                                                                                            cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTACCAGCACGGTGCCAGGCCTGAAGGCATGATGCAc                                          <  1:33668/149‑1 (MQ=255)
                                                                                                              gACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtg                                       >  8:158664/1‑149 (MQ=255)
                                                                                                                     aCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGc                                >  5:251311/1‑149 (MQ=255)
                                                                                                                     aCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGc                                >  7:201393/1‑149 (MQ=255)
                                                                                                                      cccGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCGTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCt                               >  8:39300/1‑149 (MQ=255)
                                                                                                                       ccGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                              >  7:267090/1‑149 (MQ=255)
                                                                                                                        cGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTgg                             >  2:364655/1‑149 (MQ=255)
                                                                                                                        cGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTgg                             <  7:210922/149‑1 (MQ=255)
                                                                                                                              aGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGCGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCt                       <  2:9113/149‑1 (MQ=255)
                                                                                                                               ggACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTg                      >  6:205532/1‑149 (MQ=255)
                                                                                                                               ggACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTg                      <  4:292163/149‑1 (MQ=255)
                                                                                                                                gACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTgg                     <  6:339891/149‑1 (MQ=255)
                                                                                                                                gACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTgg                     <  8:201393/149‑1 (MQ=255)
                                                                                                                                    tgcTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTg                      <  1:109182/142‑1 (MQ=255)
                                                                                                                                    tgcTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTg                      >  2:109182/3‑144 (MQ=255)
                                                                                                                                      cTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGACTGAAGGCATGATGCACCTGTGAGCGCTGGCTTGCTGGCGATAg               >  4:15564/1‑149 (MQ=255)
                                                                                                                                      cTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg               >  5:347862/1‑149 (MQ=255)
                                                                                                                                      cTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg               >  6:152283/1‑149 (MQ=255)
                                                                                                                                       tGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAgg              >  8:303081/1‑149 (MQ=255)
                                                                                                                                         ccgccgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAgg              <  4:58472/147‑1 (MQ=255)
                                                                                                                                         ccgccgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGcc            >  5:344684/1‑149 (MQ=255)
                                                                                                                                         ccgccgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGc             <  7:158664/148‑1 (MQ=255)
                                                                                                                                              cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTc         >  2:327073/1‑148 (MQ=255)
                                                                                                                                              cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCa        <  5:203669/149‑1 (MQ=255)
                                                                                                                                               ggACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCaa       <  8:303097/149‑1 (MQ=255)
                                                                                                                                                    tGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATATCa  >  5:364914/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTACAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAATATCA  >  NC_002947/5597746‑5598041

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