Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,372 A→G 7.0% K626K (AAA→AAG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3720AG 7.0% 238.4 / 6.9 87K626K (AAA→AAGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base G (3/3);  total (40/47)
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.68e-01

ACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGA  >  NC_002947/196229‑196505
                                                                                                                                                   |                                                                                                                                     
accacTACCGCGACCCTGACGGCGAGCCCGTCCGTCACCGAAGGCGCAGTGATAACTTACACAGTGACCCTCAGCAATCCTGCCCAGACGCCGGTGACACT‑‑‑GACCCTGTCCAACGGCCAAACAAATACAGTTGAAGCACGCAAGAcc                                                                                                                                    >  6:19728/1‑147 (MQ=11)
 ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                                                                                                 <  7:125708/149‑1 (MQ=17)
 ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                                                                                                 <  8:176725/149‑1 (MQ=17)
  cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                >  3:167938/1‑149 (MQ=17)
  cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAgg                                                                                                                                <  2:101687/149‑1 (MQ=37)
  cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCAAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                >  5:124020/1‑149 (MQ=17)
   acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg                                                                                                                               <  8:143229/149‑1 (MQ=17)
    ccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAggg                                                                                                                               >  4:331117/3‑148 (MQ=2)
         gcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGt                                                                                                                         >  8:217905/1‑149 (MQ=12)
             cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGac                                                                                                                     >  4:153434/1‑148 (MQ=12)
              ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAtt                                                                                                                    >  5:159268/1‑149 (MQ=34)
                tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTc                                                                                                                  >  3:265161/1‑149 (MQ=34)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                 <  2:16362/148‑1 (MQ=34)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCa                                                                                                                <  4:268158/149‑1 (MQ=34)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTcc                                                                                                                 <  2:16935/148‑1 (MQ=255)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTc                                                                                                                  >  5:359657/1‑147 (MQ=255)
                    ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                              <  1:216864/149‑1 (MQ=34)
                      cGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAAc                                                                                                                                     <  5:363233/124‑1 (MQ=35)
                      cGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAAc                                                                                                                                     >  6:363233/1‑124 (MQ=35)
                       gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                           >  5:353304/1‑149 (MQ=34)
                          cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                 >  6:32850/1‑140 (MQ=34)
                          cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                 <  5:32850/140‑1 (MQ=34)
                           ccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                       <  7:20107/149‑1 (MQ=34)
                             gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGAcc                                                                                                            >  6:250501/1‑142 (MQ=255)
                             gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGAcc                                                                                                            <  5:250501/142‑1 (MQ=255)
                                 gTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGGAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGa                                                                                                 >  5:315193/1‑149 (MQ=11)
                                   cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACg                                                                                               >  3:300735/1‑149 (MQ=255)
                                    aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                              >  7:227612/1‑149 (MQ=34)
                                    aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                              <  4:265161/149‑1 (MQ=34)
                                    aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                                                                                              >  6:274722/1‑149 (MQ=255)
                                     ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                             <  4:5990/149‑1 (MQ=34)
                                     ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                             <  3:153434/149‑1 (MQ=255)
                                      cGAAGGCGGCGTTATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGG‑‑‑GACCCTGTCCAACGGACAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCt                                                                                            <  8:310613/149‑1 (MQ=11)
                                      cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                             >  6:98259/1‑148 (MQ=34)
                                      cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                             >  2:56552/1‑148 (MQ=34)
                                      cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                             >  1:74615/1‑148 (MQ=255)
                                      cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCt                                                                                            <  2:74615/149‑1 (MQ=255)
                                       gtaGGCTGCGTGATCACTTACACCGTGACCCTGAGCATCCCTGCCCAGACGCCGTTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTa                                                                                           <  2:292637/147‑1 (MQ=11)
                                          ggcggcGTGATCACTTACAGCGTGACCCTGAGCAATCCTGGCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGTCAATACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        <  5:19728/149‑1 (MQ=255)
                                          ggcggcGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        <  6:159268/149‑1 (MQ=34)
                                           gcggcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaac                                                                                       <  6:353304/149‑1 (MQ=34)
                                              gcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac                                                                                    <  7:294069/149‑1 (MQ=34)
                                               cgTGATCACTTACACCGTGACCCTGAGCAACCCTCCCCAGACGCCGGTGACCCCTCAGGACGACACCCACTCGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTacaac                                                                                       <  1:116469/149‑1 (MQ=2)
                                                 tGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGt                                                                                 <  8:227612/149‑1 (MQ=34)
                                                  gATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGt                                                                                 <  6:123972/148‑1 (MQ=34)
                                                   aTCACTTACACCGTGTCCCTGAGCAATCCTGCGCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                <  4:204359/148‑1 (MQ=21)
                                                    tCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                               <  3:110704/148‑1 (MQ=255)
                                                      aCTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                             <  6:37168/148‑1 (MQ=34)
                                                       cTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                           <  3:331117/149‑1 (MQ=255)
                                                        ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTg                                                                          >  2:350679/1‑149 (MQ=255)
                                                        ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTg                                                                          >  3:120530/1‑149 (MQ=255)
                                                          acacCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                         >  8:149253/1‑148 (MQ=34)
                                                                tGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTc                                                                  <  4:300735/149‑1 (MQ=11)
                                                                    ccTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTaca                                                                                         <  4:55142/122‑1 (MQ=32)
                                                                    ccTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTaca                                                                                         >  3:55142/1‑122 (MQ=32)
                                                                    ccTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCa                                                              <  2:28393/149‑1 (MQ=37)
                                                                    ccTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        >  5:42074/1‑123 (MQ=255)
                                                                    ccTGAGCAACCCTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        <  6:42074/123‑1 (MQ=255)
                                                                      tGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        <  3:326413/121‑1 (MQ=32)
                                                                      tGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                        >  4:326413/1‑121 (MQ=32)
                                                                      tGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCAt                                                             >  6:70854/1‑148 (MQ=37)
                                                                          cAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                             >  3:245993/1‑112 (MQ=32)
                                                                          cAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                             <  4:245993/112‑1 (MQ=32)
                                                                           aaTCCTGCCCAGACGCCGGTGACAGT‑‑‑ACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAAc                                                      >  3:284074/1‑149 (MQ=37)
                                                                             cccTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTaca                                                                                         <  1:204218/112‑1 (MQ=255)
                                                                             cccTGCCCAGACGCCGGTGACCGT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTaca                                                                                         >  2:204218/2‑113 (MQ=255)
                                                                                   ccAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcg                                               >  8:60750/1‑149 (MQ=38)
                                                                                            ggTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcggcg                                            >  4:5425/1‑143 (MQ=39)
                                                                                            ggTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcggcg                                            <  3:5425/143‑1 (MQ=39)
                                                                                                   cgT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAGTTCCAGACCCCGGCGAATGAc                                                                                                <  6:42867/84‑1 (MQ=255)
                                                                                                   cgT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                                                                                >  5:42867/2‑85 (MQ=255)
                                                                                                     t‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                  >  8:287185/1‑145 (MQ=37)
                                                                                                     t‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                  <  7:287185/145‑1 (MQ=37)
                                                                                                          aCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACt                            >  6:21754/1‑149 (MQ=37)
                                                                                                          aCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACt                            <  7:149253/149‑1 (MQ=37)
                                                                                                                   aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                  >  8:293291/1‑134 (MQ=37)
                                                                                                                   aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                  <  7:293291/134‑1 (MQ=37)
                                                                                                                   aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                   <  1:56552/149‑1 (MQ=37)
                                                                                                                        ccAAACCATTACCTTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACTGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa              <  5:98259/149‑1 (MQ=255)
                                                                                                                           aaCCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGc           >  3:124852/1‑149 (MQ=39)
                                                                                                                                ttACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgtgt                                                                   >  3:59875/1‑88 (MQ=38)
                                                                                                                                ttACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgtgt                                                                   <  4:59875/88‑1 (MQ=38)
                                                                                                                                   ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa    <  5:21754/148‑1 (MQ=39)
                                                                                                                                    cGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGa  <  7:60750/149‑1 (MQ=39)
                                                                                                                                        gAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAAc                                                      >  2:277510/1‑93 (MQ=37)
                                                                                                                                        gAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAAc                                                      <  1:277510/93‑1 (MQ=37)
                                                                                                                                          aGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAccc                       >  4:355318/1‑122 (MQ=38)
                                                                                                                                                   |                                                                                                                                     
ACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGA  >  NC_002947/196229‑196505

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