Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 676,102 T→C 7.9% Q103Q (CAA→CAG)  PP_0578 ← putative Methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947676,1020TC 7.9% 176.0 / 5.2 63Q103Q (CAA→CAG) PP_0578putative Methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (2/3);  total (31/32)
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.12e-01

GATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGCCAT  >  NC_002947/675954‑676248
                                                                                                                                                    |                                                                                                                                                  
gATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCGCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTCGCGCTCGGTCAGt                                                                                                                                                    >  8:61736/1‑149 (MQ=255)
  tAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTg                                                                                                                                                  <  5:351088/149‑1 (MQ=255)
  tAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTg                                                                                                                                                  <  7:154064/149‑1 (MQ=255)
  tAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTg                                                                                                                                                  >  6:351088/1‑149 (MQ=255)
      ccAGGTCGTCGCGGTTGTAGAAGGCCAGCTATAACCAACCGCCGGGGGGCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCGAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCg                                                                                                                                              <  2:222184/149‑1 (MQ=255)
        aGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCgag                                                                                                                                            <  5:37454/149‑1 (MQ=255)
         gggCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGTGCTCGGCCAACTGACCgagc                                                                                                                                           <  1:99944/149‑2 (MQ=255)
         gggCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGCCAACTGACCgagc                                                                                                                                           <  6:79778/149‑2 (MQ=255)
           gCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCgaggag                                                                                                                                         <  7:183832/149‑1 (MQ=255)
                gcgGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCtt                                                                                                                                    <  2:298705/149‑1 (MQ=255)
                       tAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGGGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAgg                                                                                                                              <  7:311688/148‑1 (MQ=255)
                          aaGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCAGTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGgcgc                                                                                                                          >  7:195626/1‑149 (MQ=255)
                            ggCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCgaggag                                                                                                                                         <  8:270114/132‑1 (MQ=255)
                            ggCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCgaggag                                                                                                                                         >  7:270114/1‑132 (MQ=255)
                              ccAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTg                                                                                                                       >  4:102125/1‑148 (MQ=255)
                                        ccaaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcg                                                                                                            >  7:101612/1‑149 (MQ=255)
                                        ccaaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTCGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcg                                                                                                            >  1:117877/1‑149 (MQ=255)
                                         caaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgc                                                                                                           >  7:137777/1‑149 (MQ=255)
                                         caaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgc                                                                                                           >  1:162762/1‑149 (MQ=255)
                                          aaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgcg                                                                                                          >  3:168255/1‑149 (MQ=255)
                                          aaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgcg                                                                                                          <  2:117877/149‑1 (MQ=255)
                                              gccgGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGAt                                                                                                                    <  7:212595/135‑1 (MQ=255)
                                              gccgGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGAt                                                                                                                    >  8:212595/1‑135 (MQ=255)
                                                       cGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGt                                                                                                        >  4:70370/1‑138 (MQ=255)
                                                       cGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGt                                                                                                        <  3:70370/138‑1 (MQ=255)
                                                          gAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCgg                                                                                          <  3:23088/149‑1 (MQ=255)
                                                          gAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCgg                                                                                          >  5:260749/1‑149 (MQ=255)
                                                          gAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGCTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCgg                                                                                          >  8:214511/1‑149 (MQ=255)
                                                             cTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCa                                                                                        >  3:91362/1‑148 (MQ=255)
                                                              tGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCTGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCa                                                                                        >  8:185265/1‑147 (MQ=255)
                                                               ggTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgcg                                                                                                          >  4:350808/1‑128 (MQ=255)
                                                                gTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAg                                                                                     <  8:137777/148‑1 (MQ=255)
                                                                gTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTGCGGCAAAGc                                                                                    >  4:225188/1‑149 (MQ=255)
                                                                   cagcaCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGCTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTg                                                                                 <  7:214511/149‑1 (MQ=255)
                                                                     gcaCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGcc                                                                               >  3:369965/1‑149 (MQ=255)
                                                                      caCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGccc                                                                              >  5:192544/1‑149 (MQ=255)
                                                                      caCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGccc                                                                              >  4:345031/1‑149 (MQ=255)
                                                                       aCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGt                                                                                  >  4:112917/1‑144 (MQ=255)
                                                                       aCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGt                                                                                  <  3:112917/144‑1 (MQ=255)
                                                                                   ctctCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGa                                                                  <  6:260749/148‑1 (MQ=255)
                                                                                   ctctCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGa                                                                  <  6:192544/148‑1 (MQ=255)
                                                                                   ctctCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGa                                                                  <  4:168255/148‑1 (MQ=255)
                                                                                      tcGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCAt                                                              >  2:59720/1‑149 (MQ=255)
                                                                                        gggTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCAt                                                              <  4:369965/147‑1 (MQ=255)
                                                                                         ggTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATggg                                                           >  3:103826/1‑149 (MQ=255)
                                                                                            tCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGgcgc                                                        <  8:195626/149‑1 (MQ=255)
                                                                                               gccagccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGcgccg                                                      >  8:235335/1‑148 (MQ=255)
                                                                                               gccagccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGcgccg                                                      >  5:219747/1‑148 (MQ=255)
                                                                                                  agccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTg                                                                                                     <  6:363303/98‑1 (MQ=255)
                                                                                                  agccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTg                                                                                                     >  5:363303/1‑98 (MQ=255)
                                                                                                     caTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTc                                               <  8:23427/149‑1 (MQ=255)
                                                                                                           aGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGGCTTGCCAGGGCGCCTGGATGAAAGTCGCGGGCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAg                                         <  7:61736/149‑1 (MQ=255)
                                                                                                              aCGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAg                                                   <  2:266450/136‑1 (MQ=255)
                                                                                                              aCGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAg                                                   >  1:266450/1‑136 (MQ=255)
                                                                                                                    tGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCtg                                <  1:59720/149‑1 (MQ=255)
                                                                                                                        agagGACCAGGTCGTAGCGCTCGGTCATTTGGCCGAGGAGGTCTTGGCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCtggtgg                            <  4:105662/149‑1 (MQ=255)
                                                                                                                         gagCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGc                           >  1:320434/1‑149 (MQ=255)
                                                                                                                             aCCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCcgcg                       >  8:59548/1‑149 (MQ=255)
                                                                                                                                 ggTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCg                   <  8:25111/149‑1 (MQ=255)
                                                                                                                                        gcgcTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGgccagc             >  4:295017/1‑148 (MQ=255)
                                                                                                                                        ctccTCGGCCAACTGACCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGgccagcc            >  1:145407/4‑149 (MQ=255)
                                                                                                                                          gcTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCcac          <  3:225188/149‑1 (MQ=255)
                                                                                                                                                tCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGcc    <  3:102125/149‑1 (MQ=255)
                                                                                                                                                  aGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGCCAt  <  7:235335/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGCCAT  >  NC_002947/675954‑676248

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: ‑