Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 953,552 T→A 6.8% Y121N (TAT→AAT) ‡ cyoC → cytochrome bo terminal oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947953,5520TA 6.8% 170.1 / 3.1 59Y121N (TAT→AAT) ‡cyoCcytochrome bo terminal oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base A (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GTAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATC  >  NC_002947/953405‑953692
                                                                                                                                                    |                                                                                                                                            
gTAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTa                                                                                                                                             >  1:111420/1‑149 (MQ=255)
gTAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAACATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTa                                                                                                                                             >  4:207163/1‑149 (MQ=255)
 tAGCTGTTGAATCGCTGTTCCTGCTGGTGTCCTCGATCACGTTGGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTAt                                                                                                                                            <  7:322196/149‑1 (MQ=255)
 tAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTAt                                                                                                                                            <  8:287967/149‑1 (MQ=255)
 tAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCAGTTCTGTGCAGGCTTCATCGCGATGGAAATCTAt                                                                                                                                            >  6:230714/1‑149 (MQ=255)
     tGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGTGTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAg                                                                                                                                         >  1:101885/1‑149 (MQ=255)
      gTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGtt                                                                                                                                       >  6:219438/1‑149 (MQ=255)
           aaCGCTTTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGTTGTTCGATGGCAAGAAAGCAGGTGTACTTGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAACTCATAGATTTCCATc                                                                                                                                  <  7:357089/149‑1 (MQ=255)
           aaCGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCAt                                                                                                                                   <  3:233926/148‑1 (MQ=255)
             cGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCAc                                                                                                                                <  7:107527/149‑1 (MQ=255)
             cGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCAc                                                                                                                                <  6:182996/149‑1 (MQ=255)
              gCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAACTCATAGATTTCCATcgcg                                                                                                                               <  2:42224/149‑4 (MQ=255)
                      tgctgctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGc                                                                                                                       >  4:220246/1‑149 (MQ=255)
                        ctgctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTACAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGa                                                                                                                     <  3:377168/149‑1 (MQ=255)
                           ctgTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                    <  6:162453/147‑1 (MQ=255)
                              tCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCtt                                                                                                               >  3:110969/1‑149 (MQ=255)
                                   gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATg                                                                                                                                           <  4:362873/116‑1 (MQ=255)
                                   gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATg                                                                                                                                           >  3:362873/1‑116 (MQ=255)
                                   gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                    >  5:296404/1‑139 (MQ=255)
                                   gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAg                                                                                                                    <  6:296404/139‑1 (MQ=255)
                                   gATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGcc                                                                                                          >  8:344205/1‑149 (MQ=255)
                                    aTCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGccc                                                                                                         <  4:110969/149‑1 (MQ=255)
                                    aTCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGccc                                                                                                         >  1:71146/1‑149 (MQ=255)
                                       aCCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCAGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCcgca                                                                                                      <  5:230714/149‑1 (MQ=255)
                                              tcctCGCCATGCTGTAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGTAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGt                                                                                               <  8:262396/146‑1 (MQ=255)
                                              gCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGt                                                                                               >  4:59935/1‑149 (MQ=255)
                                                    ccATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                         <  3:105045/149‑1 (MQ=255)
                                                             aaGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCgg                                                                                  <  3:207163/148‑1 (MQ=255)
                                                                           cAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTgg                                                                   >  7:160237/1‑149 (MQ=255)
                                                                            aaGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTggt                                                                  <  7:344205/149‑1 (MQ=255)
                                                                                  gCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTAc                                                            <  2:71146/149‑1 (MQ=255)
                                                                                      gtgtACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCAc                                                        >  5:72651/1‑149 (MQ=255)
                                                                                      gtgtACTGGGCTGGTTGGCTGTGACCTTCCTGTACGGTGCACGCTTCATCGCTATGGAAATCGAGGAGTTCCATCACCTGATCAGCTAGGGCTTAGGACCGCAGCGCAGTGGCTTCATGTCGGGCTTCTTCGCCCTGGTAGTTACCCAc                                                        >  5:399675/1‑149 (MQ=255)
                                                                                         tACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACgg                                                      <  3:220246/148‑1 (MQ=255)
                                                                                         tACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGt                                                     >  6:70150/1‑149 (MQ=255)
                                                                                           cTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                         >  2:343098/1‑111 (MQ=255)
                                                                                           cTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTc                                                                                         <  1:343098/111‑1 (MQ=255)
                                                                                              ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTg                                                  <  8:272807/147‑1 (MQ=255)
                                                                                              ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCa                                                >  4:76159/1‑149 (MQ=255)
                                                                                              ggCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCa                                                <  5:219438/149‑1 (MQ=255)
                                                                                                 tGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGt                                             <  3:59935/149‑1 (MQ=255)
                                                                                                 gggTTGGCGGTGACCTTCCGGTTCGGTGCAGGCTTCATCGCTATGGAAATCTCTGAGTGCCATAACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTTCACGt                                             <  4:385172/148‑1 (MQ=255)
                                                                                                        cTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCTCTGGTAGGTACCCACGGTCTGCACGTGActgct                                      >  7:29033/1‑149 (MQ=255)
                                                                                                          gtgACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTgg                                    >  2:291014/1‑149 (MQ=255)
                                                                                                               cTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGa                               >  7:61410/1‑149 (MQ=255)
                                                                                                               cTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGa                               >  1:138575/1‑149 (MQ=255)
                                                                                                                  ccTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGAt                              <  3:76159/147‑1 (MQ=255)
                                                                                                                      ttCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGa                         <  2:111420/148‑1 (MQ=255)
                                                                                                                      ttCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGAt                        <  8:160237/149‑1 (MQ=255)
                                                                                                                        cGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATg                       <  6:98600/148‑1 (MQ=255)
                                                                                                                         ggTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATgg                      <  8:110153/148‑1 (MQ=255)
                                                                                                                             cAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAAt                  >  4:365559/1‑148 (MQ=255)
                                                                                                                             cAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATc                 <  2:259089/149‑1 (MQ=255)
                                                                                                                                    cATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCTCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCatgat            <  8:29033/147‑1 (MQ=255)
                                                                                                                                     aTCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTa         >  5:357649/1‑149 (MQ=255)
                                                                                                                                     aTCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACCGCTGGTCTGATCTGGATGGCAATCATGATGTa         >  2:166422/1‑149 (MQ=255)
                                                                                                                                       ggtgATGGAACTCATAGATTTACATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTAc        >  2:305044/4‑148 (MQ=255)
                                                                                                                                        gtgATGGAACTCATAGATTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGTTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTAc        >  4:213885/3‑147 (MQ=255)
                                                                                                                                             ggAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATc  >  7:327976/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GTAGCTGTTGAAACGCTGTTCCTGCTGCTGTCCTCGATCACCTTCGGCTTCGCCATGCTGAAGATGTTCGATGGCAAGAAAGCAGGTGTACTGGGCTGGTTGGCTGTGACCTTCCTGTTCGGTGCAGGCTTCATCGCGATGGAAATCTATGAGTTCCATCACCTGATCAGCGAGGGCTTCGGCCCGCAGCGCAGTGGCTTCCTGTCGGGCTTCTTCGCCCTGGTAGGTACCCACGGTCTGCACGTGACTGCTGGTCTGATCTGGATGGCAATCATGATGTACCAGATC  >  NC_002947/953405‑953692

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: