Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,364 T→G 8.3% intergenic (+34/+60) 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,3640TG 8.3% 176.9 / 8.5 60intergenic (+34/+60)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base G (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 = 5.52e-01

CACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCA  >  NC_002947/4988216‑4988512
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                >  5:284339/1‑149 (MQ=255)
                                      aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGcc                                                                                                                >  3:397117/1‑149 (MQ=255)
                                          gcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCtt                                                                                                            <  4:397117/149‑1 (MQ=255)
                                             gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTcgc                                                                                                         <  8:127682/149‑1 (MQ=255)
                                             gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTcgc                                                                                                         <  6:209966/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                                                   >  8:68220/1‑147 (MQ=255)
                                                                                                     aGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctgc                                                   <  7:68220/147‑1 (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)
                                                                                                                                                     tgAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCa  >  1:35330/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCA  >  NC_002947/4988216‑4988512

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