Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 943,611 C→G 9.7% G671R (GGG→CGG) ‡ PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947943,6110CG 9.7% 173.1 / 9.0 62G671R (GGG→CGG) ‡PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (3/3);  total (33/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 = 4.92e-01

GACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAG  >  NC_002947/943471‑943758
                                                                                                                                            |                                                                                                                                                     
gACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGatat                                                                                                                                               >  3:95552/1‑149 (MQ=255)
    gTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCgg                                                                                                                                           >  3:4952/1‑149 (MQ=255)
           tCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACGGGGTCGatacc                                                                                                                                              >  5:38852/1‑137 (MQ=255)
           tCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACGGGGTCGatacc                                                                                                                                              <  6:38852/139‑3 (MQ=255)
             gcgAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAAc                                                                                                                                  >  5:206335/1‑149 (MQ=255)
               gAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACGGGGTCGatacca                                                                                                                                             <  5:136277/136‑4 (MQ=255)
               gAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACGGGGTCGatacca                                                                                                                                             >  6:136277/1‑133 (MQ=255)
                  gATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTgcgc                                                                                                                             >  4:198468/1‑149 (MQ=255)
                          ctgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGc                                                                                                                     <  5:1713/149‑1 (MQ=255)
                           tgccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGc                                                                                                                     >  5:5416/1‑148 (MQ=255)
                             gactgTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGa                                                                                                                   <  2:67576/144‑1 (MQ=255)
                             ccGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGAt                                                                                                                  >  1:185371/1‑149 (MQ=255)
                                 gTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATgtgt                                                                                                              >  3:26937/1‑149 (MQ=255)
                                               cGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGACGGGGCTGCCGGTGGCGTGGTATCGACCACGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCg                                                                                                >  7:81797/1‑149 (MQ=255)
                                                   agagCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCg                                                                                            >  6:270201/1‑149 (MQ=255)
                                                       ccGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGt                                                                                        >  6:262305/1‑149 (MQ=255)
                                                        cGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTg                                                                                       <  2:185371/149‑1 (MQ=255)
                                                        cGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTg                                                                                       <  6:206335/149‑1 (MQ=255)
                                                            ttCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGtt                                                                                                             >  5:281681/1‑123 (MQ=255)
                                                            ttCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGtt                                                                                                             <  6:281681/123‑1 (MQ=255)
                                                              cGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGGTCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCtt                                                                                   <  8:193319/147‑1 (MQ=255)
                                                                ccTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTgcag                                                                               <  4:95552/149‑1 (MQ=255)
                                                                  tCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTgcagca                                                                             >  7:182655/1‑149 (MQ=255)
                                                                    gACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACt                                                                           >  1:32701/1‑149 (MQ=255)
                                                                    gACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACt                                                                           >  1:37274/1‑149 (MQ=255)
                                                                             gTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTc                                                                   >  2:58069/1‑148 (MQ=255)
                                                                              ttGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCaa                                                                 >  8:84982/1‑149 (MQ=255)
                                                                                       ggCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑g                                                      >  6:157618/1‑149 (MQ=255)
                                                                                        gCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGCTGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑gg                                                     >  8:119232/1‑149 (MQ=255)
                                                                                        gCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑gg                                                     >  1:222874/1‑149 (MQ=255)
                                                                                         cGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑ggg                                                    <  2:32701/149‑1 (MQ=255)
                                                                                         cGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑ggg                                                    <  2:37274/149‑1 (MQ=255)
                                                                                           aGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTCTCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑Gggtt                                                  <  5:262305/149‑1 (MQ=255)
                                                                                           aGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑Gggtt                                                  <  3:198468/149‑1 (MQ=255)
                                                                                                gATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑Gggttgaggt                                             >  1:181907/1‑149 (MQ=255)
                                                                                                  gcACCAAGGCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCGGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑Gggttgaggttg                                           <  4:324/148‑1 (MQ=255)
                                                                                                    aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCg                                                                                                      <  4:18258/90‑1 (MQ=255)
                                                                                                    aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCg                                                                                                      >  3:18258/1‑90 (MQ=255)
                                                                                                    aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGgttgaggttgag                                         <  6:271644/149‑1 (MQ=255)
                                                                                                    aCCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGCCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑Gggttgaggttga                                          <  1:58069/148‑1 (MQ=255)
                                                                                                     ccAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGCTGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGttgaggttgagg                                        <  7:119232/149‑1 (MQ=255)
                                                                                                     ccAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGttgaggttgagg                                        >  2:177039/1‑149 (MQ=255)
                                                                                                     ccAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGttgaggttgagg                                        <  4:26937/149‑1 (MQ=255)
                                                                                                     ccAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGttgaggttgagg                                        >  4:144428/1‑149 (MQ=255)
                                                                                                       aaGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCt                                      <  8:182655/149‑1 (MQ=255)
                                                                                                        agtcagtcGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCt                                      >  4:163836/1‑148 (MQ=255)
                                                                                                         gtcagtcGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTg                                     >  3:268755/1‑148 (MQ=255)
                                                                                                            agtcGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAaca                                 <  1:177039/149‑1 (MQ=255)
                                                                                                                      gCCGGTGGCGTGGTATCGACACCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCgg                       >  6:270094/1‑149 (MQ=255)
                                                                                                                       ccGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCggg                      <  5:157618/149‑1 (MQ=255)
                                                                                                                       ccGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCgg                       <  2:222874/148‑1 (MQ=255)
                                                                                                                            ggCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCg                                                                                                      <  2:283094/66‑1 (MQ=255)
                                                                                                                            ggCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCg                                                                                                      >  1:283094/1‑66 (MQ=255)
                                                                                                                               gTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCAc              >  5:192015/1‑149 (MQ=255)
                                                                                                                                ttgTATCGACCGCGTCGATATCGTGTGCGGTAAGTGCGGTGCTTGGCGATGTGTTGCCAGCGGCTTCGGTCGCCGTGGCTTGCAGCACTTCGGCGTCATCTTGTGGG‑‑GGGGTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACg             <  8:81797/147‑1 (MQ=18)
                                                                                                                                 ggTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAac                                  <  3:99423/127‑1 (MQ=255)
                                                                                                                                 ggTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAac                                  >  4:99423/1‑127 (MQ=255)
                                                                                                                                 gatATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACgg            >  6:5519/3‑149 (MQ=255)
                                                                                                                                  gTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCggg                      <  6:133759/138‑1 (MQ=255)
                                                                                                                                  gTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCggg                      >  5:133759/1‑138 (MQ=255)
                                                                                                                                   tATCGACGGGGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCa          <  5:5519/149‑1 (MQ=255)
                                                                                                                                        aCCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGcc     <  4:4952/149‑1 (MQ=255)
                                                                                                                                           ccGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAg  >  5:206527/1‑149 (MQ=255)
                                                                                                                                            cGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGCATGGGTTAAGGTTGAGGATCAACACGCCGTACGGCCCCACCACGCAATTGCCCa   >  1:253032/1‑149 (MQ=18)
                                                                                                                                            |                                                                                                                                                     
GACGGTGCCGGTCGCGAGGATGTTGCCTGCCGCGTCACGTACCTGCACGGTAGAGCCGGCTTCGCCTCGACCACTGAGTTGGCTGCCGGCGAGGCCGATCACCAAGTCAGTCGGGGCTGCCGGTGGCGTGGTATCGACCCCGTCGATATCCGGTGCGGTAACTGCGCTGCTTGGCGATGTGTTGCCAGCGGCATCGGTCGCCGTGGCTTGCAGCACTTCGCCGTCAACTTGTGGC‑‑GGGTTGAGGTTGAGGCTGAACACGCCGGTCGGGCCGACCACGGCAGTGCCCAG  >  NC_002947/943471‑943758

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