Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,362 G→C 8.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,3620GC 8.3% 181.2 / 6.9 60intergenic (+32/+62)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (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 = 3.35e-01

ACCACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCC  >  NC_002947/4988214‑4988508
                                                                                                                                                     |                                                                                                                                                   
accaccTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACctg                                                                                                                                                     >  7:386385/1‑149 (MQ=255)
  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                                                                                      >  7:88102/1‑149 (MQ=255)
                                                              tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTcc                                                                                      >  2:253120/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)
                                                                                                                                                     |                                                                                                                                                   
ACCACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCC  >  NC_002947/4988214‑4988508

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