Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,796 A→G 6.4% S155G (AGC→GGC)  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,7960AG 6.4% 183.3 / 2.4 62S155G (AGC→GGC) PP_3138VirK domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base G (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTGAA  >  NC_002947/3553649‑3553941
                                                                                                                                                   |                                                                                                                                                  
cTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTCGTGCTGCACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAg                                                                                                                                                   >  4:302302/1‑149 (MQ=255)
     gCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAgg                                                                                                                                             >  2:144972/1‑149 (MQ=255)
           tgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGaaa                                                                                                                                       <  4:137796/149‑1 (MQ=255)
           tgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                        <  8:130574/148‑1 (MQ=255)
              ccTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcct                                                                                                                                    <  4:113092/149‑1 (MQ=255)
               cTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAgg                                                                                                                                             <  1:254133/139‑1 (MQ=255)
               cTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAgg                                                                                                                                             >  2:254133/1‑139 (MQ=255)
                tGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcctgc                                                                                                                                  >  1:265446/1‑149 (MQ=255)
                    cAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGt                                                                                                                              >  4:95956/1‑149 (MQ=255)
                       ccACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCa                                                                                                                           >  8:213021/1‑149 (MQ=255)
                         aCTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                         >  1:247542/1‑149 (MQ=255)
                                  gCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                        <  8:79513/125‑1 (MQ=255)
                                  gCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                        >  7:79513/1‑125 (MQ=255)
                                      ttGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAgg                                                                                                            >  7:280482/1‑149 (MQ=255)
                                          cTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCAt                                                                                                        >  5:211299/1‑149 (MQ=255)
                                           tGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATc                                                                                                       >  6:231255/1‑149 (MQ=255)
                                             gTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCca                                                                                                     <  3:95956/149‑1 (MQ=255)
                                                   gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcggt                                                                                                                 <  5:33403/131‑1 (MQ=255)
                                                   gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcggt                                                                                                                 >  6:33403/1‑131 (MQ=255)
                                                   gACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAg                                                                                                <  7:213021/148‑1 (MQ=255)
                                                        gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcggt                                                                                                                 >  6:66280/1‑126 (MQ=255)
                                                        gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcggt                                                                                                                 <  5:66280/126‑1 (MQ=255)
                                                        gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATc                                                                                                       >  3:226472/1‑136 (MQ=255)
                                                        gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATc                                                                                           >  7:161275/1‑148 (MQ=255)
                                                             tCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTc                                                                                                                            >  8:101102/1‑110 (MQ=255)
                                                             tCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTc                                                                                                                            <  7:101102/110‑1 (MQ=255)
                                                             tCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGc                                                                                     >  2:250210/1‑149 (MQ=255)
                                                                   tCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAggg                                                                                                           <  2:148880/121‑1 (MQ=255)
                                                                    cAAGCGCGAAGGCGAGCTGGTGCTGAACCAGTTCCAGTCCGACCTGCGTGTGGCATACCTTGCCTTCACCATGGGCCGCAGCCAGGGTGAAACCTGCCTGATCATTGGTGCGGTGCACGGCATCAAAAAAGGTATCCACAGCTAAACCt                                                                              >  1:243212/1‑149 (MQ=255)
                                                                     aaGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTc                                                                             >  1:162565/1‑149 (MQ=255)
                                                                      aGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGtt                                                                                                                             >  1:286748/1‑100 (MQ=255)
                                                                      aGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGtt                                                                                                                             <  2:286748/100‑1 (MQ=255)
                                                                      aGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCg                                                                            >  4:32066/1‑149 (MQ=255)
                                                                       gcgcGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCATGGCATCCACAAAGGTATCGACAGCGAAACCTCGc                                                                           >  6:121216/1‑149 (MQ=255)
                                                                         gcgAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTg                                                                         <  1:144972/149‑1 (MQ=255)
                                                                             aGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCa                                                                     >  1:133178/1‑149 (MQ=255)
                                                                                  aGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTAc                                                                <  2:247542/149‑1 (MQ=255)
                                                                                       gTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtga                                                           <  7:21609/149‑1 (MQ=255)
                                                                                       gTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtga                                                           >  5:270137/1‑149 (MQ=255)
                                                                                           tGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacct                                                        >  2:134678/1‑148 (MQ=255)
                                                                                             aaCCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCATGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctga                                                      <  5:121216/148‑1 (MQ=255)
                                                                                             aaCCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctgac                                                     >  2:101020/1‑149 (MQ=255)
                                                                                                cTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacatgacc                                                    >  3:230599/1‑147 (MQ=255)
                                                                                                   ttCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGATCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAg                                                 <  1:101020/147‑1 (MQ=255)
                                                                                                      cAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACtt                                            >  7:258358/1‑149 (MQ=255)
                                                                                                               cTGCGTGTGGCATCCCTTGCCTTCTCCCTGGGCCGCAGACAGGGTTAAACCAGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCTAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGc                                    <  3:302302/148‑1 (MQ=255)
                                                                                                               cTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCt                                   >  7:261266/1‑149 (MQ=255)
                                                                                                                          aTCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACg                          >  6:213590/1‑147 (MQ=255)
                                                                                                                          aTCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACg                          >  3:182045/1‑147 (MQ=255)
                                                                                                                            cccTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATTCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAg                                       >  5:202023/1‑132 (MQ=255)
                                                                                                                            cccTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATTCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAg                                       <  6:206715/132‑1 (MQ=255)
                                                                                                                            cccTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATTCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAg                                       <  6:202023/132‑1 (MQ=255)
                                                                                                                            cccTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATTCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAg                                       >  5:206715/1‑132 (MQ=255)
                                                                                                                              gcaggtttCACCCTGGCTGCGGCCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCt                    >  7:65166/7‑149 (MQ=255)
                                                                                                                              cTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCt                    <  4:182045/149‑1 (MQ=255)
                                                                                                                                 gCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTg                   <  8:280482/147‑1 (MQ=255)
                                                                                                                                   cTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGc               >  2:278827/1‑149 (MQ=255)
                                                                                                                                    ttCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGc                      <  3:215023/141‑1 (MQ=255)
                                                                                                                                    ttCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGc                      >  4:215023/1‑141 (MQ=255)
                                                                                                                                      cACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTc             <  2:265446/148‑1 (MQ=255)
                                                                                                                                             ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGt      <  8:261266/148‑1 (MQ=255)
                                                                                                                                                cGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTGaa  <  3:32066/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGAGGCTGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTTGAA  >  NC_002947/3553649‑3553941

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: