Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,365 G→C 8.0% intergenic (+35/+59) 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,3650GC 8.0% 188.1 / 4.3 60intergenic (+35/+59)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (23/32);  new base C (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 = 1.00e+00

TGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCT  >  NC_002947/4988222‑4988509
                                                                                                                                               |                                                                                                                                                
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)
                                                                                                                                           accggAGGTGGCTGTTTCATTTGTGTTTCTGCTGTGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCTCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCt  <  8:125185/145‑1 (MQ=39)
                                                                                                                                               |                                                                                                                                                
TGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCT  >  NC_002947/4988222‑4988509

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: ‑