Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,188,703 Δ1 bp 6.9% coding (5/312 nt) PP_1038 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,188,7030A. 6.9% 180.2 / 7.4 58coding (5/312 nt)PP_1038hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base . (2/2);  total (31/27)
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

GACGGCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGCAGCCA  >  NC_002947/1188558‑1188850
                                                                                                                                                 |                                                                                                                                                   
gACGGCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACaa                                                                                                                                                  <  4:75399/149‑1 (MQ=255)
    gCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTc                                                                                                                                              <  8:357085/149‑1 (MQ=255)
    gCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTc                                                                                                                                              <  6:286686/149‑1 (MQ=255)
    gCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTc                                                                                                                                              <  6:78746/149‑1 (MQ=255)
       gcgTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGcc                                                                                                                                           <  7:231403/149‑1 (MQ=255)
        cgTGCCCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCct                                                                                                                                          >  2:184208/1‑149 (MQ=255)
           gACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCcttct                                                                                                                                       >  7:301895/1‑149 (MQ=255)
                tcatGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTc                                                                                                                                  <  8:219591/149‑1 (MQ=255)
                tcatGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTc                                                                                                                                  >  7:329107/1‑149 (MQ=255)
                    gatgCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCg                                                                                                                              >  6:393464/1‑149 (MQ=255)
                        ccGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCgcca                                                                                                                           <  1:76831/148‑1 (MQ=255)
                          gCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCgccagcc                                                                                                                        <  5:105291/149‑1 (MQ=255)
                           cACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGcc                                                                                                                                           >  8:196254/1‑129 (MQ=255)
                           cACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGcc                                                                                                                                           <  7:196254/129‑1 (MQ=255)
                                    cGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGt                                                                                                              <  3:42758/149‑1 (MQ=255)
                                        tgtTCCGCGCTGTGCAGAACTCCTGGCGCCGGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGt                                                                                                           <  1:46699/148‑1 (MQ=255)
                                        tgtTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGcc                                                                                                                                           >  5:271115/1‑116 (MQ=255)
                                        tgtTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGcc                                                                                                                                           <  6:271115/116‑1 (MQ=255)
                                           tCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAg                                                                                                       >  8:328646/1‑149 (MQ=255)
                                                 cTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCt                                                                                                 >  2:198679/1‑149 (MQ=255)
                                                 cTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCt                                                                                                 <  8:329107/149‑1 (MQ=255)
                                                 cTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCt                                                                                                 <  1:184208/149‑1 (MQ=255)
                                                  tgtgCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCtg                                                                                                >  7:95384/1‑149 (MQ=255)
                                                      cAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTgtgt                                                                                             <  2:212032/148‑1 (MQ=255)
                                                      cAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTgtgt                                                                                             <  8:301895/148‑1 (MQ=255)
                                                      cAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGtt                                                                                            >  7:248294/1‑149 (MQ=255)
                                                            tCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGT‑‑AAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgg                                                                                    >  1:145020/1‑149 (MQ=255)
                                                            tCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGT‑‑AAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgg                                                                                    >  1:139571/1‑149 (MQ=255)
                                                                 gcgcCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgggcg                                                                                 <  4:270359/149‑1 (MQ=255)
                                                                      cGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgggcggcggg                                                                            <  5:393464/149‑1 (MQ=255)
                                                                       ggACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgggcggcgggc                                                                           <  1:214402/149‑1 (MQ=255)
                                                                        gACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGGACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCgcgggcggcgggcg                                                                          >  8:135585/1‑149 (MQ=255)
                                                                           gAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCAt                                                                       >  5:374237/1‑149 (MQ=255)
                                                                                   ggAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACt                                                               <  7:135585/149‑1 (MQ=255)
                                                                                    gAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTa                                                              >  7:166753/1‑149 (MQ=255)
                                                                                      aGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGT‑‑AAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGAc                                                          <  2:139571/149‑1 (MQ=255)
                                                                                      aGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGT‑‑AAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGAc                                                          <  2:145020/149‑1 (MQ=255)
                                                                                       gATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGa                                                           <  6:306653/149‑1 (MQ=255)
                                                                                         tgctgcGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGAc                                                          <  8:166753/148‑1 (MQ=255)
                                                                                            tgcGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACt                                                      >  5:343723/1‑149 (MQ=255)
                                                                                                        tATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCaaa                                          >  4:183003/1‑149 (MQ=255)
                                                                                                           gTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGGGACTATGACCACTGCGCCTGGCAAAccc                                       >  8:243596/1‑149 (MQ=255)
                                                                                                              ccGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGAc                                                          >  4:57031/1‑127 (MQ=255)
                                                                                                              ccGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGAc                                                          <  3:57031/127‑1 (MQ=255)
                                                                                                               cGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTgg                                   >  1:60971/1‑149 (MQ=255)
                                                                                                                  aaCGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTggc                                  >  8:46924/1‑147 (MQ=255)
                                                                                                                             gtgGGTGAACTTAGGCGTGTACAATCTCCCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTcc                     <  8:248294/149‑1 (MQ=255)
                                                                                                                                 gTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATGAGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGccc                 >  1:319009/1‑149 (MQ=255)
                                                                                                                                  tGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGCCTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCg                >  8:341831/1‑149 (MQ=255)
                                                                                                                                   gAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGt               >  2:313206/1‑149 (MQ=255)
                                                                                                                                   gAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGAGCCAGTTCCGCCCGt               >  5:393508/1‑149 (MQ=255)
                                                                                                                                        aaGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCGGTGTTCGCCGCGGGCGGCGGGCGCTTCGGTGACTATGGCCACTGCGGCTGCCAGACCCTGGGCCAGGGCCAGTTCCGCCCGATGGAc          >  8:344017/1‑149 (MQ=255)
                                                                                                                                         aGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGCAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGAc          <  8:95384/148‑1 (MQ=255)
                                                                                                                                           gCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGc       >  4:280821/1‑149 (MQ=255)
                                                                                                                                               gtACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATTGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGcagc    >  8:262627/1‑148 (MQ=255)
                                                                                                                                                tACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGcagcca  >  8:97995/1‑149 (MQ=255)
                                                                                                                                                tACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGcagc    >  2:195618/1‑147 (MQ=255)
                                                                                                                                                tACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGcagc    >  6:221820/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACGGCCGCGTGACCATCATGATGCCGCACCCGGAGCGTGTGTTCCGCGCTGTGCAGAACTCCTGGCGCCCGGACGAGTGGCAGGAAGATGCTGCGCTGATGCGTATGTTCCGCAACGCGCGGGTGTGGGTGAACTAAGGCGTGTACAAGCTCGCCTTCTTCGTCCCCGCCAGCCATGTCGAAGTGGTCAAGGCCGCTGTGTTCGCCGCGGGCGGCGGGCGCATCGGTGACTATGACCACTGCGCCTGGCAAACCCTTGGCCAGGGCCAGTTCCGCCCGTTGGACGGCAGCCA  >  NC_002947/1188558‑1188850

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: