Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,435,820 C→T 7.1% A244A (GCC→GCT gpJ → baseplate/tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,435,8200CT 7.1% 163.3 / 2.4 56A244A (GCC→GCTgpJbaseplate/tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base T (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GCCATCCTGCCCAGTCTGTCCAGCACCACACCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTC  >  NC_002947/3435674‑3435968
                                                                                                                                                        |                                                                                                                                                     
gCCATCCTGCCCAGTCTGTCCAGCACCACACCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCGGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGcc                                                                                                                                                       <  8:89805/147‑1 (MQ=255)
            aGTCTGTCCAGCACCACACCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGa                                                                                                                                         <  5:233568/149‑1 (MQ=255)
                tGTCCAGCACCACACCCGTGGCCCAGCTGGTGGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCTGAGCTGCAGATCGAgg                                                                                                                                       >  1:43334/1‑148 (MQ=255)
                tGTCCAGCACCACACCCGTGGCCCAGCTGGTGGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCTGAGCTGCAGATCGAgg                                                                                                                                       >  1:48596/1‑148 (MQ=255)
                  tCCAGCACCACACCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAg                                                                                                                                    <  8:122457/148‑1 (MQ=255)
                          cacaCCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCAGCAGATCGAGGAAGGGCCAAAcc                                                                                                                           >  1:335367/1‑149 (MQ=255)
                            cacCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCt                                                                                                                         <  3:336341/149‑1 (MQ=255)
                                  tGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGTAGATCGAGGAAGGGCCAGACCCTGACgtg                                                                                                                   <  1:27724/149‑1 (MQ=255)
                                      ccAGCTGGTGGCCAAGGTCACGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGa                                                                                                               <  6:78916/149‑1 (MQ=255)
                                      ccAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtggtg                                                                                                                >  8:222137/1‑148 (MQ=255)
                                      ccAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtggtg                                                                                                                >  5:297645/1‑148 (MQ=255)
                                       cAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCCGGCCCCGCGTGGCGTTTCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtggtg                                                                                                                <  1:292647/149‑1 (MQ=255)
                                              tggCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACg                                                                                                       <  7:187207/149‑1 (MQ=255)
                                                  ccAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGgcg                                                                                                    >  5:289665/1‑149 (MQ=255)
                                                   cAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAgg                                                                                                                                   <  8:202637/117‑1 (MQ=255)
                                                   cAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAgg                                                                                                                                   >  7:202637/1‑117 (MQ=255)
                                                   cAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGgcgc                                                                                                   >  6:149931/1‑149 (MQ=255)
                                                   cAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGgcgc                                                                                                   >  5:299333/1‑149 (MQ=255)
                                                     aGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAg                                                                                                                                                    >  6:302966/1‑98 (MQ=255)
                                                     aGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAg                                                                                                                                                    <  5:302966/98‑1 (MQ=255)
                                                      ttcctcggCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtg                                                                                                                   >  1:149348/7‑130 (MQ=255)
                                                      ttcctcggCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtg                                                                                                                   <  2:149348/124‑1 (MQ=255)
                                                       gTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCggg                                                                                                <  2:335367/148‑1 (MQ=255)
                                                        tCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtg                                                                                                                   >  1:348588/1‑128 (MQ=255)
                                                        tCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACgtg                                                                                                                   <  2:348588/128‑1 (MQ=255)
                                                            ggccggccGGTCGGCCGATGTATTACCACCAAACCCCGCCCCGGCCCGCGGGGCCGGCGGGGTACCGCCCGTGTTCAATCCCATGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTg                                                                                             <  1:287862/149‑1 (MQ=255)
                                                            ggCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCGA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGcc                                                                                          >  7:195416/1‑149 (MQ=255)
                                                                 cccTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCTGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAAt                                                                                     <  2:43334/149‑1 (MQ=255)
                                                                 cccTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCTGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAAt                                                                                     <  2:48596/149‑1 (MQ=255)
                                                                 cccTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAAt                                                                                     <  6:297645/149‑1 (MQ=255)
                                                                        ggCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCa                                                                              >  5:155268/1‑149 (MQ=255)
                                                                              ggAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTg                                                                           >  3:121315/1‑146 (MQ=255)
                                                                               gAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGCCCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAAt                                                                                     >  5:362716/1‑135 (MQ=255)
                                                                               gAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAAt                                                                                     <  6:362716/135‑1 (MQ=255)
                                                                                  gTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGc                                                                    >  6:145337/1‑149 (MQ=255)
                                                                                        ccAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAAATCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCg                                                              >  1:255863/1‑149 (MQ=255)
                                                                                              accgaccgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGcc                                                                                                                                >  2:104570/1‑77 (MQ=255)
                                                                                              accgaccgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGcc                                                                                                                                <  1:104570/77‑1 (MQ=255)
                                                                                              accgaccgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGc                                                                          <  4:186431/131‑1 (MQ=255)
                                                                                              accgaccgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGc                                                                          >  3:186431/1‑131 (MQ=255)
                                                                                                cgaccgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAGGCCCGTCGCCTGCAgg                                                      <  4:121315/149‑1 (MQ=255)
                                                                                                  accgGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAgggg                                                    <  7:17783/149‑1 (MQ=255)
                                                                                                          ccGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGGCGCCTGCAGGGGCAGTTGc                                             <  6:289665/149‑1 (MQ=255)
                                                                                                            gcgTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCgg                                                                                 <  3:48483/111‑1 (MQ=255)
                                                                                                            gcgTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCgg                                                                                 >  4:48483/1‑111 (MQ=255)
                                                                                                                ggCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGt                                       >  2:215033/1‑149 (MQ=255)
                                                                                                                      ccgGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGcc                                 <  6:155268/149‑1 (MQ=255)
                                                                                                                             aCCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGc                                                                                           <  7:83438/85‑1 (MQ=255)
                                                                                                                             aCCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGc                                                                                           >  8:83438/1‑85 (MQ=255)
                                                                                                                                 aCTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCa                                                                         <  7:304024/99‑1 (MQ=255)
                                                                                                                                 aCTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCa                                                                         >  8:304024/1‑99 (MQ=255)
                                                                                                                                 aCTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCAt                               <  2:73431/140‑1 (MQ=255)
                                                                                                                                 aCTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCAt                               >  1:73431/1‑140 (MQ=255)
                                                                                                                                         cAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTacgc               <  7:172856/148‑3 (MQ=255)
                                                                                                                                         cAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACg              >  2:9283/1‑149 (MQ=255)
                                                                                                                                          aa‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACg              >  7:147253/1‑148 (MQ=255)
                                                                                                                                                    aGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCCGGCCATCTATGCCGCGTTACACGTAACggg       <  1:215033/149‑1 (MQ=255)
                                                                                                                                                       ccGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGAtt   <  6:299333/149‑1 (MQ=255)
                                                                                                                                                        cGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTc  <  8:147253/149‑1 (MQ=255)
                                                                                                                                                        cGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCAGCGTTACACGTAACGGGGATTc  >  7:70029/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GCCATCCTGCCCAGTCTGTCCAGCACCACACCCGTGGCCCAGCTGGTGGCCAAGGTCAAGGCCGCCCTGTCGGCCGAGGAAGTACGACCAATCACCGACCGGCCCCGCGTGGCCGCCGGCGTACCGACTGAGTTCAA‑‑‑GGTGGAGGCCGAGCTGCAGATCGAGGAAGGGCCAGACCCTGACGTGGTGAAGGCCACGGCGCGGGCTGGCCTGAATGCGGCCATTGCAGAAGCCCGTCGCCTGCAGGGGCAGTTGCCCCTGTCGGCCATCTATGCCGCGTTACACGTAACGGGGATTC  >  NC_002947/3435674‑3435968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: