Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,641,185 G→C 6.4% G283A (GGC→GCC) ‡ PP_4951 → EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,641,1850GC 6.4% 201.1 / 3.8 62G283A (GGC→GCC) ‡PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (2/2);  total (28/34)
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.35e-01

GCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACGGTTCGGCTATC  >  NC_002947/5641039‑5641331
                                                                                                                                                  |                                                                                                                                                  
gctgctGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAgcggcg                                                                                                                                                  >  1:163661/1‑149 (MQ=255)
       ttCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGcc                                                                                                                                           >  3:250038/1‑149 (MQ=255)
              gCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTg                                                                                                                                    >  5:254704/1‑149 (MQ=255)
               cATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTgg                                                                                                                                   <  8:155318/149‑1 (MQ=255)
                aTGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGt                                                                                                                                  <  7:221249/149‑1 (MQ=255)
                    cctgcctgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTcc                                                                                                                                      >  8:170694/1‑141 (MQ=255)
                    cctgcctgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTcc                                                                                                                                      <  7:170694/141‑1 (MQ=255)
                        cctgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCAtt                                                                                                                            <  6:102937/147‑1 (MQ=255)
                           ggCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATtgctgc                                                                                                                       >  6:87621/1‑149 (MQ=255)
                            gCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATtgctgct                                                                                                                      <  3:332215/149‑1 (MQ=255)
                                  tACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGt                                                                                                                <  4:250038/149‑1 (MQ=255)
                                   aCTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGt                                                                                                                                  <  5:13566/130‑1 (MQ=255)
                                   aCTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGt                                                                                                                                  >  6:13566/1‑130 (MQ=255)
                                     tGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTgc                                                                                                             >  8:93773/1‑149 (MQ=255)
                                     tGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTgc                                                                                                             >  5:13407/1‑149 (MQ=255)
                                     tGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTgc                                                                                                             >  8:273616/1‑149 (MQ=255)
                                       gCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTgcgc                                                                                                           <  1:354041/149‑1 (MQ=255)
                                         tGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCGGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTg                                                                                                         >  7:184072/1‑149 (MQ=255)
                                          gcgcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTgg                                                                                                        <  8:125183/149‑1 (MQ=255)
                                          gcgcgcCGGGCGCGACCCGGAGCCGCCGTTCTCGGGAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCGGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTgg                                                                                                        <  8:301811/149‑1 (MQ=255)
                                           cgcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCCGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTggg                                                                                                       <  4:7526/149‑1 (MQ=255)
                                            gcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGc                                                                                                      >  3:187470/1‑149 (MQ=255)
                                            gcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGc                                                                                                      <  5:109678/149‑1 (MQ=255)
                                            gcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGc                                                                                                      >  6:109678/1‑149 (MQ=255)
                                              gcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGCGCCGCTGccc                                                                                                                                          <  1:311521/111‑2 (MQ=255)
                                              gcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGCGCCGCTGccc                                                                                                                                          >  2:311521/1‑110 (MQ=255)
                                                cGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACt                                                                                                  <  6:135759/149‑1 (MQ=255)
                                                     gcgACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCAccc                                                                                             <  7:93773/149‑1 (MQ=255)
                                                                ccGGTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCgg                                                                                  <  8:269476/149‑1 (MQ=255)
                                                                 cGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCgg                                                                                  >  6:79185/1‑148 (MQ=255)
                                                                    tGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGATCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcat                                                                               >  3:38419/1‑148 (MQ=255)
                                                                        ctcGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagca                                                                          >  6:184167/1‑149 (MQ=255)
                                                                        ctcGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagca                                                                          >  6:302260/1‑149 (MQ=255)
                                                                        ctcGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagca                                                                          >  6:277396/1‑149 (MQ=255)
                                                                         tcGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatga                                                                             <  5:310925/145‑1 (MQ=255)
                                                                         tcGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatga                                                                             >  6:310925/1‑145 (MQ=255)
                                                                          cGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagcat                                                                         >  2:316119/1‑148 (MQ=255)
                                                                           gCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagcatg                                                                        <  3:287244/148‑1 (MQ=255)
                                                                           gCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagca                                                                          <  1:14804/146‑1 (MQ=255)
                                                                            ccAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGAt                                                                      >  3:207024/1‑149 (MQ=255)
                                                                            ccAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGAt                                                                      >  6:294317/1‑149 (MQ=255)
                                                                              aaGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATcc                                                                    >  4:47903/1‑149 (MQ=255)
                                                                                 ccTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGGTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGc                                                                 <  4:187470/149‑1 (MQ=255)
                                                                                             gcgtACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTggcggcggcgg                                                     <  3:47903/149‑1 (MQ=255)
                                                                                                      cGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTg                                             <  5:302260/148‑1 (MQ=255)
                                                                                                        ccATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGcc                                          <  5:87621/149‑1 (MQ=255)
                                                                                                          aTCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagcatga                                                                       >  8:299321/1‑118 (MQ=255)
                                                                                                          aTCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGgcatgagcatga                                                                       <  7:299321/118‑1 (MQ=255)
                                                                                                             ggCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAAt                                      <  1:316119/148‑1 (MQ=255)
                                                                                                             ggCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAAt                                      <  5:184167/148‑1 (MQ=255)
                                                                                                             ggCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGACCATGCTGGCCAAAt                                      <  7:273616/148‑1 (MQ=255)
                                                                                                               cATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCa                                   <  6:254704/149‑1 (MQ=255)
                                                                                                                      ggCAAACTGTTCGCCCGCCTGGGCAGCGGCGCTCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGccc                             >  2:196489/1‑148 (MQ=255)
                                                                                                                       gCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGa                           <  4:207024/149‑1 (MQ=255)
                                                                                                                       gCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGa                           >  6:83707/1‑149 (MQ=255)
                                                                                                                            cTGTGCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCTCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCt                      <  3:278054/149‑1 (MQ=255)
                                                                                                                                 cgcccgccTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATc                  <  6:13407/148‑1 (MQ=255)
                                                                                                                                 cgcccgccTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATc                  >  6:329500/1‑148 (MQ=255)
                                                                                                                                    caggaacggCAGCGCGCCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAAc              >  7:102254/8‑149 (MQ=255)
                                                                                                                                     cgccTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACg             <  5:83707/149‑1 (MQ=255)
                                                                                                                                      gccTGGGCAGCGGCGCGCTGCCGATCCTGGTGGCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGATCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAActg            <  4:38419/149‑3 (MQ=255)
                                                                                                                                         acggCAGCGCGCCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACGGTTc         <  8:102254/147‑1 (MQ=255)
                                                                                                                                              aGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACGGTTCGGCt     >  6:239889/1‑148 (MQ=255)
                                                                                                                                                cggcgCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACGGTTCGGCTATc  <  5:277396/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGGCTACTCACCCGCCCAGGCCGGCATGAGCATGATCCCGCTGGCGGCGGCGGCCATGCTGGCCAAATCCATTGCCCGACCGCTGATCGAACGGTTCGGCTATC  >  NC_002947/5641039‑5641331

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