Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,367,110 G→A 5.3% intergenic (‑93/+17) pmrD ← / ← menE inactive two‑component system connector protein/O‑succinylbenzoate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,367,1100GA 5.3% 1179.4 / 32.0 322intergenic (‑93/+17)pmrD/menEinactive two‑component system connector protein/O‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (171/134);  new base A (10/7);  total (181/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

AGTACATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCACCCATTC  >  CP009273/2366962‑2367258
                                                                                                                                                    |                                                                                                                                                    
aGTACATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAg                                                                                                                                                      <  1:3934631/149‑1 (MQ=255)
  tACATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTg                                                                                                                                                    >  1:142722/1‑149 (MQ=255)
   aCATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGa                                                                                                                                                   >  1:3219161/1‑149 (MQ=255)
    cATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGaa                                                                                                                                                  >  2:4732550/1‑149 (MQ=255)
    cATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAATGaa                                                                                                                                                  >  2:2167131/1‑149 (MQ=255)
     aTGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGaaa                                                                                                                                                 >  2:852317/1‑149 (MQ=255)
     aTGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGaaa                                                                                                                                                 >  2:2798962/1‑149 (MQ=255)
     aTGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGaaa                                                                                                                                                 <  2:700100/149‑1 (MQ=255)
     aTGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGaaa                                                                                                                                                 <  1:119227/149‑1 (MQ=255)
        tCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACgg                                                                                                                                              <  2:275697/149‑1 (MQ=255)
           gTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAgg                                                                                                                                           >  1:760803/1‑149 (MQ=255)
            ttATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAATGAAACGGAggg                                                                                                                                          >  2:3477138/1‑149 (MQ=255)
            ttATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAAAAGATAAACTCTTTATAAACGGCAATGGAATGGCAGAAATGAGAAGTGAAAAGGAggg                                                                                                                                          >  2:3695912/1‑149 (MQ=255)
            ttATCGTGTTTAGTGCAGCACGATTTTTTGACCAGCCATTCCAGTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAggg                                                                                                                                          <  2:3406229/149‑1 (MQ=255)
             tATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGATTTTTCACATCAATTTACAGATTAACAGTTAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGt                                                                                                                                         <  2:4803873/149‑1 (MQ=255)
              aTCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGACGTGAAACGGAGGGtt                                                                                                                                        >  1:889057/1‑149 (MQ=255)
               tCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTa                                                                                                                                       >  1:985870/1‑149 (MQ=255)
               tCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATAAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTa                                                                                                                                       >  1:3287936/1‑149 (MQ=255)
                cTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTaa                                                                                                                                      <  2:1623646/149‑1 (MQ=255)
                  tGTTTATTGCAGCACTATATTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTACACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAatt                                                                                                                                    <  2:3139386/149‑1 (MQ=255)
                  tGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAatt                                                                                                                                    <  2:2691920/149‑1 (MQ=255)
                  tGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAatt                                                                                                                                    >  2:3482289/1‑149 (MQ=255)
                  tGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAatt                                                                                                                                    >  2:4428601/1‑149 (MQ=255)
                  tGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAatt                                                                                                                                    >  1:4097172/1‑149 (MQ=255)
                    tttATTGCAGCACGATTTTTTGACCAGCCATTACATTGCATTATCCTGTTTGCTAAGAGTTTTAACATAAATTTACATATTAACATATAAACTCTTTATAAACGGCAATGGATTGGCAGAAATGAGAAGTGAAACGGAGGGTTAattat                                                                                                                                  <  1:3808158/149‑1 (MQ=255)
                     ttATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAATGAAACGGAGGGTTAattatt                                                                                                                                 <  2:1401335/149‑1 (MQ=255)
                     ttATTGCAGCACGATTTTTTGACCAGCCATTCCAGTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAattatt                                                                                                                                 <  1:1892004/149‑1 (MQ=255)
                      tATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTg                                                                                                                                <  2:2514469/149‑1 (MQ=255)
                       aTTGCCGCACGATTTTTTGACCAGCCTTTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGc                                                                                                                               <  2:1941051/149‑1 (MQ=255)
                         tGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCAAATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGAAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTg                                                                                                                             <  1:488676/149‑1 (MQ=255)
                           cagcaCTATTTTTAGACCAGCCATTCCATTGCATTATCCTGTTTGCTAATAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGAc                                                                                                                           <  2:4644449/149‑1 (MQ=255)
                           cagcaCGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGAc                                                                                                                           <  1:4596595/149‑1 (MQ=255)
                                  aTTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAc                                                                                                                    >  2:1638085/1‑149 (MQ=255)
                                  aTTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAc                                                                                                                    <  1:207429/149‑1 (MQ=255)
                                   ttttttGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAcc                                                                                                                   >  1:3110767/1‑149 (MQ=255)
                                   ttttttGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCAcc                                                                                                                   >  2:4303809/1‑149 (MQ=255)
                                        tGACCAGCAATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCActc                                                                                                              <  2:4357131/149‑1 (MQ=255)
                                         gACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCActct                                                                                                             <  1:2139044/149‑1 (MQ=255)
                                         gACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCActct                                                                                                             <  1:465554/149‑1 (MQ=255)
                                         gACCACCAATTCCATTGCATTATCATGTTTGCTAAGAGTTTACACATCAATTTACAGATTAACATATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCActct                                                                                                             <  1:4960264/149‑1 (MQ=255)
                                          aCCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGACTGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCtt                                                                                                            >  1:4420049/1‑149 (MQ=255)
                                          aCCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGACATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCtt                                                                                                            >  1:1825561/1‑149 (MQ=255)
                                                cattccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgtg                                                                                                      >  1:2583369/1‑147 (MQ=255)
                                                cattccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgcg                                                                                                      >  1:4984569/1‑149 (MQ=255)
                                                 attccattGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAgcgc                                                                                                     >  1:1579725/1‑149 (MQ=255)
                                                       ttGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACg                                                                                               >  1:62959/1‑149 (MQ=255)
                                                           aTTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTCATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGaa                                                                                           >  2:5112691/1‑149 (MQ=255)
                                                            ttATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGaaa                                                                                          <  2:1282510/149‑1 (MQ=255)
                                                              aTCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAtt                                                                                        >  1:1013126/1‑149 (MQ=255)
                                                                ccTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAAtttt                                                                                      >  2:289060/1‑149 (MQ=255)
                                                                   gTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAAt                                                                                   <  1:2798962/149‑1 (MQ=255)
                                                                        cTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccgc                                                                              >  1:4054122/1‑149 (MQ=255)
                                                                        cTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccgc                                                                              >  1:3637764/1‑149 (MQ=255)
                                                                        cTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccgc                                                                              >  1:3648140/1‑149 (MQ=255)
                                                                          aaGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATAccgccg                                                                            >  1:3726280/1‑149 (MQ=255)
                                                                            gagTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtt                                                                          >  1:3112940/1‑149 (MQ=255)
                                                                            gagTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAAATATTGCTGACGATGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtt                                                                          >  1:3024734/1‑149 (MQ=255)
                                                                            gagTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAATGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtt                                                                          <  2:3222103/149‑1 (MQ=255)
                                                                              gTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtttt                                                                        <  2:3581149/149‑1 (MQ=255)
                                                                              gTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtttt                                                                        >  1:4490923/1‑149 (MQ=255)
                                                                              gTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGtttt                                                                        <  2:53350/149‑1 (MQ=255)
                                                                                   cacaTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGc                                                                   >  1:4258442/1‑149 (MQ=255)
                                                                                      aTCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTcc                                                                >  1:137087/1‑149 (MQ=255)
                                                                                         aaTTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCcggc                                                             >  1:2872174/1‑149 (MQ=255)
                                                                                          aTTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCcggcg                                                            >  2:3366818/1‑149 (MQ=255)
                                                                                           tttACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCcggcgg                                                           >  1:2195876/1‑149 (MQ=255)
                                                                                              aCAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCag                                                        <  1:5096735/149‑1 (MQ=255)
                                                                                                    tAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACAGCCGTTTTTCAGCTCAGGCGGCAGAGTTAg                                                  >  2:4258273/1‑149 (MQ=255)
                                                                                                     aaCAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGc                                                 <  2:1102993/149‑1 (MQ=255)
                                                                                                        aGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAg                                              >  1:2415639/1‑149 (MQ=255)
                                                                                                         gATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAgc                                             >  1:900462/1‑149 (MQ=255)
                                                                                                            aaaCTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCa                                          >  2:2964284/1‑149 (MQ=255)
                                                                                                              aCTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCAcc                                        >  2:3849879/1‑149 (MQ=255)
                                                                                                              aCTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCAcc                                        <  2:3103421/149‑1 (MQ=255)
                                                                                                              aCTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCAcc                                        <  2:3266913/149‑1 (MQ=255)
                                                                                                                tctTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCgg                                      >  2:2357926/1‑149 (MQ=255)
                                                                                                                 ctTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgt                                     <  2:2429426/149‑1 (MQ=255)
                                                                                                                 ctTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgt                                     <  2:4962248/149‑1 (MQ=255)
                                                                                                                  tttATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgtt                                    >  2:3971016/1‑149 (MQ=255)
                                                                                                                     aTCAACGGCATTGGAATGGAAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgttgtt                                 <  1:473065/149‑1 (MQ=255)
                                                                                                                     aTCAACGGCAATGGAATGGCAGAACTGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgttgtt                                 <  2:3518239/149‑1 (MQ=255)
                                                                                                                     aTCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGgttgtt                                 <  2:4144467/149‑1 (MQ=255)
                                                                                                                         aCGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGaaa                             >  1:3830917/1‑149 (MQ=255)
                                                                                                                           ggCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAAc                           <  1:3000567/149‑1 (MQ=255)
                                                                                                                                   aatggCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGc                   >  2:260922/1‑149 (MQ=255)
                                                                                                                                     tggCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGTGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCtt                 <  1:852220/149‑1 (MQ=255)
                                                                                                                                     tggCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCCCCGGTTGTTGAAAACGGGCCAGCtt                 >  1:2299495/1‑149 (MQ=255)
                                                                                                                                      ggCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTa                >  2:2979628/1‑149 (MQ=255)
                                                                                                                                      ggCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTa                <  1:4386522/149‑1 (MQ=255)
                                                                                                                                       gCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGGTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTAt               >  1:4934141/1‑149 (MQ=255)
                                                                                                                                        cAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATc              <  2:1860601/149‑1 (MQ=255)
                                                                                                                                           aaaTGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCttt           <  2:2279720/149‑1 (MQ=255)
                                                                                                                                            aaTGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTc          >  2:258435/1‑149 (MQ=255)
                                                                                                                                            aaTGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTc          >  1:1476702/1‑149 (MQ=255)
                                                                                                                                             aTGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCa         >  1:1075148/1‑149 (MQ=255)
                                                                                                                                                agaAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCAccc      <  2:4097172/149‑1 (MQ=255)
                                                                                                                                                  aaGTGAAACGGAGGGTTTATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCCGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCACCCAt    <  2:103914/149‑1 (MQ=255)
                                                                                                                                                   aGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGATAACGGGCCAGCTTATCTTTCACCCAtt   <  2:2091972/149‑1 (MQ=255)
                                                                                                                                                   aGTGAAACAGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCACCCAtt   >  1:1690692/1‑149 (MQ=255)
                                                                                                                                                    gTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCACCCATTc  >  1:1275326/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTACATGTCTGTTATCTTGTTTATTGCAGCACGATTTTTTGACCAGCCATTCCATTGCATTATCCTGTTTGCTAAGAGTTTTCACATCAATTTACAGATTAACAGATAAACTCTTTATCAACGGCAATGGAATGGCAGAAATGAGAAGTGAAACGGAGGGTTAATTATTGCTGACGTTGCACCCACTCTTTTAGCGCCTGACGTGAAATTTTAATACCGCCGTTTTTCAGCTCCGGCGGCAGAGTTAGCCAGCGCACCGGTTGTTGAAAACGGGCCAGCTTATCTTTCACCCATTC  >  CP009273/2366962‑2367258

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