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% 182.9 / 7.0 60intergenic (+33/+61)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/3);  total (29/31)
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.75e-01

CACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGC  >  NC_002947/4988216‑4988511
                                                                                                                                                    |                                                                                                                                                     
caccTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACctgct                                                                                                                                                      >  1:241956/1‑149 (MQ=255)
 accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACGAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACctgctg                                                                                                                                                     <  8:245456/149‑1 (MQ=255)
 accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACctgct                                                                                                                                                      >  5:209966/1‑148 (MQ=255)
    tttGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAgg                                                                                                                                                  >  2:379877/1‑149 (MQ=255)
         cGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTCAGCAGGTGGCg                                                                                                                                             <  5:409127/149‑1 (MQ=255)
          ggCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGt                                                                                                                                            <  2:241956/149‑1 (MQ=255)
           gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGtt                                                                                                                                           <  1:373595/149‑1 (MQ=255)
           gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGtt                                                                                                                                           <  1:64141/149‑1 (MQ=255)
                cgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCa                                                                                                                                       <  7:238256/148‑1 (MQ=255)
                 gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCAtt                                                                                                                                     >  1:271327/1‑149 (MQ=255)
                 gattCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCAtt                                                                                                                                     <  2:284688/147‑1 (MQ=255)
                  tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                    <  8:386385/149‑1 (MQ=255)
                      cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACATAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCAt                                                                                                                                      <  4:24573/143‑1 (MQ=255)
                      cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑CCCTGCTGAGGTGGCGTGTTCATTTgtgt                                                                                                                                >  4:249168/1‑149 (MQ=255)
                      cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCAt                                                                                                                                      >  3:24573/1‑143 (MQ=255)
                         aaCCTGAGCGGGCAGGCACAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTGTCATTTGTGTTTc                                                                                                                             >  8:134963/1‑149 (MQ=255)
                                gCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTgct‑ggg                                                                                                                     <  2:271327/149‑1 (MQ=255)
                                    gCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGGCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTgct‑gggctgg                                                                                                                 <  1:294407/149‑1 (MQ=255)
                                     cAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGc                                                                                                                <  1:379877/149‑1 (MQ=255)
                                      aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGcc                                                                                                               >  3:397117/1‑149 (MQ=255)
                                      aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGcc                                                                                                               >  5:284339/1‑149 (MQ=255)
                                          gcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCtt                                                                                                           <  4:397117/149‑1 (MQ=255)
                                             gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTcgc                                                                                                        <  6:209966/149‑1 (MQ=255)
                                             gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTcgc                                                                                                        <  8:127682/149‑1 (MQ=255)
                                                       gagaGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAAccc                                                                                              >  1:62688/1‑149 (MQ=255)
                                                           gTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCAGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCt                                                                                           <  3:109798/148‑1 (MQ=255)
                                                            tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTcc                                                                                         >  2:253120/1‑149 (MQ=255)
                                                            tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTcc                                                                                         >  7:88102/1‑149 (MQ=255)
                                                                     ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGa                                                                                 <  1:253120/148‑1 (MQ=255)
                                                                           gCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCt                                                                           >  4:112963/1‑148 (MQ=255)
                                                                            cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTg                                                                          <  6:284339/148‑1 (MQ=255)
                                                                                   aGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTc                                                                  <  3:112963/149‑1 (MQ=255)
                                                                                     gcgcGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAg                                                                <  8:88102/149‑1 (MQ=255)
                                                                                       gcggcgGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgc                                                              >  2:218232/1‑149 (MQ=255)
                                                                                        cggcggTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgcg                                                             >  6:268553/1‑149 (MQ=255)
                                                                                                  aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGTGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                   >  5:256938/1‑149 (MQ=255)
                                                                                                  aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAACGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                   >  6:362376/1‑149 (MQ=255)
                                                                                                  aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCCGTTTGATTTGTGATTCTGCT‑GGGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGACTTGCTGACCTCTTCAGGCGCCTGTtgct                                                    >  4:379054/1‑148 (MQ=39)
                                                                                                    aaGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCGACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                                  >  1:285647/1‑149 (MQ=255)
                                                                                                    aaGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcg                                                 <  2:285647/149‑1 (MQ=255)
                                                                                                     aGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                                  <  7:68220/147‑1 (MQ=255)
                                                                                                     aGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                                  >  8:68220/1‑147 (MQ=255)
                                                                                                       gACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                              >  2:315747/1‑149 (MQ=255)
                                                                                                          gATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCt                                           >  6:74519/1‑149 (MQ=255)
                                                                                                               ttGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGTGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                      <  6:256938/149‑1 (MQ=255)
                                                                                                               ttGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTctct                                      <  1:218232/149‑1 (MQ=255)
                                                                                                                  aTCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTc                                    <  5:74519/148‑1 (MQ=255)
                                                                                                                      tcggtcggTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCt                               >  7:128643/1‑149 (MQ=255)
                                                                                                                        ggtcggtAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                             <  5:268553/149‑1 (MQ=255)
                                                                                                                         gtcggtAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                             >  1:117045/1‑149 (MQ=255)
                                                                                                                          tcggtAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgt                           >  5:57925/1‑149 (MQ=255)
                                                                                                                          tcggtAAAACAGAAAACGCC‑ACCTGCTGAGGGGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgt                           <  8:147674/149‑1 (MQ=255)
                                                                                                                                  aCAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAg                   >  8:332368/1‑149 (MQ=255)
                                                                                                                                  aCAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCAGCT‑GCGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTGTTCAg                   >  2:242182/1‑149 (MQ=255)
                                                                                                                                   cAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGccc                                            <  1:199099/123‑1 (MQ=255)
                                                                                                                                    aGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTTCACCGACTTGCTGACCTCTTCAGGCGCCTGTTGGTGCGCGCCCTTCTCTTCATACTCGGTGTCTTCAgcc                 >  7:146817/1‑149 (MQ=255)
                                                                                                                                    aGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgcc                 <  3:9504/149‑1 (MQ=255)
                                                                                                                                     gAAAACGCC‑ACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGccc                                            >  3:154929/1‑121 (MQ=255)
                                                                                                                                     gAAAACGCC‑ACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGccc                                            <  4:154929/121‑1 (MQ=255)
                                                                                                                                      aaaaCGCC‑ACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTc                                    >  2:224198/1‑128 (MQ=255)
                                                                                                                                      aaaaCGCC‑ACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTc                                    <  1:224198/128‑1 (MQ=255)
                                                                                                                                               aCCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGgcg       >  6:221010/1‑149 (MQ=255)
                                                                                                                                                ccTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGgcgc      <  2:62688/149‑1 (MQ=255)
                                                                                                                                                 ctgctgAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGcc     >  6:223850/1‑149 (MQ=255)
                                                                                                                                                    ctgAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCt    <  5:221010/147‑1 (MQ=255)
                                                                                                                                                    ctgAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGc  <  6:57925/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGC  >  NC_002947/4988216‑4988511

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