Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 887,770 G→A 6.4% R268R (CGG→CGA)  PP_0768 → TPR domain‑containing response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947887,7700GA 6.4% 229.7 / 3.1 77R268R (CGG→CGA) PP_0768TPR domain‑containing response regulator
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base A (3/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

GTGTTTACGAGCAGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAG  >  NC_002947/887622‑887913
                                                                                                                                                    |                                                                                                                                               
gtgtTTACGAGCAGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCgg                                                                                                                                                 <  1:41773/149‑1 (MQ=255)
   tttACGAGCAGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCa                                                                                                                                               <  5:190003/148‑1 (MQ=255)
       cGAGCAGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGa                                                                                                                                          >  6:47927/1‑149 (MQ=255)
           cAGGCGCTCCAGGCGTTGCCGATGATGCCAGGGGTTTAGGATGGCATGGCCGAGGTGATGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTGGATGAAGCCGTGCGCCTGTCGCGTTTGTCGGTGCGTCGGCAGTCGACGCt                                                                                                                                      >  4:151994/1‑149 (MQ=255)
            aGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTg                                                                                                                                     >  4:277139/1‑149 (MQ=255)
             ggCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTgg                                                                                                                                    <  1:269986/149‑1 (MQ=255)
             ggCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGACTGCCGACGcacc                                                                                                                                    <  5:273314/149‑4 (MQ=255)
               cgcTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGc                                                                                                                                  <  3:179569/149‑1 (MQ=255)
                    aaGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTc                                                                                                                                            >  3:125485/1‑134 (MQ=255)
                    aaGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTc                                                                                                                                            <  4:125485/134‑1 (MQ=255)
                     aGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCaa                                                                                                                                <  3:267571/145‑1 (MQ=255)
                     aGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCaa                                                                                                                                >  4:267571/1‑145 (MQ=255)
                     aGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTg                                                                                                                            <  2:275121/149‑1 (MQ=255)
                     aGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTg                                                                                                                            <  4:227816/149‑1 (MQ=255)
                           tCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGTGGGATACACGGCGCGCGCAGAACCTGCTTGAATAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGAAGCTGTGCAAGCTGGCCCTg                                                                                                                      >  7:289236/1‑149 (MQ=255)
                            cccGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTg                                                                                                                      >  5:20165/1‑148 (MQ=255)
                            cccGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTg                                                                                                                      <  7:112737/148‑1 (MQ=255)
                             ccGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGga                                                                                                                    <  6:20165/149‑1 (MQ=255)
                              cGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTg                                                                                                                      >  3:225025/1‑146 (MQ=255)
                                aTCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGaa                                                                                                                 <  2:66224/149‑1 (MQ=255)
                                    tGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACg                                                                                                               >  3:314362/1‑147 (MQ=255)
                                    tGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGaa                                                                                                             <  3:277139/149‑1 (MQ=255)
                                       cAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGAc                                                                                                          <  7:109199/149‑1 (MQ=255)
                                           ccTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCg                                                                                                      <  2:273221/149‑1 (MQ=255)
                                            cTTTACTCTGTCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCTCGCGCAGAACCTGCTTGAAGAAGCGGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCga                                                                                                     <  8:289236/149‑1 (MQ=255)
                                                  gATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGaa                                                                                                             >  7:69098/1‑135 (MQ=255)
                                                  gATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGaa                                                                                                             <  8:69098/135‑1 (MQ=255)
                                                  gATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAgcg                                                                                                >  5:158895/1‑148 (MQ=255)
                                                   atggcatggcCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCa                                                                                              <  4:256405/149‑1 (MQ=255)
                                                      gcatggcCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATc                                                                                            <  4:284634/148‑1 (MQ=255)
                                                        atggcCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCg                                                                                           >  5:110304/1‑147 (MQ=255)
                                                         tggcCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAg                                                                                        <  6:110304/149‑1 (MQ=255)
                                                                     tggtggCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCAt                                                                            >  6:25547/1‑149 (MQ=255)
                                                                          gCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAg                                                                                        >  7:46538/1‑132 (MQ=255)
                                                                          gCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAg                                                                                        <  8:46538/132‑1 (MQ=255)
                                                                           cTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTg                                                                      <  8:8383/149‑1 (MQ=255)
                                                                            tCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGaa                                                                                                             <  5:108871/109‑1 (MQ=255)
                                                                            tCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGaa                                                                                                             >  6:108871/1‑109 (MQ=255)
                                                                            tCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTg                                                                      >  5:329622/1‑148 (MQ=255)
                                                                              aGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAAc                                                                   <  6:124162/149‑1 (MQ=255)
                                                                              aGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAAc                                                                   <  5:47927/149‑1 (MQ=255)
                                                                               gggggATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGTCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAAcc                                                                  <  1:46718/149‑1 (MQ=255)
                                                                                 gggATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAg                                                                <  4:225025/149‑1 (MQ=255)
                                                                                  ggATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCTAAGGCGTTCCGTCATGCAGTGAACCAgg                                                               >  5:131055/1‑149 (MQ=255)
                                                                                   gATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCCGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAggg                                                              >  5:231113/1‑149 (MQ=255)
                                                                                       cccGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAgggg                                                             >  4:239057/1‑146 (MQ=255)
                                                                                        ccGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCaga                                                         <  3:192547/149‑1 (MQ=255)
                                                                                         cGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCagag                                                        >  7:89320/1‑149 (MQ=255)
                                                                                           gcgcgcgcAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGtt                                                      <  6:158895/149‑1 (MQ=255)
                                                                                                          gCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATg                                       >  7:251784/1‑149 (MQ=255)
                                                                                                            ttGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGt                                                       >  4:37357/1‑131 (MQ=255)
                                                                                                            ttGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGt                                                       <  3:37357/131‑1 (MQ=255)
                                                                                                              gaagaaGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCag                                                          <  8:170202/126‑1 (MQ=255)
                                                                                                              gaagaaGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCag                                                          >  7:170202/1‑126 (MQ=255)
                                                                                                               aagaagCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCga                                                                                                     >  5:87427/1‑82 (MQ=255)
                                                                                                               aagaagCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCga                                                                                                     <  6:87427/82‑1 (MQ=255)
                                                                                                                agaagCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGaaa                                 <  5:25547/149‑1 (MQ=255)
                                                                                                                 gaagCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGaaaa                                >  3:116268/1‑149 (MQ=255)
                                                                                                                   agCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAaaca                              <  3:151994/149‑1 (MQ=255)
                                                                                                                     ccGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGtt                                                      >  5:123632/1‑123 (MQ=255)
                                                                                                                     ccGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGtt                                                      <  6:123632/123‑1 (MQ=255)
                                                                                                                     ccGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAaacaac                            >  5:9690/1‑149 (MQ=255)
                                                                                                                      cGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGTGCAAGCTGGAGCTGGAGAACGAAGACTTCGATAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGATTTCGCGCTACAAGGATGCCGAAaacaac                            <  1:297630/148‑1 (MQ=255)
                                                                                                                              ggTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGTCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTg                   <  6:231113/148‑1 (MQ=255)
                                                                                                                                tCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTg                   <  8:89320/147‑1 (MQ=255)
                                                                                                                                  gCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAaacaa                             >  6:327140/1‑135 (MQ=255)
                                                                                                                                  gCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTg                >  4:107851/1‑148 (MQ=255)
                                                                                                                                        gTCGGTGCGTCGGCAGTCGACGCTGGGCGAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGt                 >  2:81644/1‑141 (MQ=255)
                                                                                                                                        gTCGGTGCGTCGGCAGTCGACGCTGGGCGAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGt                 <  1:81644/141‑1 (MQ=255)
                                                                                                                                        gTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTa                                               >  7:47629/1‑111 (MQ=255)
                                                                                                                                        gTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTa                                               <  8:47629/111‑1 (MQ=255)
                                                                                                                                          cccagCGTCGACTGCCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCCTGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGtt        >  5:235576/5‑148 (MQ=255)
                                                                                                                                           ccagCGTCGACTGCCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGaaaa                                >  6:4230/4‑123 (MQ=255)
                                                                                                                                           ccagCGTCGACTGCCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGaaaa                                <  5:4230/120‑1 (MQ=255)
                                                                                                                                            gTGGGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTtgat     <  6:9690/149‑1 (MQ=255)
                                                                                                                                            gTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCAGCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTtgat     <  6:329622/149‑1 (MQ=255)
                                                                                                                                            cagCGTCGACTGCCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAAACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTAGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTtgat     >  4:172945/3‑149 (MQ=255)
                                                                                                                                                gTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAg  >  3:135184/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTGTTTACGAGCAGGCGCTCAAGGCGTTCCCGATCATGCCAGGCCTTTACGATGGCATGGCCGAGGTGCTGGTGGCTCAGGGGGATACCCGGCGCGCGCAGAACCTGCTTGAAGAAGCCGTGCGCCTGTCGCCTTTGTCGGTGCGTCGGCAGTCGACGCTGGGCAAGCTGGCGCTGGAGAACGAAGACTTCGAGAGCGCATCGAAGGCGTTCCGTCATGCAGTGAACCAGGGGCAGAGTTCGCGCTACAAGGATGCCGAAAACAACCTTGGCCTGGTGCAGGCGTTGATGAG  >  NC_002947/887622‑887913

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