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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,699,8350GC 5.7% 221.2 / 3.7 70P226A (CCC→GCC) ‡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 = 7.12e-01

GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATGGCAG  >  NC_002947/2699688‑2699981
                                                                                                                                                   |                                                                                                                                                  
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)
                  tCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAg                                                                                                                                 >  1:346730/1‑149 (MQ=255)
                  tCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAg                                                                                                                                 >  5:138494/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                                                                                                               <  5:326941/145‑1 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTgg                                                                                                               >  6:326941/1‑145 (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                                                                                         <  7:115503/134‑1 (MQ=255)
                                                                         gTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCg                                                                                         >  8:115503/1‑134 (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                                                                         >  1:414996/1‑146 (MQ=255)
                                                                             aCGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCAc                                                                         >  2:129722/1‑146 (MQ=255)
                                                                             aCGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAg                                                                      <  5:444560/149‑1 (MQ=255)
                                                                                 aCTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAca                                                                  <  4:295718/149‑1 (MQ=255)
                                                                                 aCTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAca                                                                  <  1:436280/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                                                                >  8:219702/1‑135 (MQ=255)
                                                                                                 cACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCACAGTAGCACCAGAAcagc                                                                <  7:219702/135‑1 (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                               <  8:354530/149‑1 (MQ=255)
                                                                                                                    gcgcgcCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGc                               <  6:9203/149‑1 (MQ=255)
                                                                                                                         ccACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGtt                          <  5:101350/149‑1 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                        >  2:240817/1‑149 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                        >  3:279023/1‑149 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                        >  4:360131/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)
                                                                                                                                                   gggCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATGgcag  >  5:61443/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAATGGCAG  >  NC_002947/2699688‑2699981

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