Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,468 C→G 7.2% R153P (CGC→CCC) ‡ PP_4671 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,4680CG 7.2% 268.1 / 8.6 83R153P (CGC→CCC) ‡PP_4671hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/38);  new base G (2/4);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

GCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCA  >  NC_002947/5297321‑5297602
                                                                                                                                                   |                                                                                                                                      
gcgcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTcgcg                                                                                                                                       <  8:159827/149‑1 (MQ=255)
  gcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGa                                                                                                                                     >  6:15306/1‑149 (MQ=255)
  gcAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGa                                                                                                                                     <  4:55407/149‑1 (MQ=255)
       gTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCg                                                                                                                                 >  5:252005/1‑148 (MQ=255)
          gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcagg                                                                                                                             >  2:223808/1‑149 (MQ=255)
          gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  1:71545/149‑1 (MQ=255)
            gcgcACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggca                                                                                                                           <  5:124752/149‑1 (MQ=255)
                aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAAc                                                                                                                       <  7:90839/149‑1 (MQ=255)
                aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcaggtgacct                                                                                                                       <  8:131444/149‑7 (MQ=255)
                      cGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTg                                                                                                                 <  3:98248/149‑1 (MQ=255)
                       gCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGc                                                                                                                >  6:52628/1‑149 (MQ=255)
                         ttGCGGGTC‑GCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGa                                                                                                             >  5:40306/1‑149 (MQ=255)
                            cGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGAcc                                                                                                           <  6:369271/149‑1 (MQ=255)
                              ggTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGt                                                                                                         <  2:340828/149‑1 (MQ=255)
                                tCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTg                                                                                                        >  6:129849/1‑148 (MQ=255)
                                    cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                    >  5:43950/1‑132 (MQ=255)
                                    cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                    <  6:43950/132‑1 (MQ=255)
                                    cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  6:40306/149‑1 (MQ=255)
                                    cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  6:252005/149‑1 (MQ=255)
                                    cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                   <  1:223808/149‑1 (MQ=255)
                                        aTGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCa                                                                                               >  5:62210/1‑149 (MQ=255)
                                              gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGc                                                                                          >  2:189237/1‑148 (MQ=255)
                                              gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  6:367316/113‑1 (MQ=255)
                                              gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             <  6:362078/113‑1 (MQ=255)
                                              gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             >  5:367316/1‑113 (MQ=255)
                                              gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                             >  5:362078/1‑113 (MQ=255)
                                                tCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCAGcgag                                                                                       >  5:76276/1‑146 (MQ=255)
                                                  cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     <  4:134900/149‑1 (MQ=255)
                                                  cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     >  7:51874/1‑149 (MQ=255)
                                                  cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                     <  2:267547/149‑1 (MQ=255)
                                                       ggCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATg                                                                                <  5:52628/149‑1 (MQ=255)
                                                          gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCt                                                                              >  3:276124/1‑148 (MQ=255)
                                                          gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             <  6:171720/149‑1 (MQ=255)
                                                           ccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCa                                                                            <  4:353808/149‑1 (MQ=255)
                                                               gACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCAc                                                                        >  5:74253/1‑149 (MQ=255)
                                                                    ttGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGt                                                                   <  8:13366/149‑1 (MQ=255)
                                                                        agcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgcc                                                               <  6:320672/149‑1 (MQ=255)
                                                                         gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             <  4:60449/134‑1 (MQ=255)
                                                                         gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                             >  3:60449/1‑134 (MQ=255)
                                                                         gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccg                                                              >  4:151736/1‑149 (MQ=255)
                                                                          caCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCACAGCACCGGGTGgccgc                                                             >  5:305333/1‑149 (MQ=255)
                                                                           aCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccgcc                                                            <  8:51874/149‑1 (MQ=255)
                                                                             gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          <  6:208112/149‑1 (MQ=255)
                                                                             gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          <  6:213236/149‑1 (MQ=255)
                                                                             gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                          >  8:98927/1‑149 (MQ=255)
                                                                               tGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATa                                                        <  6:74253/149‑1 (MQ=255)
                                                                                 gggTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGATGATGCTCAGCACCGGGTGGCCGCCAATACg                                                      <  3:52133/149‑1 (MQ=255)
                                                                                   gTGAATGCGCACGCCAGAGCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGCCCGGGCATGGGTCATCCGCCTAGCTGATGCTCAGCACCGGGTGGCCTCCACTACGCt                                                    >  8:79097/1‑149 (MQ=255)
                                                                                   gTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCt                                                    >  4:234366/1‑149 (MQ=255)
                                                                                    tGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTg                                                   <  3:151736/149‑1 (MQ=255)
                                                                                           gcACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGc                                            >  7:124427/1‑149 (MQ=255)
                                                                                            cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                           <  3:208747/149‑1 (MQ=255)
                                                                                            cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                           <  2:313043/149‑1 (MQ=255)
                                                                                              cGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                         <  5:15306/149‑1 (MQ=255)
                                                                                              cGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                         <  6:62210/149‑1 (MQ=255)
                                                                                                 cAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcg                                      >  6:333990/1‑149 (MQ=255)
                                                                                                  agagCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGCCCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgc                                     >  4:156291/1‑149 (MQ=255)
                                                                                                   gagCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgcg                                    >  1:163117/1‑149 (MQ=255)
                                                                                                      cGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCAGCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTg                                 <  7:79097/149‑1 (MQ=255)
                                                                                                        tGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGc                               >  4:131786/1‑149 (MQ=255)
                                                                                                          ccACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCAACCGTGCATGGGCCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGccc                             >  4:90506/1‑149 (MQ=255)
                                                                                                           cACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCt                            >  7:285108/1‑149 (MQ=255)
                                                                                                                      gaTTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCTGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTg                 <  7:314789/149‑1 (MQ=255)
                                                                                                                          gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                            >  2:236687/1‑38 (MQ=255)
                                                                                                                          gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                            <  1:236687/38‑1 (MQ=255)
                                                                                                                          gACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCaa             >  8:82612/1‑149 (MQ=255)
                                                                                                                            cTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg           >  1:120344/1‑149 (MQ=255)
                                                                                                                             ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCTTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg           >  7:171400/1‑148 (MQ=255)
                                                                                                                             ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg           >  5:33470/1‑148 (MQ=255)
                                                                                                                             ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGAAGAACAGGCTGCAGACCGTGCATGGGTCATACGCCTCGCTGATGATCAGCACCGTGTGGCCGCCAATACGCTGGAACAGCCGATAGCGCTGGCCCTCGTCAAGCTTGCCAAtc           >  4:238422/1‑147 (MQ=255)
                                                                                                                              tCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg           >  3:242810/1‑147 (MQ=255)
                                                                                                                               cAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGt          <  3:238422/147‑1 (MQ=255)
                                                                                                                                aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTaga       <  2:120344/149‑1 (MQ=255)
                                                                                                                                aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTa         >  5:275699/1‑147 (MQ=255)
                                                                                                                                 aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag      <  3:353608/149‑1 (MQ=255)
                                                                                                                                 aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag      <  6:33470/149‑1 (MQ=255)
                                                                                                                                    gTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAggg    >  7:248779/1‑148 (MQ=255)
                                                                                                                                     tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCTCTGGTCAAGCTTGCCAATGTAGAGGGCa  >  1:359187/1‑149 (MQ=255)
                                                                                                                                     tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGc   >  4:81272/1‑148 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                    <  8:150774/65‑1 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                    >  7:150774/1‑65 (MQ=255)
                                                                                                                                          tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                <  6:120290/114‑1 (MQ=255)
                                                                                                                                          tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                >  5:120290/1‑114 (MQ=255)
                                                                                                                                            ctgGTCGCGGACCCGCAGGCAGAACAGGCTTCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGTCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGcc               <  3:367777/129‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                      
GCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCA  >  NC_002947/5297321‑5297602

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