Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,285,233 A→T 5.8% D213V (GAC→GTC)  hbd → 3‑hydroxybutyryl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,285,2330AT 5.8% 207.2 / 2.3 69D213V (GAC→GTC) hbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TACGCCCAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCC  >  NC_002947/4285085‑4285367
                                                                                                                                                    |                                                                                                                                      
tACGCCCAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGa                                                                                                                                        >  7:37267/1‑149 (MQ=255)
  cGCCCAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCGCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGAc                                                                                                                                       <  5:193332/148‑1 (MQ=255)
     ccAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGTCATCg                                                                                                                                   <  6:238573/149‑1 (MQ=255)
     ccAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATACTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATc                                                                                                                                    >  3:258196/1‑148 (MQ=255)
        ggCGCGGCTGGTGACCGAAAAAGTCGTCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGaca                                                                                                                                <  4:258196/149‑1 (MQ=255)
           gcTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCg                                                                                                                             >  8:265740/1‑149 (MQ=255)
           gcTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCg                                                                                                                             <  7:282547/149‑1 (MQ=255)
               aTGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGCTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGTCCAGTGCCGAGGACATCGACACCCGCAt                                                                                                                         >  8:181047/1‑149 (MQ=255)
                   tGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATgcgc                                                                                                                     >  6:221336/1‑149 (MQ=255)
                      ccGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCt                                                                                                                   <  4:282608/148‑1 (MQ=255)
                        gAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCg                                                                                                                                   <  8:102492/130‑1 (MQ=255)
                        gAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCg                                                                                                                                   >  7:102492/1‑130 (MQ=255)
                          aaaaGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCt                                                                                                              <  1:330280/149‑1 (MQ=255)
                                            gATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGcc                                                                                            >  4:263843/1‑149 (MQ=255)
                                                aCCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCt                                                                                        >  2:71651/1‑149 (MQ=255)
                                                aCCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCt                                                                                        >  3:263696/1‑149 (MQ=255)
                                                  cgccgGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTgg                                                                                      >  1:17382/1‑149 (MQ=255)
                                                        cAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTgg                                                                                <  8:37267/149‑1 (MQ=255)
                                                         aaCCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGc                                                                               >  7:100589/1‑149 (MQ=255)
                                                         aaCCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGc                                                                               >  6:37748/1‑149 (MQ=255)
                                                            cGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGa                                                                            >  7:108652/1‑149 (MQ=255)
                                                             gCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCTATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTgg                                                                                <  7:365699/144‑1 (MQ=255)
                                                             gCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGAc                                                                           <  3:210147/149‑1 (MQ=255)
                                                             gCCCGGGCTTCCTGGTCAACCGCATCCTGGTGCCAATGATAAACGACGCGATCATCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGCCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCCCTTGCCTTgg                                                                                >  8:365699/1‑144 (MQ=255)
                                                              cccGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGAcc                                                                          <  2:39824/149‑1 (MQ=255)
                                                              cccGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGAcc                                                                          <  5:214215/149‑1 (MQ=255)
                                                                    cTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCg                                                                    <  3:140169/149‑1 (MQ=255)
                                                                         tgctCAACCGCATCCTGGTGCCAATGATCAACGAGGCTATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCATCCGATCGGCCCGCTGGCCTTGGCCGATCAGGTCGGccaa                                                               <  4:233252/146‑3 (MQ=255)
                                                                           gTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGcc                                                                                            <  7:37388/118‑1 (MQ=255)
                                                                           gTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGcc                                                                                            >  8:37388/1‑118 (MQ=255)
                                                                            tCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGa                                                             >  4:326781/1‑148 (MQ=255)
                                                                              aaCCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGacac                                                          >  2:313627/1‑149 (MQ=255)
                                                                               aCCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACAcc                                                         <  6:127759/149‑1 (MQ=255)
                                                                                    aTCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCctgct                                                    <  4:263696/149‑1 (MQ=255)
                                                                                       ttGGTGCCAATGATCAACGAGGCGATCTTCGCGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGc                                                 <  3:70018/148‑1 (MQ=255)
                                                                                       cTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTgg                                                  <  5:221336/148‑1 (MQ=255)
                                                                                       cTGGTGCCAATGATCAACGAGGCGATCTTCATGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTgg                                                  >  2:43162/1‑148 (MQ=255)
                                                                                         ggTGCCAATGATCAACGAGGCGATCTTCGCGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGc                                                 >  4:70018/1‑147 (MQ=255)
                                                                                          gTGCCATTGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCTAGTACATCGACACCGGCATGCGCCTGGGCTGCTACAAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAAt                                              <  7:181047/149‑1 (MQ=255)
                                                                                          gTGCCAATGGTCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACACTGCTGGCAAt                                              >  8:345673/1‑149 (MQ=255)
                                                                                           tGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATc                                             <  5:37748/149‑1 (MQ=255)
                                                                                              cAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATg                                          <  2:184435/149‑1 (MQ=255)
                                                                                               aaTGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATg                                          >  5:142591/1‑148 (MQ=255)
                                                                                                aTGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGacac                                                          <  7:15639/131‑1 (MQ=255)
                                                                                                aTGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGacac                                                          >  8:15639/1‑131 (MQ=255)
                                                                                                aTGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGa                                        <  6:176854/149‑1 (MQ=255)
                                                                                                 tGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAg                                       <  3:161571/149‑1 (MQ=255)
                                                                                                  gATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAgg                                      >  2:376576/1‑149 (MQ=255)
                                                                                                    tCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGcc                                    >  3:322402/1‑149 (MQ=255)
                                                                                                      aaCGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCtt                                  >  6:30414/1‑149 (MQ=255)
                                                                                                      aaCGAGGCGATCTTCATGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCtt                                  <  1:43162/149‑1 (MQ=255)
                                                                                                                    cGTGCGCCAGGAAGGCCTTGCCTGGGCCGAGGACATCGACACCGGCTTGCGCCGGGGCTGCAACCAGCCGATAGGCCCGCTGGCCTTGGCCGACCTGATCGGGCTGGACACCCTGCTTGCAATCATGGAGGCCTTCCACGAAGGCTTc                     <  6:103192/148‑1 (MQ=255)
                                                                                                                        cgcCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGa                <  1:376576/149‑1 (MQ=255)
                                                                                                                            aGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGc            <  7:265740/149‑1 (MQ=255)
                                                                                                                                 ggCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATa       >  1:244013/1‑149 (MQ=255)
                                                                                                                                 gcgcatgCCGGTGTCGATGTCCTCGGCACTGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATa       >  8:359142/7‑149 (MQ=255)
                                                                                                                                  gCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGa                                        <  6:381392/115‑1 (MQ=255)
                                                                                                                                  gCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGa                                        >  5:381392/1‑115 (MQ=255)
                                                                                                                                  gCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATAc      <  2:244013/149‑1 (MQ=255)
                                                                                                                                   ccTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACACTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATAcc     >  8:325306/1‑149 (MQ=255)
                                                                                                                                    cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGTCTTCCAc                              >  4:117809/1‑123 (MQ=255)
                                                                                                                                    cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGTCTTCCAc                              <  3:117809/123‑1 (MQ=255)
                                                                                                                                    cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAAc                  >  4:70068/1‑135 (MQ=255)
                                                                                                                                    cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAAc                  <  3:70068/135‑1 (MQ=255)
                                                                                                                                    cTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCg    >  6:349474/1‑149 (MQ=255)
                                                                                                                                      tgCCTGTGCCGAGGACATCGACACCGGCATGCGGCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCGACGACAGCAAATACCGcc  <  6:260438/148‑1 (MQ=255)
                                                                                                                                       gCCGGTGTCGATGTCCTCGGCACTGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATa       <  2:231272/143‑1 (MQ=255)
                                                                                                                                       gCCGGTGTCGATGTCCTCGGCACTGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATa       >  1:231272/1‑143 (MQ=255)
                                                                                                                                          aGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTcc                                <  6:145897/115‑1 (MQ=255)
                                                                                                                                          aGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTcc                                >  5:145897/1‑115 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
TACGCCCAGGCGCTGATGGTGACCGAAAAAGTCGGCAAGACCCCGATCACCGCCGGCAACCGCCCGGGCTTCGTGGTCAACCGCATCCTGGTGCCAATGATCAACGAGGCGATCTTCGTGCGCCAGGAAGGCCTGGCCAGTGCCGAGGACATCGACACCGGCATGCGCCTGGGCTGCAACCAGCCGATCGGCCCGCTGGCCTTGGCCGACCTGATCGGCCTGGACACCCTGCTGGCAATCATGGAGGCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCC  >  NC_002947/4285085‑4285367

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