Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,446,628 C→T 8.6% P341L (CCG→CTG)  gpD → tail formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,446,6280CT 8.6% 161.7 / 4.8 59P341L (CCG→CTG) gpDtail formation
Reads supporting (aligned to +/- strand):  ref base C (31/22);  new base T (2/3);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

TTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATG  >  NC_002947/3446483‑3446776
                                                                                                                                                 |                                                                                                                                                     
ttCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCAAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCt                                                                                                                                                    >  4:284329/1‑149 (MQ=255)
 tCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGCCACCAAGCAGCt                                                                                                                                                    >  4:116379/1‑148 (MQ=255)
           gCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGagc                                                                                                                                         <  3:219366/149‑1 (MQ=255)
           gCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGagc                                                                                                                                         <  6:160359/149‑1 (MQ=255)
            cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGa                                                                                                                                           <  3:99200/146‑1 (MQ=255)
            cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGa                                                                                                                                           >  4:99200/1‑146 (MQ=255)
            cTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTCag                                                                                                                                          >  1:137214/1‑147 (MQ=255)
               gATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCt                                                                                                                                     <  8:28675/149‑1 (MQ=255)
                 tGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCTGCTCAACCTTGAGCGGCttg                                                                                                                                   <  7:213699/149‑3 (MQ=255)
                          ccTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACg                                                                                                                          >  2:294376/1‑149 (MQ=255)
                           cTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACgg                                                                                                                         >  3:105089/1‑149 (MQ=255)
                           cTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACgg                                                                                                                         <  1:186629/149‑1 (MQ=255)
                               tggACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGt                                                                                                                     <  2:137214/149‑1 (MQ=255)
                                    ggggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAACGACCCCAAGCCGCTCAAGCTTGAGAAGCTGTAACCCCCAGTCCGTacacc                                                                                                                >  8:295267/1‑148 (MQ=255)
                                      ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                              >  8:200796/1‑149 (MQ=255)
                                      ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                              >  8:251278/1‑149 (MQ=255)
                                          gAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCAt                                                                                                           >  1:100041/1‑148 (MQ=255)
                                            gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                        >  7:194595/1‑149 (MQ=255)
                                            gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                        >  4:282951/1‑149 (MQ=255)
                                            gCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCAc                                                                                                        >  4:19085/1‑149 (MQ=255)
                                              aGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTg                                                                                                      <  3:38464/149‑1 (MQ=255)
                                                gTATTTACCCAGGCCGGCTGGACTACGACCGTCGAGTCCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCCGGCAAGAAAAACAAAGACCCCACGCCGCGCACGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga                                                                                                    <  8:357450/149‑1 (MQ=255)
                                                  aTTTACACAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACACCACGGCCGTACACGGCCATCACTGgag                                                                                                   >  2:249108/1‑148 (MQ=255)
                                                        ccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa                                                                                            >  3:80784/1‑149 (MQ=255)
                                                               ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGaa                                                                                                 >  2:93703/1‑137 (MQ=255)
                                                               ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGaa                                                                                                 <  1:93703/137‑1 (MQ=255)
                                                               ggCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTAt                                                                                     <  3:19085/149‑1 (MQ=255)
                                                                            ccGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGAAGCTGTAACCAAACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCAATATCAGTTCAACAg                                                                        >  5:236580/1‑149 (MQ=255)
                                                                               tCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTg                                                                     >  2:1131/1‑149 (MQ=255)
                                                                                 gAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCCCGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGc                                                                    >  6:112140/1‑148 (MQ=255)
                                                                                        aCGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGAtc                                                            <  1:294376/149‑1 (MQ=255)
                                                                                           gcggcAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGAtcctc                                                         <  3:284329/149‑1 (MQ=255)
                                                                                                aagaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGCCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCaa                                                   >  1:333273/1‑149 (MQ=255)
                                                                                                 agaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAAc                                                  <  4:105089/149‑1 (MQ=255)
                                                                                                  gaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACg                                                 >  6:219340/1‑149 (MQ=255)
                                                                                                  gaagGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACg                                                 >  6:224426/1‑149 (MQ=255)
                                                                                                    agGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACgcc                                               >  1:298634/1‑149 (MQ=255)
                                                                                                       gCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCAACACTGGAGAAACCAATGGCTATCTCAGTTCACCAGCTGCAACACATCCTCCCCAACgccggc                                            >  6:306483/1‑149 (MQ=255)
                                                                                                        caaggccaaggcCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACgccggcc                                           <  2:100041/149‑1 (MQ=255)
                                                                                                          aggccaaggccTCGGGCAAGAAAAAGAAAGAGACCAAGCCTCTCAAGGTTGAGCAGCTGTAACCCCACGGGCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGc                                         <  8:340209/149‑1 (MQ=255)
                                                                                                             ccaaggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGt                                                                               >  4:52708/1‑109 (MQ=255)
                                                                                                             ccaaggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGt                                                                               <  3:52708/109‑1 (MQ=255)
                                                                                                                    cTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCg                               >  6:33466/1‑149 (MQ=255)
                                                                                                                    cTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGACCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCg                               >  2:304212/1‑149 (MQ=255)
                                                                                                                       gggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCTCAAAGCCGGCGttt                            <  5:112140/149‑1 (MQ=255)
                                                                                                                        ggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTa                           <  3:282951/149‑1 (MQ=255)
                                                                                                                               aaaagaaaGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                    >  1:21080/1‑149 (MQ=255)
                                                                                                                                aacgaaagaCACCAAGCCGCTCAAGGATGAGCAGCTGTAACCCCACCGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCg                   <  8:356075/146‑1 (MQ=255)
                                                                                                                                 aagaaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCg                   <  1:249108/149‑1 (MQ=255)
                                                                                                                                  agaaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGa                 <  1:1131/149‑1 (MQ=255)
                                                                                                                                      agaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCa             <  8:194595/149‑1 (MQ=255)
                                                                                                                                          ttacaGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCgg                                 <  6:178450/121‑1 (MQ=255)
                                                                                                                                          ttacaGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCgg                                 >  5:178450/5‑125 (MQ=255)
                                                                                                                                            cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGcta       >  1:231678/1‑149 (MQ=255)
                                                                                                                                            cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGcta       >  1:228836/1‑149 (MQ=255)
                                                                                                                                            cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGct        >  5:225925/1‑148 (MQ=255)
                                                                                                                                            cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGct        >  5:179803/1‑148 (MQ=255)
                                                                                                                                              aGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                    <  1:354775/134‑1 (MQ=255)
                                                                                                                                              aGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTccc                    >  2:354775/1‑134 (MQ=255)
                                                                                                                                                 cGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATg  <  2:298634/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTCAAGCCGGGGCTCGATGGCGAGTACCTGGTGGACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACAGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATG  >  NC_002947/3446483‑3446776

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