Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,226,653 T→A 5.1% intergenic (+15/+15) 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,6530TA 5.1% 293.3 / 3.6 97intergenic (+15/+15)PP_4605/PP_4606AraC family transcriptional regulator/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (40/50);  new base A (3/2);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTG  >  NC_002947/5226510‑5226799
                                                                                                                                               |                                                                                                                                                  
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                                                                                                                             >  3:48096/1‑149 (MQ=255)
                  ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTc                                                                                                                             >  5:442042/1‑149 (MQ=255)
                  ggAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCAGTTCCGCAACCGGTTTCGCCAAGCCCTGGGGGCCCCCCCTGGCGAAGAGCGCCAGGCGTTTTTCGTCTGACGCCGGGGCCGCTGGGCGGCCGGATTTGGCTAAAACTc                                                                                                                             >  6:379378/1‑149 (MQ=255)
                    agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCg                                                                                                                            >  5:66/1‑148 (MQ=255)
                    agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGa                                                                                                                           <  8:39150/149‑1 (MQ=255)
                    agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGa                                                                                                                           <  1:378847/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                                                                                                           >  3:302714/1‑148 (MQ=255)
                                     gCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACt                                                                                                           >  2:159520/1‑148 (MQ=255)
                                       gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                        >  4:84934/1‑149 (MQ=255)
                                       gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                        >  1:153755/1‑149 (MQ=255)
                                       gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                        <  6:66372/149‑1 (MQ=255)
                                       gggTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGAc                                                                                                        <  1:158986/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                                                                                                             >  6:376349/1‑131 (MQ=255)
                                                    ccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACaaa                                                                                                             <  5:376349/131‑1 (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                                                      <  4:48096/149‑1 (MQ=255)
                                                                                         ggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGt                                                      <  6:442042/149‑1 (MQ=255)
                                                                                           ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                    >  1:222653/1‑149 (MQ=255)
                                                                                           ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                    >  8:271376/1‑149 (MQ=255)
                                                                                           ccACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                    <  6:114385/149‑1 (MQ=255)
                                                                                                ccTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATacac                                               <  8:258521/149‑1 (MQ=255)
                                                                                                 cTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACAcc                                              >  7:233460/1‑149 (MQ=255)
                                                                                                 cTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACAcc                                              <  8:292970/149‑1 (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                                                    <  1:385008/123‑1 (MQ=255)
                                                                                                                     tttttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                    >  2:385008/1‑123 (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                 <  3:414463/147‑1 (MQ=255)
                                                                                                                                aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                 >  3:422238/1‑147 (MQ=255)
                                                                                                                                aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                 <  4:280780/147‑1 (MQ=255)
                                                                                                                                aaaCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTg                 <  5:388123/147‑1 (MQ=255)
                                                                                                                                 aaatggggCCGCAGAGCGGCCCCGTTTGCCTAAAACTCGAAGGTGCTCGACAACCTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAg              >  1:300409/5‑149 (MQ=255)
                                                                                                                                 aaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACg                         >  3:132061/5‑138 (MQ=255)
                                                                                                                                 aaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACg                         <  4:132061/134‑1 (MQ=255)
                                                                                                                                  aCCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAgg             <  7:128408/149‑1 (MQ=255)
                                                                                                                                     ggggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTGGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAACGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGgcgc          >  3:87489/1‑149 (MQ=255)
                                                                                                                                      gggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCa         >  2:255189/1‑149 (MQ=255)
                                                                                                                                        gCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                   >  7:164643/1‑73 (MQ=255)
                                                                                                                                        gCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                   <  8:164643/73‑1 (MQ=255)
                                                                                                                                           gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtg    >  1:244278/1‑149 (MQ=255)
                                                                                                                                            ctctGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGCTTCTTCACGTTCAGTTGCAGGCGCACGTtgt   <  6:19984/149‑1 (MQ=255)
                                                                                                                                               tGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtgtg  <  2:222653/147‑1 (MQ=255)
                                                                                                                                               tGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtgtg  <  4:3565/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTG  >  NC_002947/5226510‑5226799

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: