Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,793 C→T 8.0% R154C (CGC→TGC)  PP_3138 → VirK domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,553,7930CT 8.0% 178.8 / 3.0 63R154C (CGC→TGC) PP_3138VirK domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTG  >  NC_002947/3553647‑3553939
                                                                                                                                                  |                                                                                                                                                        
aCCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGc                                                                                                                                                        <  5:111423/149‑1 (MQ=255)
 ccTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCa                                                                                                                                                       <  2:155698/149‑1 (MQ=255)
 ccTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCa                                                                                                                                                       <  8:154724/149‑1 (MQ=255)
           tgtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGt                                                                                                                                              <  1:122231/147‑1 (MQ=255)
           tgtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGt                                                                                                                                              <  1:323688/147‑1 (MQ=255)
                ccTGTCGAGCCACTGGCCGGGCGGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGGGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaaggc                                                                                                                                        <  8:120400/148‑4 (MQ=255)
                 cTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcctg                                                                                                                                      >  1:25855/1‑149 (MQ=255)
                   gTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCTGCAGCCAGGGTGAAAcctgcc                                                                                                                                    >  7:150075/1‑149 (MQ=255)
                    tCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcctgc                                                                                                                                     <  5:260298/147‑1 (MQ=255)
                           aCTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGTTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                            >  7:242352/1‑149 (MQ=255)
                            cTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTg                                                                                                                           <  1:65358/149‑1 (MQ=255)
                            cTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTg                                                                                                                           <  3:85871/149‑1 (MQ=255)
                               gCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtg                                                                                                                        <  5:149985/149‑1 (MQ=255)
                                  gggCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcgg                                                                                                                     <  3:267410/149‑1 (MQ=255)
                                      tGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGTTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCa                                                                                                                 <  8:242352/149‑1 (MQ=255)
                                      tGTAGCCTGGTGCTGGACCGGCCAATCTGGATCAAGCGCGAACGAGAGCTGGTGCTGAACCTGTTCCAGGCCGAGCTGCGTGTGGCATCCCTTGACTTCACCCTGGGCCGCAGCCAGGGGGAAAACTGCCTGTTCATTGGTGCGGTCCa                                                                                                                 >  8:246565/1‑149 (MQ=255)
                                                 gCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAgg                                                                                                                                                <  3:322117/107‑1 (MQ=255)
                                                 gCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAgg                                                                                                                                                >  4:322117/1‑107 (MQ=255)
                                                   tGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCcaca                                                                                                     >  5:141340/1‑147 (MQ=255)
                                                     gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACaa‑‑‑a                                                                                                <  2:169387/147‑1 (MQ=255)
                                                     gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑Agg                                                                                              >  1:264295/1‑149 (MQ=255)
                                                      aCCGGCCAATCTGGTTCAAGCGCGAAGGCGTGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGt                                                                                            <  1:133179/149‑1 (MQ=255)
                                                         ggCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATc                                                                                         <  7:53490/149‑1 (MQ=255)
                                                          gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATc                                                                                         >  2:255075/1‑148 (MQ=255)
                                                          gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATc                                                                                         >  2:106366/1‑148 (MQ=255)
                                                             aaTCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                            >  7:249311/1‑115 (MQ=255)
                                                             aaTCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                            <  8:249311/115‑1 (MQ=255)
                                                             aaTCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACa                                                                                     >  4:183795/1‑149 (MQ=255)
                                                               tCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCCTGTGGCATCCCTTGCCTTCACCCTGGGACGCAGCAAGGGAGCAACCTGCCTGGTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACacc                                                                                   >  3:215539/1‑147 (MQ=255)
                                                                         gcgcGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGc                                                                         >  1:323588/1‑149 (MQ=255)
                                                                               aGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCa                                                                   <  5:68176/149‑1 (MQ=255)
                                                                                ggCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGt                                                                                                                                 >  6:42371/1‑91 (MQ=255)
                                                                                ggCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGt                                                                                                                                 <  5:42371/91‑1 (MQ=255)
                                                                                              gAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctga                                                    >  1:174364/1‑149 (MQ=255)
                                                                                               aaCCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctgac                                                   <  2:264295/149‑1 (MQ=255)
                                                                                               aaCCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctgac                                                   >  8:290521/1‑149 (MQ=255)
                                                                                                       ccAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACt                                           >  5:44048/1‑149 (MQ=255)
                                                                                                         aGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGCTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACtt                                          <  1:255075/149‑1 (MQ=255)
                                                                                                         aGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTc                                         >  3:190446/1‑149 (MQ=255)
                                                                                                         aGGCCGACCTGCGTGGGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTc                                         <  3:183795/149‑1 (MQ=255)
                                                                                                           accGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcg                                                                                                                      >  2:31965/2‑75 (MQ=255)
                                                                                                            ccGACCTGCGTGTGGCATACGTTGCCTTCAGCGTTGGCAGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGGGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCGATCTACCTTTACCTGACCAAGGACTTCGaa                                      <  4:215539/149‑1 (MQ=255)
                                                                                                                 cTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCt                                 >  8:221549/1‑149 (MQ=255)
                                                                                                                   gCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGc                               >  4:732/1‑149 (MQ=255)
                                                                                                                     gggtgGCATCCCTTGCATTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCgg                             <  8:249685/147‑1 (MQ=255)
                                                                                                                       gttgCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGc                            >  1:67869/4‑148 (MQ=255)
                                                                                                                       gtgGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGc                            >  8:237313/1‑148 (MQ=255)
                                                                                                                          gCATCCCTTGCCTTCACCCTGGGCTGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCAc                         <  8:150075/148‑1 (MQ=255)
                                                                                                                                cTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCt                  >  8:52026/1‑149 (MQ=255)
                                                                                                                                   gCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCATACACAACCACAACCTCAGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGc                    <  1:106366/148‑1 (MQ=255)
                                                                                                                                   gCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGtt               >  1:221928/1‑149 (MQ=255)
                                                                                                                                   gCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGt                <  6:44048/148‑1 (MQ=255)
                                                                                                                                    gtttCACCCTGGCTGCGGCCCAGGGTGAAGCCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGaa                                      <  3:139245/123‑1 (MQ=255)
                                                                                                                                    gtttCACCCTGGCTGCGGCCCAGGGTGAAGCCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGaa                                      >  4:139245/3‑125 (MQ=255)
                                                                                                                                     cTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGc             >  8:190330/1‑149 (MQ=255)
                                                                                                                                      ttCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCt            >  3:109813/1‑149 (MQ=255)
                                                                                                                                      ttCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCt            >  3:89225/1‑149 (MQ=255)
                                                                                                                                      ttCACACTGGGAAGCAGCCAGGGTCAAACCTGCCTGTTCATTGGTGCGGTGCAGGCCATCCACAA‑‑‑AGGTATCGACAGCGAAACctctct                                                                        >  7:266625/1‑89 (MQ=255)
                                                                                                                                         aCCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGa         <  2:174364/149‑1 (MQ=255)
                                                                                                                                          cccTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCg          >  2:279219/1‑147 (MQ=255)
                                                                                                                                           ccTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAg       >  4:213719/1‑149 (MQ=255)
                                                                                                                                               ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGtt   <  7:290521/149‑1 (MQ=255)
                                                                                                                                               ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGt    <  3:732/148‑1 (MQ=255)
                                                                                                                                                 ccGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTg  >  5:339006/1‑148 (MQ=255)
                                                                                                                                                 ccGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTg  >  1:212938/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                        
ACCTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAA‑‑‑AGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTG  >  NC_002947/3553647‑3553939

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: