Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,284,687 C→T 6.3% N16N (AAC→AAT hbd → 3‑hydroxybutyryl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,284,6870CT 6.3% 235.1 / 4.2 79N16N (AAC→AAThbd3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base T (2/3);  total (36/43)
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.85e-01

CTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCG  >  NC_002947/4284540‑4284832
                                                                                                                                                   |                                                                                                                                                 
cTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACgg                                                                                                                                                 <  1:85021/150‑1 (MQ=255)
         aTCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCa                                                                                                                                        >  2:892718/1‑150 (MQ=255)
          tCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAATGCCGTTGCCCat                                                                                                                                       <  1:254649/150‑2 (MQ=255)
          tCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAg                                                                                                                                       >  2:978079/1‑150 (MQ=255)
           ccAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAATGCCGTTGCCCatg                                                                                                                                      <  2:935476/150‑3 (MQ=255)
           ccAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAgg                                                                                                                                      <  2:1179362/150‑1 (MQ=255)
           ccAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAgg                                                                                                                                      <  2:1048060/150‑1 (MQ=255)
              gggCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGGTGACCCATGAACAGGTTAAAGAACTCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGgtgt                                                                                                                                   >  1:1236904/1‑150 (MQ=255)
               ggCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGtgtg                                                                                                                                  >  1:840157/1‑150 (MQ=255)
                  cGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTgcgc                                                                                                                               <  2:511551/150‑1 (MQ=255)
                  cGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTgcgc                                                                                                                               <  1:95393/150‑1 (MQ=255)
                         cccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGcc                                                                                                                        >  1:789374/1‑150 (MQ=255)
                          ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAATGCCGTTGCCCa                                                                                                                                        <  1:989463/133‑1 (MQ=255)
                          ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAATGCCGTTGCCCa                                                                                                                                        >  2:989463/1‑133 (MQ=255)
                          ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                       >  2:854515/1‑150 (MQ=255)
                          ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                       <  2:42559/150‑1 (MQ=255)
                          ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                       <  2:487261/150‑1 (MQ=255)
                                  cGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAg                                                                                                               <  1:1305308/150‑1 (MQ=255)
                                    aTGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGt                                                                                                             >  2:21719/1‑150 (MQ=255)
                                           tCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGctgctg                                                                                                      <  1:892718/150‑1 (MQ=255)
                                             ggcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggt                                                                                                    <  2:789374/150‑1 (MQ=255)
                                             ggcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggt                                                                                                    <  1:854515/150‑1 (MQ=255)
                                              gcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggtg                                                                                                   <  1:950725/150‑1 (MQ=255)
                                              gcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggtg                                                                                                   >  1:627586/1‑150 (MQ=255)
                                              gcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCtggtg                                                                                                   <  1:978079/150‑1 (MQ=255)
                                                      cAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTc                                                                                           >  1:445795/1‑150 (MQ=255)
                                                      cAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGACCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTc                                                                                           >  1:804083/1‑150 (MQ=255)
                                                       aaGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTcc                                                                                          <  2:1236904/150‑1 (MQ=255)
                                                            aTCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGt                                                                                              >  1:740864/1‑141 (MQ=255)
                                                            aTCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGt                                                                                              <  2:740864/141‑1 (MQ=255)
                                                             tCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                       <  2:653456/115‑1 (MQ=255)
                                                             tCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                       >  1:653456/1‑115 (MQ=255)
                                                               gCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAgc                                                                                  <  2:242192/150‑1 (MQ=255)
                                                               gCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAgc                                                                                  >  2:8803/1‑150 (MQ=255)
                                                                       gtTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAgc                                                                          >  2:175079/1‑150 (MQ=255)
                                                                          cGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGcgcg                                                                       <  1:175079/150‑1 (MQ=255)
                                                                          cGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGcgcg                                                                       <  1:21719/150‑1 (MQ=255)
                                                                           gAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCgg                                                                      >  2:490716/1‑150 (MQ=255)
                                                                           gAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCgg                                                                      <  1:8803/150‑1 (MQ=255)
                                                                            aGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgt                                                                     <  2:627586/150‑1 (MQ=255)
                                                                               gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                                  <  2:942632/150‑1 (MQ=255)
                                                                               gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                                  <  2:1037001/150‑1 (MQ=255)
                                                                                 gTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGc                                                                >  2:886630/1‑150 (MQ=255)
                                                                                     gagggaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCAcc                                                            <  2:212798/150‑1 (MQ=255)
                                                                                       gggaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGAt                                                                                                <  1:549078/112‑1 (MQ=255)
                                                                                       gggaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGAt                                                                                                >  2:549078/1‑112 (MQ=255)
                                                                                       gggaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCtt                                                          <  1:316670/150‑1 (MQ=255)
                                                                                          aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAg                                                       >  2:484411/1‑150 (MQ=255)
                                                                                            gCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTAc                                                                                                                  >  2:340572/1‑89 (MQ=255)
                                                                                            gCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTAc                                                                                                                  <  1:340572/89‑1 (MQ=255)
                                                                                             cTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCaa                                                    <  1:1326118/150‑1 (MQ=255)
                                                                                             cTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCaa                                                    <  1:1326234/150‑1 (MQ=255)
                                                                                                acacATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTg                                                                                                                                            >  1:1080471/1‑59 (MQ=255)
                                                                                                acacATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTg                                                                                                                                            <  2:1080471/59‑1 (MQ=255)
                                                                                                    ctGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTATACGACGCCGTGTTCGATCGTGGGGTTGCCACCTTGAGGAAGTACGTc                                             >  2:1225074/2‑150 (MQ=255)
                                                                                                    aTGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTc                                             <  1:886630/150‑1 (MQ=255)
                                                                                                      gAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGgt                                                                                                                                     >  2:760019/1‑60 (MQ=255)
                                                                                                      gAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGgt                                                                                                                                     <  1:760019/60‑1 (MQ=255)
                                                                                                       aGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAg                                          <  2:445795/150‑1 (MQ=255)
                                                                                                             gAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAg                                    <  2:840157/150‑1 (MQ=255)
                                                                                                             gAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAg                                    >  1:71817/1‑150 (MQ=255)
                                                                                                              aaaaGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGa                                           <  1:949146/142‑1 (MQ=255)
                                                                                                              aaaaGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGa                                           >  2:949146/1‑142 (MQ=255)
                                                                                                                   aTCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCACTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCCCGATCGCGGTGTGGCCACATTGAGCAAGAACCTCGAGCGCCAGGTCAga                              <  2:794520/150‑2 (MQ=255)
                                                                                                                    tCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCa                             >  1:1389889/1‑150 (MQ=255)
                                                                                                                      gCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAg                           <  2:330411/150‑1 (MQ=255)
                                                                                                                       ccGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGAGGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAgg                          >  1:1240808/1‑150 (MQ=255)
                                                                                                                         gTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGc                        <  1:472323/150‑1 (MQ=255)
                                                                                                                             tCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTgg                                                                 >  2:934917/1‑105 (MQ=255)
                                                                                                                             tCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTgg                                                                 <  1:934917/105‑1 (MQ=255)
                                                                                                                                 cgcgGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGa                >  1:1212845/1‑150 (MQ=255)
                                                                                                                                    gggCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGa                                                                            <  1:451330/87‑1 (MQ=255)
                                                                                                                                    gggCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGa                                                                            >  2:451330/1‑87 (MQ=255)
                                                                                                                                       cACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGa          >  2:501013/1‑150 (MQ=255)
                                                                                                                                        aCCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGAc         >  1:704678/1‑150 (MQ=255)
                                                                                                                                        aCCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCAGGCTACAAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGAc         >  1:274062/1‑150 (MQ=255)
                                                                                                                                           ctGGGCAATGCCGTTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAg      >  2:676326/2‑150 (MQ=255)
                                                                                                                                               gCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCg  >  1:47269/1‑150 (MQ=255)
                                                                                                                                               gCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCg  >  1:134457/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGCAGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACCATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTGCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCG  >  NC_002947/4284540‑4284832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: