Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,363 C→A 8.3% intergenic (+33/+61) 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,3630CA 8.3% 185.4 / 6.2 60intergenic (+33/+61)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base A (2/3);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCG  >  NC_002947/4988222‑4988506
                                                                                                                                             |                                                                                                                                               
tGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAgg                                                                                                                                            <  2:226348/147‑1 (MQ=255)
    ggCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGt                                                                                                                                      >  7:281513/1‑149 (MQ=255)
     gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                     >  3:162195/1‑149 (MQ=255)
     gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                     <  3:280270/149‑1 (MQ=255)
      cgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGttt                                                                                                                                    <  8:23020/149‑1 (MQ=255)
        tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTc                                                                                                                                  <  4:230932/149‑1 (MQ=255)
        tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTc                                                                                                                                  <  1:195985/149‑1 (MQ=255)
          cgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCa                                                                                                                                 <  7:31996/148‑1 (MQ=255)
            tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                              <  5:64774/149‑1 (MQ=255)
                cacaGCCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgtg                                                                                                                           <  1:82016/148‑1 (MQ=255)
                 acaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGtt                                                                                                                         >  1:248593/1‑149 (MQ=255)
                    aCCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTct                                                                                                                      <  7:58000/149‑1 (MQ=255)
                     ccTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTct                                                                                                                      >  5:50636/1‑148 (MQ=255)
                             ggCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTgctgggctg                                                                                                             >  3:69057/1‑149 (MQ=255)
                              gCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctg                                                                                                                     <  5:24595/140‑1 (MQ=255)
                              gCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTACGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctg                                                                                                                     >  6:24595/1‑140 (MQ=255)
                                       gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcgc                                                                                                   >  4:143930/1‑149 (MQ=255)
                                        cgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcgcg                                                                                                  <  3:143930/149‑1 (MQ=255)
                                         gcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                 <  2:248593/149‑1 (MQ=255)
                                         gcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                 <  2:132367/149‑1 (MQ=255)
                                          caAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCggg                                                                                                <  4:144888/149‑1 (MQ=255)
                                           aaGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCggg                                                                                                >  2:203951/1‑148 (MQ=255)
                                              aTCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTaaa                                                                                            >  4:237471/1‑149 (MQ=255)
                                               tCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGAATATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAc                                                                                           >  3:258732/1‑149 (MQ=255)
                                                     gTTTCGTGTACCAGCTGGATCGTGATGCGCGGCGGTACGAAAAGGACGTTTATTGTTCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCA‑TTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                     <  2:277089/148‑1 (MQ=255)
                                                           tgtACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTacac                                                                               <  1:4802/149‑1 (MQ=255)
                                                            gtACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAcc                                                                              >  6:239186/1‑149 (MQ=255)
                                                             tACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAcc                                                                              <  4:162195/148‑1 (MQ=255)
                                                                       aGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                   <  8:120647/149‑1 (MQ=255)
                                                                        gCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAcc                                                                  >  5:31320/1‑149 (MQ=255)
                                                                               gcgcGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAg                                                           <  4:69057/149‑1 (MQ=255)
                                                                                 gcggcgGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgc                                                         <  2:290944/149‑1 (MQ=255)
                                                                                       gTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGt                                                   >  2:38886/1‑149 (MQ=255)
                                                                                         tCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACct                                                                 >  4:217861/1‑133 (MQ=255)
                                                                                         tCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACct                                                                 <  3:217861/133‑1 (MQ=255)
                                                                                          cGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAACTCTTGATGCGCCTGTtg                                                 <  6:50636/148‑1 (MQ=255)
                                                                                            aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgcgc                                                       <  5:115402/140‑1 (MQ=255)
                                                                                            aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgcgc                                                       >  6:115402/1‑140 (MQ=255)
                                                                                             cAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                     >  5:14749/1‑109 (MQ=255)
                                                                                             cAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                     <  6:14749/109‑1 (MQ=255)
                                                                                             cAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                             <  8:40340/149‑1 (MQ=255)
                                                                                                  aCGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGcc                                        <  1:38886/149‑1 (MQ=255)
                                                                                                    gATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCt                                      >  6:224560/1‑149 (MQ=255)
                                                                                                      ttattGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTtc                                    <  1:203951/149‑1 (MQ=255)
                                                                                                             tcgtcgGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCAt                             >  7:49965/1‑149 (MQ=255)
                                                                                                                  ggtcggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                        <  8:49965/149‑1 (MQ=255)
                                                                                                                  agccggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                        <  3:221676/146‑1 (MQ=255)
                                                                                                                      ggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgtgt                    <  4:258732/149‑1 (MQ=255)
                                                                                                                               gAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                             <  8:91471/115‑1 (MQ=39)
                                                                                                                               gAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                             >  7:91471/1‑115 (MQ=39)
                                                                                                                                aaaaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCAGCCTCGGTGTCTTCAgccgc          >  4:1235/1‑149 (MQ=255)
                                                                                                                                 aaaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgccgc          <  5:224560/148‑1 (MQ=255)
                                                                                                                                  aaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCa        >  2:157551/1‑149 (MQ=255)
                                                                                                                                  aaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCa        >  3:284689/1‑149 (MQ=255)
                                                                                                                                  aaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACctc                                                                >  3:227810/1‑93 (MQ=37)
                                                                                                                                  aaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACctc                                                                <  4:227810/93‑1 (MQ=37)
                                                                                                                                   aCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCAc       <  5:239186/149‑1 (MQ=255)
                                                                                                                                      ccACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAg              >  2:279740/1‑139 (MQ=255)
                                                                                                                                      ccACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAg              <  1:279740/139‑1 (MQ=255)
                                                                                                                                        aCCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCTCGACCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGgcg  >  1:185544/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                               
TGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCG  >  NC_002947/4988222‑4988506

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