Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,793 C→T 7.8% 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 7.8% 185.6 / 4.2 64R154C (CGC→TGC) PP_3138VirK domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (2/3);  total (31/33)
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.24e-01

TGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTT  >  NC_002947/3553656‑3553938
                                                                                                                                         |                                                                                                                                                  
tgtgtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGt                                                                                                                                        >  1:366203/1‑149 (MQ=255)
tgtgtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGt                                                                                                                                        <  8:237025/149‑1 (MQ=255)
   gtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGcc                                                                                                                                                    >  1:49298/1‑135 (MQ=255)
   gtgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGcc                                                                                                                                                    <  2:49298/135‑1 (MQ=255)
    tgACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                     <  7:337321/148‑1 (MQ=255)
             cAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGt                                                                                                                           <  3:174139/149‑1 (MQ=255)
             cAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAACGCGAGCTGGTGCTGAACCTGTTCCAGGCCGCCATGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCGGCgtgt                                                                                                                           >  2:225302/1‑149 (MQ=255)
              aGCCACGCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGTCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGtt                                                                                                                          <  1:372189/149‑1 (MQ=255)
                  aCTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                      <  4:143620/149‑1 (MQ=255)
                  aCTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCAtt                                                                                                                      >  6:246173/1‑149 (MQ=255)
                     cGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggt                                                                                                                   <  2:89292/149‑1 (MQ=255)
                       ccGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgc                                                                                                                 <  8:313821/149‑1 (MQ=255)
                         gggCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggtgcgg                                                                                                               >  3:75032/1‑149 (MQ=255)
                             tGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGTCCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCa                                                                                                           <  5:100/149‑1 (MQ=255)
                               ttGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAgg                                                                                                         >  8:73207/1‑149 (MQ=255)
                               ttGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCATCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggt                                                                                                                   >  4:246886/1‑139 (MQ=255)
                               ttGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCATCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTggt                                                                                                                   <  3:246886/139‑1 (MQ=255)
                                      gTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATc                                                                                                    <  7:73207/147‑1 (MQ=255)
                                          tGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCcaca                                                                                                >  7:258293/1‑147 (MQ=255)
                                             aCCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                     >  7:76273/1‑107 (MQ=38)
                                             aCCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGCTGCGGCCCAGGGTGaa                                                                                                                                     <  8:76273/107‑1 (MQ=38)
                                              ccGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTa                                                                                          >  4:289078/1‑149 (MQ=255)
                                              ccGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTa                                                                                          <  4:75032/149‑1 (MQ=255)
                                                 gCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATc                                                                                        >  8:81014/1‑148 (MQ=255)
                                                  ccAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAgg                                                                                                         >  2:235002/1‑130 (MQ=255)
                                                  ccAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAgg                                                                                                         <  1:235002/130‑1 (MQ=255)
                                                     aTCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAg                                                                                   <  7:81014/149‑1 (MQ=255)
                                                              aaGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTc                                                                          >  3:329547/1‑149 (MQ=255)
                                                                 cgcgAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCt                                                                       >  5:214529/1‑149 (MQ=255)
                                                                  gcgAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTg                                                                      >  5:201486/1‑149 (MQ=255)
                                                                    gAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGc                                                                    >  7:154467/1‑149 (MQ=255)
                                                                       gtcGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCTCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCAt                                                                 <  1:225302/147‑1 (MQ=255)
                                                                        gCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcctgc                                                                                                                               >  2:373156/1‑86 (MQ=255)
                                                                        gCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAAcctgc                                                                                                                               <  1:373156/86‑1 (MQ=255)
                                                                           aGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTAc                                                             <  2:213937/149‑1 (MQ=255)
                                                                              tGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGt                                                          >  5:242233/1‑149 (MQ=255)
                                                                                gTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtga                                                        <  6:201486/149‑1 (MQ=255)
                                                                                gTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtg                                                         <  8:154467/148‑1 (MQ=255)
                                                                                 tGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgac                                                       <  5:246173/149‑1 (MQ=255)
                                                                                 tGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgac                                                       <  6:107563/149‑1 (MQ=255)
                                                                                      aaCCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGtgacctgac                                                  <  6:214529/149‑1 (MQ=255)
                                                                                        ccTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTg                                                                      >  6:120643/1‑127 (MQ=255)
                                                                                        ccTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTg                                                                      <  5:120643/127‑1 (MQ=255)
                                                                                         cTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCaa                                               <  8:258293/149‑1 (MQ=255)
                                                                                           gTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACACTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAgg                                             >  4:297165/1‑149 (MQ=255)
                                                                                             tCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGa                                            <  6:408976/148‑1 (MQ=255)
                                                                                              ccAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGAc                                           <  2:366203/148‑1 (MQ=255)
                                                                                               cAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACtt                                         >  4:145450/1‑149 (MQ=255)
                                                                                                  gCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCg                                       >  8:240851/1‑148 (MQ=255)
                                                                                                       catGCGTGTGGAATCCCTTGCCTTCACCATGGGACGCAGCCACGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATACACAAAGGTATCGACAGCGAAACCTAGCTGGGCATCTACCGTGACCTGACCAAGTACTTCGaattgc                                 >  7:333879/3‑145 (MQ=255)
                                                                                                             tgtgGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGc                           >  5:31593/1‑149 (MQ=255)
                                                                                                              gtgGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGc                           >  6:47979/1‑148 (MQ=255)
                                                                                                                       cTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTgg                                                                     <  6:93291/97‑1 (MQ=255)
                                                                                                                       cTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTgg                                                                     >  5:93291/1‑97 (MQ=255)
                                                                                                                          gCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGtt              >  1:145209/1‑149 (MQ=255)
                                                                                                                            cTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGc            >  3:399699/1‑149 (MQ=255)
                                                                                                                            cTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGc            >  8:151965/1‑149 (MQ=255)
                                                                                                                             ttCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTAGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCt           >  1:77137/1‑149 (MQ=255)
                                                                                                                             ttCACCCTGGCTGCGGCCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATc                                                                <  3:298594/96‑1 (MQ=38)
                                                                                                                             ttCACCCTGGCTGCGGCCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATc                                                                >  4:298594/1‑96 (MQ=38)
                                                                                                                                   cTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGc     <  4:329547/149‑1 (MQ=255)
                                                                                                                                   cTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGc     <  6:242233/149‑1 (MQ=255)
                                                                                                                                     gggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGt   >  1:66866/1‑149 (MQ=255)
                                                                                                                                      ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGtt  <  7:240851/149‑1 (MQ=255)
                                                                                                                                      ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGtt  <  5:47979/149‑1 (MQ=255)
                                                                                                                                      ggCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGt   <  3:289078/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
TGTGTGACCTGTCCAGCCACTCGCCGGGCTGTTGCCTGGTGCTGGACCGGCCAATCTGGTTCAAGCGCGAAGGCGAGCTGGTGCTGAACCTGTTCCAGGCCGACCTGCGTGTGGCATCCCTTGCCTTCACCCTGGGCCGCAGCCAGGGTGAAACCTGCCTGTTCATTGGTGCGGTGCAGGGCATCCACAAAGGTATCGACAGCGAAACCTCGCTGGCCATCTACCGTGACCTGACCAAGGACTTCGAAGGGCTGCGGCCACGCAGCCTGTTGCTCGAAGCGTT  >  NC_002947/3553656‑3553938

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: