Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,765,218 Δ1 bp 12.9% noncoding (6/199 nt) PP_mr33 ← Cobalamin regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,765,2180C.12.9% 182.0 / 22.0 62noncoding (6/199 nt)PP_mr33Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base . (4/4);  total (30/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 = 7.64e-01

GCCCGGCTCCGGCCTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACA  >  NC_002947/2765073‑2765361
                                                                                                                                                 |                                                                                                                                                
gCCCGGCTCCGGCCTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATg                                                                                                                                               >  8:46957/1‑149 (MQ=255)
          ggCCTTCCCGGGCGTCAGGCACAGTGGCATCGATGGAGAAGGCTCGGCACCTACAGTTGCGGGGGCAGTTACGATTGGCCTCTATTGGCAGGCATGGGATTCCCTGGTAGTTCCGTGGTTGAAACCGGCGGCCCGAATGGTAGACCATc                                                                                                                                    >  2:283057/1‑149 (MQ=255)
             cTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATcgcg                                                                                                                                 >  6:234710/1‑149 (MQ=255)
                  cGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCGCGATGGTCTACCATGCGggccgc                                                                                                                             >  7:206741/1‑148 (MQ=255)
                  cGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCGCGATGGTCTACCATGCGggccgc                                                                                                                             <  8:206741/148‑1 (MQ=255)
                      cggcAGGCCCAGTGGCATCGATGTAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCAGATTGGCCTGTATTTGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGt                                                                                                                        <  8:5547/149‑1 (MQ=255)
                                aGTGGCATCGTTGGAGCAGGCTCGGCACCTACAGTTCCGGGGGCAGTTCCGATTGGCCTGTATTGGCATGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCtt                                                                                                              <  2:291106/149‑1 (MQ=255)
                                     cATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                         <  7:46957/149‑1 (MQ=255)
                                      aTCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                         >  5:153907/1‑148 (MQ=255)
                                       tCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGaa                                                                                                       <  2:240723/149‑1 (MQ=255)
                                                  ggCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGc                                                                                             >  2:266178/1‑148 (MQ=255)
                                                           ccTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGTAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGt                                                                                   <  5:23033/149‑1 (MQ=255)
                                                            cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     >  1:339983/1‑146 (MQ=255)
                                                            cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     <  2:339983/146‑1 (MQ=255)
                                                            cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTc                                                                                  >  1:207974/1‑149 (MQ=255)
                                                             tACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     <  3:239114/145‑1 (MQ=255)
                                                             tACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     >  4:239114/1‑145 (MQ=255)
                                                                     cGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     <  3:134757/137‑1 (MQ=255)
                                                                     cGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     >  4:134757/1‑137 (MQ=255)
                                                                       ggggCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                         <  7:221650/115‑1 (MQ=255)
                                                                         ggCAGTTCGGATTGGCATGTATTGGCAGGCATGGGATTCCCTGTTATTTCCGTGGTGGCCCACTGCGTCCCGCCCGTTAGAACATCGCGCCGCCAGGTGAATACGCTTTTCAGAAGTACTTCAGCCAACTGATGTCGCGGCGACgctgc                                                                     <  1:283057/149‑1 (MQ=255)
                                                                         ggCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACgctgc                                                                     <  2:207974/149‑1 (MQ=255)
                                                                             gTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCAGCtt                                                                   >  7:297213/1‑147 (MQ=255)
                                                                              ttCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     >  7:251227/1‑128 (MQ=255)
                                                                              ttCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                     <  8:251227/128‑1 (MQ=255)
                                                                                  gATTGGCCTGTATTGGCAGGCCTGGGATTGCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCa                                                                                               <  8:313970/114‑1 (MQ=255)
                                                                                  gATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgg                                                            <  5:234710/149‑1 (MQ=255)
                                                                                  gATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgg                                                            <  4:189257/149‑1 (MQ=255)
                                                                                            tATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATg                                                                                    >  4:43340/1‑115 (MQ=255)
                                                                                            tATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATg                                                                                    <  3:43340/115‑1 (MQ=255)
                                                                                                 gCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTATCTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCAcc                                             <  1:266178/149‑1 (MQ=255)
                                                                                                         gggATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAg                                        >  6:44097/1‑146 (MQ=255)
                                                                                                         gggATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAg                                        <  5:44097/146‑1 (MQ=255)
                                                                                                            aTTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCAGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAg                                  <  8:297213/149‑1 (MQ=255)
                                                                                                             ttCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgc                                 >  7:43093/1‑149 (MQ=255)
                                                                                                                ccTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgcgc                               <  8:43093/148‑1 (MQ=255)
                                                                                                                  tGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAc                            >  3:75048/1‑149 (MQ=255)
                                                                                                                  tGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAc                            >  3:321144/1‑149 (MQ=255)
                                                                                                                    ttAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                           <  4:321144/148‑1 (MQ=255)
                                                                                                                          ccGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcag                      >  5:89551/1‑147 (MQ=255)
                                                                                                                          ccGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcag                      <  6:89551/147‑1 (MQ=255)
                                                                                                                           cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                           <  3:211898/109‑1 (MQ=255)
                                                                                                                           cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                           >  4:211898/1‑109 (MQ=255)
                                                                                                                           cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcagcag                   >  5:230960/1‑149 (MQ=255)
                                                                                                                                 ggAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAgcagca                 <  3:304329/145‑1 (MQ=255)
                                                                                                                                 ggAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACcgc             >  6:37695/1‑149 (MQ=255)
                                                                                                                                  gAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACcgc             <  6:153907/148‑1 (MQ=255)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                               >  8:234698/3‑94 (MQ=25)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                               <  1:121510/92‑1 (MQ=25)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                               <  7:234698/92‑1 (MQ=25)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                               >  2:121510/3‑94 (MQ=25)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                           <  4:266411/128‑1 (MQ=38)
                                                                                                                                      ctggcggcGCGATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                           >  3:266411/3‑130 (MQ=38)
                                                                                                                                       cggcggCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg       >  5:293113/1‑149 (MQ=255)
                                                                                                                                        ggcggcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGc      >  8:145662/1‑149 (MQ=255)
                                                                                                                                         gcggcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGcc          <  4:224077/144‑1 (MQ=255)
                                                                                                                                         gcggcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGcc          >  3:224077/1‑144 (MQ=255)
                                                                                                                                           ggcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                           <  2:142469/93‑1 (MQ=255)
                                                                                                                                           ggcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                           >  1:142469/1‑93 (MQ=255)
                                                                                                                                            gcCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCcaca  >  4:265936/1‑149 (MQ=255)
                                                                                                                                               cGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg       <  4:148536/141‑1 (MQ=255)
                                                                                                                                               cGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg       >  3:148536/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCCCGGCTCCGGCCTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACA  >  NC_002947/2765073‑2765361

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: