Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,365 G→C 8.3% 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.3% 186.5 / 7.0 60intergenic (+35/+59)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 = 8.24e-01

ACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCA  >  NC_002947/4988217‑4988512
                                                                                                                                                     |                                                                                                                                                    
accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACGAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACctgctg                                                                                                                                                      <  8:245456/149‑1 (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:64141/149‑1 (MQ=255)
          gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGtt                                                                                                                                            <  1:373595/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                                                                                                         <  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)
                                                                                                                                                     |                                                                                                                                                    
ACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCC‑ACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCT‑GGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCA  >  NC_002947/4988217‑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: ‑