Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,285,296 T→A 6.3% L234Q (CTG→CAG)  hbd → 3‑hydroxybutyryl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,285,2960TA 6.3% 193.1 / 2.5 63L234Q (CTG→CAG) hbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base A (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACT  >  NC_002947/4285148‑4285442
                                                                                                                                                     |                                                                                                                                                  
ccGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCt                                                                                                                                                    >  3:8678/1‑149 (MQ=255)
 cGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCt                                                                                                                                                    >  5:344526/1‑148 (MQ=255)
                        cTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTgg                                                                                                                             <  4:166023/148‑1 (MQ=255)
                              ccAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCaa                                                                                                                          <  7:190815/145‑1 (MQ=255)
                              ccAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCaa                                                                                                                          >  8:190815/1‑145 (MQ=255)
                               cAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATg                                                                                                                     <  3:23736/149‑1 (MQ=255)
                                        aCGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGATCAGGTCGGccaa                                                                                                                                          >  7:113052/1‑117 (MQ=255)
                                        aCGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGATCAGGTCGGccaa                                                                                                                                          <  8:113052/119‑3 (MQ=255)
                                             gCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGa                                                                                                       <  3:236207/149‑1 (MQ=255)
                                              cGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGaa                                                                                                      <  6:229305/149‑1 (MQ=255)
                                                    tCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAgg                                                                                                                 >  8:346923/1‑132 (MQ=255)
                                                    tCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAgg                                                                                                                 <  7:346923/132‑1 (MQ=255)
                                                     cGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGCCCTCCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCtt                                                                                                 <  1:307790/148‑1 (MQ=255)
                                                     cGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGGCCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTc                                                                                                <  1:169250/148‑1 (MQ=255)
                                                     cGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCtt                                                                                                 <  8:194356/147‑1 (MQ=255)
                                                     cGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTc                                                                                                <  7:23327/148‑1 (MQ=255)
                                                       tGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAAc                                                                                             <  6:344526/149‑1 (MQ=255)
                                                          gcCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGa                                                                                           >  3:50093/1‑148 (MQ=255)
                                                          gcCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGa                                                                                           >  4:362111/1‑148 (MQ=255)
                                                           ccAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGAGACCCTGCTGGCAAGCATGGAGGGCTTCAAGGAAGGCttaatcgc                                                                                           >  6:260424/1‑141 (MQ=255)
                                                                aaGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCaaa                                                                                    >  8:254545/1‑149 (MQ=255)
                                                                  ggCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATa                                                                                  <  6:102242/149‑1 (MQ=255)
                                                                    ccTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATAcc                                                                                <  8:347263/149‑1 (MQ=255)
                                                                     cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATgg                                                                                                                    >  2:128290/1‑112 (MQ=255)
                                                                               ccGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCa                                                                     >  3:255249/1‑149 (MQ=255)
                                                                                       aTCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGTCAGCAAATACCGCCCTGCTCCACTGCTCaa                                                             <  5:149680/149‑1 (MQ=255)
                                                                                       aTCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGcc                                                                             >  6:166785/1‑133 (MQ=255)
                                                                                       aTCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGcc                                                                             <  5:166785/133‑1 (MQ=255)
                                                                                       aTCGACACCGGCATGCGCCTGGGCTGCAACCAACCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTc                                                               >  3:19485/1‑147 (MQ=255)
                                                                                        tCGTCACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGGCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAg                                                            <  3:13329/149‑1 (MQ=255)
                                                                                        tCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCGACGAAGGCTTCCACGACAg                                                                                        >  8:241388/1‑121 (MQ=255)
                                                                                        tCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGaa                                                                                                      <  7:48637/107‑1 (MQ=255)
                                                                                        tCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGaa                                                                                                      >  8:48637/1‑107 (MQ=255)
                                                                                          gACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGc                                                                                       >  5:266610/1‑120 (MQ=255)
                                                                                          gACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGc                                                                                       <  6:266610/120‑1 (MQ=255)
                                                                                              acCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAAt                                                       >  5:60764/1‑149 (MQ=255)
                                                                                                cGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATgg                                                     >  8:185080/1‑149 (MQ=255)
                                                                                                  gCATGCGCCTGGGCTGCAACCAACCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTc                                                   <  4:19485/149‑1 (MQ=255)
                                                                                                        gcCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAGGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGcggcc                                             >  8:289247/1‑149 (MQ=255)
                                                                                                            gggCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCg                                                                               >  8:155784/1‑111 (MQ=255)
                                                                                                            gggCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCg                                                                               <  7:155784/111‑1 (MQ=255)
                                                                                                             taggctaaCCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCggctg                                        <  8:7828/143‑1 (MQ=255)
                                                                                                                   aaCCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGcggc                                              <  3:22684/137‑1 (MQ=255)
                                                                                                                   aaCCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGcggc                                              >  4:22684/1‑137 (MQ=255)
                                                                                                                   aaCCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTggg                                  <  6:60764/149‑1 (MQ=255)
                                                                                                                   aaCCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCGCAAGGAAATGGTCGCGGCCggctgggtcgg                                  >  4:257499/1‑144 (MQ=255)
                                                                                                                     ccAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGgc                                <  4:8678/149‑1 (MQ=255)
                                                                                                                      cAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCa                                                                     <  1:269462/111‑1 (MQ=255)
                                                                                                                      cAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCa                                                                     >  2:269462/1‑111 (MQ=255)
                                                                                                                         ccGAGCGGCCCGCTGGCCTTGGCCGACCTGCTCGGCCTGGACACCCTGCAGGCAATCAAGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCAGCTCCACTGCTCAACGAAATGGTCGCGGCCGGCTGGCTGGGGCGCaa                            <  5:260424/149‑1 (MQ=255)
                                                                                                                                 cccGCTGGCCTTGGCCGACCTGATCTGCCTTGACACCCTGCGGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATTCCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCggctggc                                      <  5:235785/131‑1 (MQ=255)
                                                                                                                                 cccGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTAGCGGCAggctggc                                      >  6:235785/1‑131 (MQ=255)
                                                                                                                                  ccGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCggctgg                                       >  7:267695/1‑129 (MQ=255)
                                                                                                                                  ccGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCggctgg                                       <  8:267695/129‑1 (MQ=255)
                                                                                                                                   cGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTcgc                  >  1:346852/1‑149 (MQ=255)
                                                                                                                                    gCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTcgcg                 >  8:66497/1‑149 (MQ=255)
                                                                                                                                     cTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCgg                >  7:128764/1‑149 (MQ=255)
                                                                                                                                        gCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGtt              >  8:117719/1‑148 (MQ=255)
                                                                                                                                        gCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGtt              <  3:362111/148‑1 (MQ=255)
                                                                                                                                        gCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGt               <  7:289247/147‑1 (MQ=255)
                                                                                                                                          cTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTc                    <  3:58057/140‑1 (MQ=255)
                                                                                                                                          cTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTc                    >  4:58057/1‑140 (MQ=255)
                                                                                                                                            aggCCGATCAGGTCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttcttc         <  3:98020/148‑1 (MQ=255)
                                                                                                                                            aggCCGATCAGGTCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTttctt          >  4:98020/2‑148 (MQ=255)
                                                                                                                                                   ccTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGTGGTCGCGGTTTCTTCACCTACt  >  4:116627/1‑149 (MQ=255)
                                                                                                                                                     tGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACt  >  1:39302/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACT  >  NC_002947/4285148‑4285442

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: