Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 932,607 T→C 8.6% N4339D (AAC→GAC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947932,6070TC 8.6% 159.8 / 5.5 58N4339D (AAC→GAC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base C (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

ACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGTGCTGCCGGC  >  NC_002947/932459‑932754
                                                                                                                                                    |                                                                                                                                                   
aCGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGtt                                                                                                                                                     <  8:201920/149‑1 (MQ=255)
 cGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTg                                                                                                                                                    <  6:206899/149‑1 (MQ=255)
   gTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGcc                                                                                                                                                  <  6:55953/149‑1 (MQ=255)
    tACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCg                                                                                                                                                 >  8:154812/1‑149 (MQ=255)
    tACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCg                                                                                                                                                 <  7:2794/149‑1 (MQ=255)
        gTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCAGCCGGCAACGTGTCGCCGGccg                                                                                                                                             <  6:59340/149‑3 (MQ=255)
        gTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCAGCCGGCAACGTGTCGCCGGccg                                                                                                                                             <  8:170020/149‑3 (MQ=255)
           cGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCAGCCGGCAACGTGTCGCCGGCCGca                                                                                                                                           <  6:306737/148‑2 (MQ=255)
             gTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCg                                                                                                                                        >  3:242441/1‑149 (MQ=255)
             gTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGGGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCg                                                                                                                                        >  3:48120/1‑149 (MQ=255)
                   gTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTgg                                                                                                                                    <  7:154812/147‑1 (MQ=255)
                   gTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCg                                                                                                                                  >  3:295368/1‑149 (MQ=255)
                   gTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCg                                                                                                                                  <  3:285150/149‑1 (MQ=255)
                           aCGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCt                                                                                                                           <  6:138540/148‑1 (MQ=255)
                                   gcCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGGATCt                                                                                                                  >  4:45063/1‑149 (MQ=255)
                                   gcCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCAGCTGCGGTGAAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGGCGGTGGCGGTCACCTGCAGTATCt                                                                                                                  >  3:192447/1‑149 (MQ=255)
                                    ccTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTg                                                                                                                 <  3:222880/149‑1 (MQ=255)
                                     cTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGTCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTgg                                                                                                                <  1:223578/149‑1 (MQ=255)
                                      tGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                               >  3:43071/1‑149 (MQ=255)
                                       gCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCGGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                               >  1:108962/1‑148 (MQ=255)
                                       gCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                               >  8:62732/1‑148 (MQ=255)
                                        cTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCa                                                                                                             <  7:131969/149‑1 (MQ=255)
                                         ttCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCAt                                                                                                            >  2:222020/1‑149 (MQ=255)
                                              cATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                               >  2:54024/1‑141 (MQ=255)
                                              cATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGc                                                                                                               <  1:54024/141‑1 (MQ=255)
                                               aTCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCt                                                                                                      <  7:36132/149‑1 (MQ=255)
                                                       tCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACt                                                                                               <  8:102049/148‑1 (MQ=255)
                                                         aGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGtt                                                                                            >  2:229951/1‑149 (MQ=255)
                                                                cgcCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGt                                                                                      >  8:227519/1‑148 (MQ=255)
                                                                        gCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGaa                                                                             <  4:43071/149‑1 (MQ=255)
                                                                                aGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGt                                                                     >  8:236383/1‑149 (MQ=255)
                                                                                 gCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGtt                                                                    <  3:23964/149‑1 (MQ=255)
                                                                                  ccAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTg                                                                   <  7:62732/149‑1 (MQ=255)
                                                                                     ggTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCa                                                                <  1:229951/149‑1 (MQ=255)
                                                                                          gCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCgg                                                            <  4:242441/148‑1 (MQ=255)
                                                                                            cACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCACCCCCTTGCCGGCTGCGTGGGGGGCGGTGCCCGGCAGGATCTGGCCATTGGTGTGGGGACTGTTGAAGGTGGCCTGGCAGGTGCCGTTTCCCCCGGGTa                                                         >  1:80428/1‑149 (MQ=255)
                                                                                                 cggccggcggcGTTGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGTCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCt                                                    <  4:48120/149‑1 (MQ=255)
                                                                                                    ccggcggcGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACTCGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCAGTGTTCAAGGTGACCTGGAAGGTGCCGTTGGCACCGGCTACCGCtgtg                                                 <  4:192447/149‑1 (MQ=255)
                                                                                                    ccggcggcGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCtgtg                                                 >  3:285124/1‑149 (MQ=255)
                                                                                                          gcgTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGc                                           <  4:285124/149‑1 (MQ=255)
                                                                                                              ggtgTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGt                                       >  8:89067/1‑149 (MQ=255)
                                                                                                               gtgtCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGtt                                      <  5:298286/149‑1 (MQ=255)
                                                                                                                            tgtAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATcg                         >  4:93606/1‑149 (MQ=255)
                                                                                                                                 ggTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcac                    >  1:146386/1‑149 (MQ=255)
                                                                                                                                 ggTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcac                    >  4:112773/1‑149 (MQ=255)
                                                                                                                                  gTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTg                                                                   >  4:67420/1‑101 (MQ=255)
                                                                                                                                  gTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTg                                                                   <  3:67420/101‑1 (MQ=255)
                                                                                                                                  gTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacg                   <  2:108962/149‑1 (MQ=255)
                                                                                                                                     cagccggcAACGTGTCGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggt                 >  3:156059/3‑148 (MQ=255)
                                                                                                                                      agccggcAACGTGTCGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtca               >  6:28915/2‑149 (MQ=255)
                                                                                                                                       gccggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcac              >  3:317289/1‑149 (MQ=255)
                                                                                                                                        ccggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcacg             >  6:314688/1‑149 (MQ=255)
                                                                                                                                        ccggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGGCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcacg             <  1:222020/149‑1 (MQ=255)
                                                                                                                                         cggcGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGcacggtcacgg            <  7:89067/149‑1 (MQ=255)
                                                                                                                                             gACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGtgct        <  2:80428/149‑1 (MQ=255)
                                                                                                                                             gACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGtgct        <  2:159003/149‑1 (MQ=255)
                                                                                                                                               cacGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGtgctgc      <  7:227519/149‑1 (MQ=255)
                                                                                                                                                acGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGCCCTGGAAGGTGCCGTTGCCACCGGATACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGTGCTGcc     >  5:251937/1‑149 (MQ=255)
                                                                                                                                                  gTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGTGCTGCCgg   >  4:320073/1‑149 (MQ=255)
                                                                                                                                                   ttGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGTGCTGCCGGc  >  3:135971/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGCGGTGTAAGGTGCGGCCGGCGACACGTTGCCGGCTGCGTCGGTGGCGGTCACCTGCAGTATCTGGCCATTGGTCTGGGCACTGTTCAAGGTGACCTGGAAGGTGCCGTTGCCACCGGCTACCGCTGTGCCCAGCACGTTGCCGCTGGCATCGCGCACGGTCACGGTGCTGCCGGC  >  NC_002947/932459‑932754

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: