Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,699,834 G→C 5.7% P226R (CCC→CGC) ‡ fsr‑II ← fosmidomycin efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,699,8340GC 5.7% 226.3 / 3.6 70P226R (CCC→CGC) ‡fsr‑IIfosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATGG  >  NC_002947/2699688‑2699978
                                                                                                                                                  |                                                                                                                                                
gTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgggg                                                                                                                                                <  6:429558/149‑1 (MQ=255)
 ttGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCggggg                                                                                                                                               >  3:295718/1‑149 (MQ=255)
    gCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCt                                                                                                                                            <  8:16467/149‑1 (MQ=255)
    gCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGACACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCt                                                                                                                                            >  7:179642/1‑149 (MQ=255)
         gctgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGcc                                                                                                                                             >  4:101449/1‑143 (MQ=255)
         gctgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGcc                                                                                                                                             <  3:101449/143‑1 (MQ=255)
         gctgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTgg                                                                                                                                       >  3:303318/1‑149 (MQ=255)
            gcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGCTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgg                                                                                                                                                  >  4:158934/1‑135 (MQ=255)
            gcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGCTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgg                                                                                                                                                  <  3:158934/135‑1 (MQ=255)
                  tCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAg                                                                                                                              >  5:138494/1‑149 (MQ=255)
                  tCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAg                                                                                                                              >  1:346730/1‑149 (MQ=255)
                        cGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCAGGGGGCCTGTTGGCCCGCGAAGGTCtt                                                                                                                         >  8:391436/1‑148 (MQ=255)
                          aaaTGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAg                                                                                                                      <  7:359048/149‑1 (MQ=255)
                             tGCTCGCTCATGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTCGGCCCGCGAAGGTCTTGAGCCg                                                                                                                   <  7:273332/149‑1 (MQ=255)
                               cTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCggg                                                                                                                 <  6:232306/149‑1 (MQ=255)
                                  gATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCt                                                                                                              >  2:153754/1‑149 (MQ=255)
                                   aTCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGACCTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTg                                                                                                             <  4:74917/149‑1 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTgg                                                                                                            >  6:326941/1‑145 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTgg                                                                                                            <  5:326941/145‑1 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCCCCGGCCTGtcgca                                                                                                                                      >  5:106090/1‑115 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCCCCGGCCTGtcgca                                                                                                                                      <  6:106090/119‑5 (MQ=255)
                                           gAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGgc                                                                                                     <  8:179642/149‑1 (MQ=255)
                                             aGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGgcg                                                                                                    >  3:79359/1‑148 (MQ=255)
                                              gttGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCa                                                                                                  <  5:284606/147‑1 (MQ=255)
                                                  aagaagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCt                                                                                              >  2:41602/1‑149 (MQ=255)
                                                    gaagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTcc                                                                                            >  5:9203/1‑149 (MQ=255)
                                                                cGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCAccc                                                                                >  7:354530/1‑149 (MQ=255)
                                                                cGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCAccc                                                                                >  4:196164/1‑149 (MQ=255)
                                                                   aGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCa                                                                             <  5:102361/149‑1 (MQ=255)
                                                                     cGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGt                                                                           <  8:363714/149‑1 (MQ=255)
                                                                         gTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCg                                                                                      >  8:115503/1‑134 (MQ=255)
                                                                         gTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCg                                                                                      <  7:115503/134‑1 (MQ=255)
                                                                          tAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCGGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCa                                                                       <  1:307785/148‑1 (MQ=255)
                                                                          tAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCNCACGCTCCAGCCCGTCACCCGCAGTAGCAc                                                                      <  1:41602/149‑1 (MQ=255)
                                                                          tAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCa                                                                       <  4:20314/148‑1 (MQ=255)
                                                                            gACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCAc                                                                      <  2:346730/147‑1 (MQ=255)
                                                                             aCGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCAc                                                                      >  2:129722/1‑146 (MQ=255)
                                                                             aCGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCAc                                                                      >  1:414996/1‑146 (MQ=255)
                                                                             aCGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAg                                                                   <  5:444560/149‑1 (MQ=255)
                                                                                 aCTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAca                                                               <  1:436280/149‑1 (MQ=255)
                                                                                 aCTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAca                                                               <  4:295718/149‑1 (MQ=255)
                                                                                  cTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCCCCCGCAGTAGCCCCAGAAcag                                                              >  7:358828/1‑149 (MQ=255)
                                                                                     ggCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAcagc                                                             <  3:196164/147‑1 (MQ=255)
                                                                                             tCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGCCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCAGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCcagcgcagcg                                                   >  3:306105/1‑149 (MQ=255)
                                                                                                gCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGcc                                                >  7:220726/1‑149 (MQ=255)
                                                                                                gCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGc                                                 >  4:22470/1‑148 (MQ=255)
                                                                                                 cACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCa                                               <  2:79179/149‑1 (MQ=255)
                                                                                                 cACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCACAGTAGCACCAGAAcagc                                                             <  7:219702/135‑1 (MQ=255)
                                                                                                 cACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCACAGTAGCACCAGAAcagc                                                             >  8:219702/1‑135 (MQ=255)
                                                                                                   cGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAg                                              >  3:333278/1‑148 (MQ=255)
                                                                                                    ggCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCa                                            <  7:391436/149‑1 (MQ=255)
                                                                                                      cgctCACAGCCACAGCGCGCCAGATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGGCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATg                                          <  6:279123/146‑1 (MQ=255)
                                                                                                         tCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAAc                                       >  1:291971/1‑149 (MQ=255)
                                                                                                            ccGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAg                                    >  6:416897/1‑149 (MQ=255)
                                                                                                             cGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAgg                                   <  4:303318/149‑1 (MQ=255)
                                                                                                                  caGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATg                              <  4:79359/149‑1 (MQ=255)
                                                                                                                    gcgcgcCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGc                            <  6:9203/149‑1 (MQ=255)
                                                                                                                    gcgcgcCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGc                            <  8:354530/149‑1 (MQ=255)
                                                                                                                         ccACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGtt                       <  5:101350/149‑1 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                     >  4:360131/1‑149 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                     >  2:240817/1‑149 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                     >  3:279023/1‑149 (MQ=255)
                                                                                                                            caTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGc                    >  4:271768/1‑149 (MQ=255)
                                                                                                                             aTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGcc                   >  8:4376/1‑149 (MQ=255)
                                                                                                                               ggCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGc                            <  3:86892/138‑1 (MQ=255)
                                                                                                                               ggCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGc                            >  4:86892/1‑138 (MQ=255)
                                                                                                                                 cACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGTAGGTCTTGAGGCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGc               <  4:306105/149‑1 (MQ=255)
                                                                                                                                 cACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGc               >  7:412053/1‑149 (MQ=255)
                                                                                                                                      gccaaCAGGCCCCCGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGtt                       <  2:99439/132‑1 (MQ=255)
                                                                                                                                      gccaaCAGGCCCCCGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGtt                       >  1:99439/5‑136 (MQ=255)
                                                                                                                                              gCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCATCCCGTCACCCGCAGTAGCACCAGACCAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATgg  <  7:19209/149‑1 (MQ=255)
                                                                                                                                              gCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATgg  >  1:259771/1‑149 (MQ=255)
                                                                                                                                                cGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGcc                   <  2:52682/130‑1 (MQ=255)
                                                                                                                                                cGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGcc                   >  1:52682/1‑130 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATGG  >  NC_002947/2699688‑2699978

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