Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,362 G→C 10.3% intergenic (+32/+62) ycgR → / ← PP_4398 flagellar brake protein YcgR/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,988,3620GC10.3% 255.4 / 16.6 87intergenic (+32/+62)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (37/41);  new base C (4/5);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

ACCACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCC  >  NC_002947/4988214‑4988508
                                                                                                                                                    |                                                                                                                                                   
accaccTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACctg                                                                                                                                                     >  3:333725/1‑149 (MQ=255)
     cTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAg                                                                                                                                               >  5:57745/1‑149 (MQ=255)
            ggCGTGCGTTTCCACAACCTCAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGAGAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGt                                                                                                                                        <  3:25847/149‑1 (MQ=255)
             gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                       <  6:57745/149‑1 (MQ=255)
             gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                       <  8:243891/149‑1 (MQ=255)
             gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                       <  2:264447/149‑1 (MQ=255)
             gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                       >  7:243891/1‑149 (MQ=255)
                tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCTTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTc                                                                                                                                    <  1:68265/149‑1 (MQ=255)
                tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCATCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTTTTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTc                                                                                                                                    <  7:157440/149‑1 (MQ=255)
                tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTc                                                                                                                                    <  4:232704/149‑1 (MQ=255)
                  cgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAt                                                                                                                                  <  2:183973/149‑1 (MQ=255)
                    tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCTGTAAAACAGAAAACGCCACCTGCTGAg                                                                                                                                               <  4:219678/134‑1 (MQ=255)
                    tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAg                                                                                                                                               >  3:219678/1‑134 (MQ=255)
                      tCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTctgttt                                                                                                                               <  1:114686/148‑6 (MQ=39)
                       ccACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgtg                                                                                                                             >  8:101192/1‑149 (MQ=255)
                       ccACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgt                                                                                                                              <  3:298762/148‑1 (MQ=255)
                       ccACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgt                                                                                                                              >  4:298762/1‑148 (MQ=255)
                        cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtttt                                                                                                                                     >  7:204874/1‑140 (MQ=255)
                        cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtttt                                                                                                                                     <  8:204874/140‑1 (MQ=255)
                        cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGtttt                                                                                                                                     >  3:118177/1‑140 (MQ=255)
                        cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGtttt                                                                                                                                     <  4:118177/140‑1 (MQ=255)
                           aaCCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTc                                                                                                                         <  4:224648/149‑1 (MQ=255)
                           aaCCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTc                                                                                                                         >  6:176401/1‑149 (MQ=255)
                              cTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctgc                                                                                                                      >  1:215441/1‑149 (MQ=255)
                                     ggCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTgctgggctg                                                                                                               <  2:215441/149‑1 (MQ=255)
                                     ggCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTgctgggctg                                                                                                               >  3:224208/1‑149 (MQ=255)
                                           cgcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGACGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCt                                                                                                         >  1:113037/1‑149 (MQ=255)
                                            gcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCtt                                                                                                        >  5:128996/1‑149 (MQ=255)
                                            gcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCtt                                                                                                        <  7:101192/149‑1 (MQ=255)
                                               gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGGTCGACAAGGACGATTATTGATCGGCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcgc                                                                                                     >  1:266493/1‑149 (MQ=255)
                                                cgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcgcg                                                                                                    >  6:336355/1‑149 (MQ=255)
                                                      aTCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTaaa                                                                                              >  7:106260/1‑149 (MQ=255)
                                                       tCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTaaa                                                                                              >  6:337097/1‑148 (MQ=255)
                                                        cGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAcc                                                                                            >  2:228760/1‑149 (MQ=255)
                                                          agagTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTAGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGAGGCGTTTTCATTTGTGTTTCTGCTGGGCTGTCCGCTTCGCGGGTAAACCCg                                                                                          >  7:132857/1‑149 (MQ=255)
                                                                       ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                   <  3:335016/127‑1 (MQ=255)
                                                                       ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                   >  4:335016/1‑127 (MQ=255)
                                                                       ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGAc                                                                             <  1:228760/149‑1 (MQ=255)
                                                                          gCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTg                                                                          >  6:239012/1‑149 (MQ=255)
                                                                           cTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGc                                                                         >  1:10784/1‑149 (MQ=255)
                                                                              cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCt                                                                        <  6:128996/147‑1 (MQ=255)
                                                                              cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCt                                                                        <  2:58051/147‑1 (MQ=255)
                                                                              cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTg                                                                       <  4:224208/148‑1 (MQ=255)
                                                                               aGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                     <  2:10784/149‑1 (MQ=255)
                                                                               aGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                     <  4:121200/149‑1 (MQ=255)
                                                                                     aGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTc                                                               <  6:139665/149‑1 (MQ=255)
                                                                                      ggCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTc                                                               <  6:307/148‑1 (MQ=255)
                                                                                       gcgcGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAg                                                             >  3:153418/1‑149 (MQ=255)
                                                                                       gcgcGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAg                                                             <  4:162110/149‑1 (MQ=255)
                                                                                         gcggcgGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgc                                                           <  5:155/149‑1 (MQ=255)
                                                                                             cggTTCGACAAGGAGGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCt                                                       <  5:337097/149‑1 (MQ=255)
                                                                                              ggTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTg                                                      >  8:288970/1‑149 (MQ=255)
                                                                                               gTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGt                                                     <  8:106260/149‑1 (MQ=255)
                                                                                                    aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                >  3:194219/1‑149 (MQ=255)
                                                                                                    aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                >  1:238495/1‑149 (MQ=255)
                                                                                                    aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                >  3:31388/1‑149 (MQ=255)
                                                                                                        ggACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcg                                            <  5:336355/149‑1 (MQ=255)
                                                                                                         gACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTTGGATGGCCGCTTCGCGGGTAACCCCGCTCCTACACGGACTAGCTGACCTCTTCAGGCGCCTGTtgctg                                                >  8:252395/1‑144 (MQ=255)
                                                                                                         gACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                <  7:252395/144‑1 (MQ=255)
                                                                                                         gACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                           <  6:245332/149‑1 (MQ=255)
                                                                                                            gATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGccc                                         >  1:254754/1‑149 (MQ=255)
                                                                                                                 ttGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTGATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGAGTTGCTGACCTCTTCAGGCGCCTGTTGCTGGGCGCCCTtctc                                    <  4:333725/148‑1 (MQ=255)
                                                                                                                 ttGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                   >  3:155792/1‑149 (MQ=255)
                                                                                                                 ttGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                   >  3:150464/1‑149 (MQ=255)
                                                                                                                       gtcggtcgGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTcc                                                                                      <  3:85729/92‑1 (MQ=255)
                                                                                                                       gtcggtcgGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTcc                                                                                      >  4:85729/1‑92 (MQ=255)
                                                                                                                        tcggtcggTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                   <  1:339316/126‑1 (MQ=255)
                                                                                                                        tcggtcggTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                   >  2:339316/1‑126 (MQ=255)
                                                                                                                        tcggtcggTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCt                            <  5:239012/149‑1 (MQ=255)
                                                                                                                            tcggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgt                        <  5:176401/149‑1 (MQ=255)
                                                                                                                            tcggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgt                        <  4:31388/149‑1 (MQ=255)
                                                                                                                               gtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTc                     >  8:176502/1‑149 (MQ=255)
                                                                                                                                cacaaatgAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCt                    >  6:351580/8‑149 (MQ=39)
                                                                                                                                 aaaaCAGAAAACGCCACCTGCTGAGGTTGCGTATTCATTTGTGTTTCTGCTGGGCTGGACGCTTCGCGTGTAAACCCCCTCCAAAACCGACTTGCTGACCACTTCCTGCGCATGTTGCTGCGCGACCTTGTATTCATACTCGCTGTCt                    >  6:226237/1‑148 (MQ=255)
                                                                                                                                  aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAccc                                                                                           >  6:247685/1‑76 (MQ=255)
                                                                                                                                  aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAccc                                                                                           <  5:247685/76‑1 (MQ=255)
                                                                                                                                  aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                           <  7:100178/124‑1 (MQ=255)
                                                                                                                                  aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                           >  8:100178/1‑124 (MQ=255)
                                                                                                                                   aaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCa                 <  4:150464/149‑1 (MQ=255)
                                                                                                                                   aaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCa                 <  2:113037/149‑1 (MQ=255)
                                                                                                                                   aaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCa                 <  4:155792/149‑1 (MQ=255)
                                                                                                                                      aGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgcc              <  7:288970/149‑1 (MQ=255)
                                                                                                                                        aaaaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgccgc            >  3:69415/1‑149 (MQ=255)
                                                                                                                                          aaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCa          >  3:246066/1‑149 (MQ=255)
                                                                                                                                          aaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                   >  8:121313/1‑124 (MQ=255)
                                                                                                                                          aaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                   <  7:121313/124‑1 (MQ=255)
                                                                                                                                           aCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCAc         >  7:300679/1‑149 (MQ=255)
                                                                                                                                                 ccTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGgcgc   <  2:238495/149‑1 (MQ=255)
                                                                                                                                                  ctgctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGcc  <  8:223777/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCC  >  NC_002947/4988214‑4988508

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: