Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,351 C→G 30.0% T619T (ACC→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3510CG30.0% 121.4 / 76.1 90T619T (ACC→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (14/13);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGC  >  NC_002947/196204‑196496
                                                                                                                                                       |                                                                                                                                                 
gACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCa                                                                                                                                                  >  1:28393/1‑149 (MQ=39)
     cGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTAcc                                                                                                                                             <  5:44885/149‑1 (MQ=39)
      gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTAcc                                                                                                                                             >  8:125708/1‑148 (MQ=39)
      gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCg                                                                                                                                            >  3:268158/1‑149 (MQ=39)
       aTCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGt                                                                                                                                           >  5:37168/1‑149 (MQ=39)
        tCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                                                          >  1:101687/1‑149 (MQ=39)
           aCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGaa                                                                                                                                       <  8:55732/149‑1 (MQ=39)
            cGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCg                                                                                                                                            <  6:268563/143‑1 (MQ=39)
             gACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCCACCCAATACCGTTGAAGc                                                                                                                                     >  8:294069/1‑149 (MQ=255)
                    tCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaag                                                                                                                              >  3:330081/1‑148 (MQ=38)
                         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                                                                                                        <  5:32850/140‑1 (MQ=34)
                                                   cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                        >  6:32850/1‑140 (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                                                                                     <  4:265161/149‑1 (MQ=34)
                                                             aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                     >  7:227612/1‑149 (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                                   <  3:5425/143‑1 (MQ=39)
                                                                                                                     ggTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcggcg                                   >  4:5425/1‑143 (MQ=39)
                                                                                                                            cgT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAGTTCCAGACCCCGGCGAATGAc                                                                                       <  6:42867/84‑1 (MQ=255)
                                                                                                                            cgT‑‑‑GACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                                                                       >  5:42867/2‑85 (MQ=255)
                                                                                                                              t‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                         <  7:287185/145‑1 (MQ=37)
                                                                                                                              t‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                         >  8:287185/1‑145 (MQ=37)
                                                                                                                                   aCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACt                   >  6:21754/1‑149 (MQ=37)
                                                                                                                                   aCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACt                   <  7:149253/149‑1 (MQ=37)
                                                                                                                                            aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                         <  7:293291/134‑1 (MQ=37)
                                                                                                                                            aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                         >  8:293291/1‑134 (MQ=37)
                                                                                                                                            aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa          <  1:56552/149‑1 (MQ=37)
                                                                                                                                                 ccAAACCATTACCTTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACTGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa     <  5:98259/149‑1 (MQ=255)
                                                                                                                                                    aaCCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGc  >  3:124852/1‑149 (MQ=39)
                                                                                                                                                       |                                                                                                                                                 
GACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGT‑‑‑GACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGC  >  NC_002947/196204‑196496

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: