Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,226,651 T→A 5.1% intergenic (+13/+17) PP_4605 → / ← PP_4606 AraC family transcriptional regulator/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,226,6510TA 5.1% 305.3 / 4.1 97intergenic (+13/+17)PP_4605/PP_4606AraC family transcriptional regulator/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (42/50);  new base A (3/2);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

GTTGCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGT  >  NC_002947/5226503‑5226798
                                                                                                                                                    |                                                                                                                                                   
gttgCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGct                                                                                                                                                     >  5:128866/1‑149 (MQ=255)
gttgCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGct                                                                                                                                                     >  1:237056/1‑149 (MQ=255)
       cTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGc                                                                                                                                              <  7:401639/149‑1 (MQ=255)
          ggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGcc                                                                                                                                             <  5:428932/147‑1 (MQ=255)
           gAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACATGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGccccg                                                                                                                                          >  1:168789/1‑148 (MQ=255)
            aGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTGTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGccccgg                                                                                                                                         <  4:129068/149‑3 (MQ=255)
               aaTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACTTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                       <  6:363582/148‑1 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTg                                                                                                                                      <  6:165473/148‑1 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGc                                                                                                                                     <  7:305017/149‑1 (MQ=255)
                    cGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTa                                                                                                                                  >  6:101732/1‑148 (MQ=255)
                         ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAAc                                                                                                                              <  2:237056/147‑1 (MQ=255)
                         ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTc                                                                                                                            >  5:19984/1‑149 (MQ=255)
                         ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTc                                                                                                                            >  5:442042/1‑149 (MQ=255)
                         ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTc                                                                                                                            >  3:48096/1‑149 (MQ=255)
                         ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCAGTTCCGCAACCGGTTTCGCCAAGCCCTGGGGGCCCCCCCTGGCGAAGAGCGCCAGGCGTTTTTCGTCTGACGCCGGGGCCGCTGGGCGGCCGGATTTGGCTAAAACTc                                                                                                                            >  6:379378/1‑149 (MQ=255)
                           agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCg                                                                                                                           >  5:66/1‑148 (MQ=255)
                           agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGa                                                                                                                          <  1:378847/149‑1 (MQ=255)
                           agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGa                                                                                                                          <  8:39150/149‑1 (MQ=255)
                            gagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGaa                                                                                                                         <  5:379378/149‑1 (MQ=255)
                             agaTTTCCTTATTGTGCGGGTTTGACAGTGCCAGCCATTTCATCAACCGGTTTCGCCAAGCACTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAg                                                                                                                        <  4:265284/149‑1 (MQ=255)
                             agaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGaa                                                                                                                         <  6:290472/148‑1 (MQ=255)
                                 tGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGc                                                                                                                    >  7:82209/1‑149 (MQ=255)
                                           tgCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAc                                                                                                           <  2:412089/148‑1 (MQ=255)
                                           tgCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCCGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAc                                                                                                           <  6:212909/148‑1 (MQ=255)
                                            gCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACt                                                                                                          >  2:159520/1‑148 (MQ=255)
                                            gCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACt                                                                                                          >  3:302714/1‑148 (MQ=255)
                                              gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                       >  1:153755/1‑149 (MQ=255)
                                              gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                       >  4:84934/1‑149 (MQ=255)
                                              gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                       <  1:158986/149‑1 (MQ=255)
                                              gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                       <  6:66372/149‑1 (MQ=255)
                                                gTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCt                                                                                                     <  4:7290/149‑1 (MQ=255)
                                                gTTTGCCAGTGCCAGCCATTTCAGCAAACGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACAGGGGCCGCTCTTCGGCCCCAATTGCCTAAAAATCGAAGGTGCAGCACAAACTTACCt                                                                                                     >  4:4240/1‑149 (MQ=255)
                                                 tttGCCAGTGGCAGCCATTTCAGCAACGTGTTTCGCCAAGCCGTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCGTAAAACTCGAAGGTGCTCGACAAACTGACCTg                                                                                                    <  4:302714/149‑1 (MQ=255)
                                                 tttGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCGGCGGCCCCATTTGCCTAAAACTCGAAGTTGGTCGACAAACTGACCTg                                                                                                    >  3:280780/1‑149 (MQ=255)
                                                    gccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGcc                                                                                                  <  5:18586/148‑1 (MQ=255)
                                                      cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcgcg                                                                                               >  7:294798/1‑149 (MQ=255)
                                                           ccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACaaa                                                                                                            <  5:376349/131‑1 (MQ=255)
                                                           ccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACaaa                                                                                                            >  6:376349/1‑131 (MQ=255)
                                                           ccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcg                                                                                          >  3:407366/1‑149 (MQ=255)
                                                             agccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgc                                                                                         >  7:292970/1‑148 (MQ=255)
                                                                  tttCAGCAACCGGTTTCGCCACGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcg                                                                                   >  5:114385/1‑149 (MQ=255)
                                                                   ttCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                  >  4:414463/1‑149 (MQ=255)
                                                                     cAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGAGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCAGCCGCGcatcgcccatcgcc                                                                                 >  4:60892/1‑148 (MQ=255)
                                                                      aGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgcca                                                                                >  1:205941/1‑148 (MQ=255)
                                                                       gCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCGCCATTTGCCTAAAACTCGAAGGTGCTCGACAACCTGACCTGCCGCGCATCGCCAATCGCcac                                                                               >  3:127535/1‑148 (MQ=255)
                                                                        cAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAACCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcaca                                                                              >  7:239918/1‑148 (MQ=255)
                                                                        cAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcaca                                                                              >  7:258521/1‑148 (MQ=255)
                                                                             ggTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTAc                                                                        <  5:282831/149‑1 (MQ=255)
                                                                              gTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACt                                                                       >  3:121215/1‑149 (MQ=255)
                                                                                 tCGCCAAGCCCTGGGGGCCACACCTGGCTAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTgg                                                                     <  3:235028/148‑1 (MQ=255)
                                                                                 tCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTgg                                                                     <  5:101732/148‑1 (MQ=255)
                                                                                     cAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                >  5:260650/1‑149 (MQ=255)
                                                                                         cccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAccgcc                                                            >  8:15479/1‑149 (MQ=255)
                                                                                              ggggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCt                                                       >  6:388123/1‑149 (MQ=255)
                                                                                              ggggCCACACCTCGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAg                                                         <  6:128866/147‑1 (MQ=255)
                                                                                               gggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTg                                                      <  5:398295/149‑1 (MQ=255)
                                                                                                ggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGt                                                     <  6:442042/149‑1 (MQ=255)
                                                                                                ggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGt                                                     <  4:48096/149‑1 (MQ=255)
                                                                                                  ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                   >  8:271376/1‑149 (MQ=255)
                                                                                                  ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                   >  1:222653/1‑149 (MQ=255)
                                                                                                  ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                   <  6:114385/149‑1 (MQ=255)
                                                                                                       ccTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATacac                                              <  8:258521/149‑1 (MQ=255)
                                                                                                        cTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACAcc                                             <  8:292970/149‑1 (MQ=255)
                                                                                                        cTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACAcc                                             >  7:233460/1‑149 (MQ=255)
                                                                                                           gCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTg                                          >  1:257292/1‑149 (MQ=255)
                                                                                                                tACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGaa                                     <  3:84934/149‑1 (MQ=255)
                                                                                                                  ccGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGaa                                     >  8:128408/1‑147 (MQ=255)
                                                                                                                  ccGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAAc                                    <  8:294798/148‑1 (MQ=255)
                                                                                                                   cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGt                                                     >  7:377986/1‑130 (MQ=255)
                                                                                                                   cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGt                                                     <  8:377986/130‑1 (MQ=255)
                                                                                                                     cAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGt                                <  2:257292/149‑1 (MQ=255)
                                                                                                                      aaGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAccgcc                                                            <  4:316964/120‑1 (MQ=255)
                                                                                                                          gcTTTATCATCTGAAACCGGCGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                           <  4:18726/148‑1 (MQ=255)
                                                                                                                          ccTTTTTCATCTGAAACCGGGGCCGCTCTGCTGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                           <  2:60302/149‑1 (MQ=255)
                                                                                                                          ccTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                           <  8:233460/149‑1 (MQ=255)
                                                                                                                          ccTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                           <  2:153755/149‑1 (MQ=255)
                                                                                                                           cTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCa                          >  4:181189/1‑149 (MQ=255)
                                                                                                                            tttttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                   >  2:385008/1‑123 (MQ=255)
                                                                                                                            tttttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                   <  1:385008/123‑1 (MQ=255)
                                                                                                                              tttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCAc                         <  4:407366/147‑1 (MQ=255)
                                                                                                                               ttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCGAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGt                       <  3:60892/148‑1 (MQ=255)
                                                                                                                               ttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGt                       >  8:214156/1‑148 (MQ=255)
                                                                                                                                tCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGtt                      <  6:260650/148‑1 (MQ=255)
                                                                                                                                  aTCTGAAACCGTGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAg                   <  1:189910/149‑1 (MQ=255)
                                                                                                                                  aTCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCa                    <  8:82209/148‑1 (MQ=255)
                                                                                                                                  aTCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAg                   <  4:121215/149‑1 (MQ=255)
                                                                                                                                       aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                <  4:280780/147‑1 (MQ=255)
                                                                                                                                       aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                <  3:414463/147‑1 (MQ=255)
                                                                                                                                       aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                <  5:388123/147‑1 (MQ=255)
                                                                                                                                       aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                >  3:422238/1‑147 (MQ=255)
                                                                                                                                        aaatggggCCGCAGAGCGGCCCCGTTTGCCTAAAACTCGAAGGTGCTCGACAACCTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAg             >  1:300409/5‑149 (MQ=255)
                                                                                                                                        aaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACg                        <  4:132061/134‑1 (MQ=255)
                                                                                                                                        aaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACg                        >  3:132061/5‑138 (MQ=255)
                                                                                                                                         aCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAgg            <  7:128408/149‑1 (MQ=255)
                                                                                                                                            ggggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTGGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAACGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGgcgc         >  3:87489/1‑149 (MQ=255)
                                                                                                                                             gggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCa        >  2:255189/1‑149 (MQ=255)
                                                                                                                                               gCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                  <  8:164643/73‑1 (MQ=255)
                                                                                                                                               gCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                  >  7:164643/1‑73 (MQ=255)
                                                                                                                                                  gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtg   >  1:244278/1‑149 (MQ=255)
                                                                                                                                                   ctctGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGCTTCTTCACGTTCAGTTGCAGGCGCACGTtgt  <  6:19984/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGT  >  NC_002947/5226503‑5226798

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