Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,844,929 C→G 7.9% A235A (GCG→GCC PP_3394 ← 3‑hydroxy‑3‑methylglutaryl‑CoA lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,844,9290CG 7.9% 187.2 / 7.0 63A235A (GCG→GCCPP_33943‑hydroxy‑3‑methylglutaryl‑CoA lyase
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base G (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

CGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGG  >  NC_002947/3844785‑3845076
                                                                                                                                                  |                                                                                                                                                     
cgcttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGc                                                                                                                                                     >  6:55855/1‑147 (MQ=255)
  cttcgcttGGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAgg                                                                                                                                                 <  4:77123/149‑1 (MQ=255)
  cttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAgg                                                                                                                                                 <  6:124906/149‑1 (MQ=255)
  cttcgcttAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAgg                                                                                                                                                 <  5:55855/149‑1 (MQ=255)
     cgcttAGTACCTGCCGGGTTTCGATTAGCTCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCTGGCGGGCCTGGGa                                                                                                                                               <  1:227917/148‑1 (MQ=255)
          aGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATc                                                                                                                                         >  1:64912/1‑149 (MQ=255)
             aCCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAg                                                                                                                                      >  6:274885/1‑149 (MQ=255)
                  ccGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACAGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGtt                                                                                                                                    <  1:64594/147‑1 (MQ=255)
                  ccGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTccc                                                                                                                                 <  1:303401/149‑1 (MQ=255)
                          cGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGcc                                                                                                                          <  4:205920/148‑1 (MQ=255)
                                gcgcAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCctg‑‑ctg                                                                                                                  <  5:274885/148‑1 (MQ=255)
                                      gCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGACCTG‑‑CTGCAACGc                                                                                                            >  6:98770/1‑148 (MQ=255)
                                        aaGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGccac                                                                                                         >  1:166160/1‑149 (MQ=255)
                                        aaGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGccac                                                                                                         >  8:111228/1‑149 (MQ=255)
                                          ggcGATTCCGGTGGCGTAGCCCATGCTTTGGAGCATGAACACGAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCCGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTT‑‑CTGCACCGccacca                                                                                                       <  4:262949/147‑1 (MQ=255)
                                             gATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAgg                                                                                                     <  2:166160/148‑1 (MQ=255)
                                                  cGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGcc                                                                                                           <  1:49833/137‑1 (MQ=255)
                                                  cGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGcc                                                                                                           >  2:49833/1‑137 (MQ=255)
                                                         gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGCGGCAGCCTCCCAGGCCCGCCAGGGAGGCATc                                                                                                                                         <  1:37764/103‑1 (MQ=255)
                                                         gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATc                                                                                                                                         >  2:37764/1‑102 (MQ=255)
                                                         gTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAag                                                                                         <  1:33497/148‑1 (MQ=255)
                                                           aGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAaggc                                                                                       >  2:286196/1‑148 (MQ=255)
                                                            gCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAaggc                                                                                       >  4:227257/1‑147 (MQ=255)
                                                              ccATGCTTTGCAGCATGAACACGGGGTACTCGGTCAGGTTATAACCAGACGCGCCAGGGGCGAACGGGCACCCACCCCGGCCCGCCAGGGAGGCAGGCAGGGCCACAAGGCGGAG‑‑CAGCACCGCCCCCAGGGTGGTGGCCAaggccggg                                                                                   >  7:54792/1‑149 (MQ=255)
                                                              ccATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAaggccagg                                                                                   >  1:233139/1‑149 (MQ=255)
                                                               cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAaggccaggc                                                                                  >  4:267307/1‑149 (MQ=255)
                                                               cATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCACGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCACTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTCTTGGCCAaggccaggc                                                                                  >  6:183784/1‑149 (MQ=255)
                                                                      tGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCgg                                                                             <  1:156935/147‑1 (MQ=255)
                                                                         agcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTc                                                                        <  2:64912/149‑1 (MQ=255)
                                                                         agcaTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTc                                                                        <  2:233139/149‑1 (MQ=255)
                                                                           caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGgt                                                                           >  6:289591/1‑144 (MQ=255)
                                                                           caTGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGgt                                                                           <  5:289591/144‑1 (MQ=255)
                                                                             tGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATg                                                                     >  2:353648/1‑148 (MQ=255)
                                                                                  acCAGGTCCTCGGTCACGGTATTACAAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTg                                                               >  1:129450/1‑149 (MQ=255)
                                                                                   ccAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTc                                                                        >  4:105842/1‑139 (MQ=255)
                                                                                   ccAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTc                                                                        <  3:105842/139‑1 (MQ=255)
                                                                                    cAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTgcg                                                             >  5:306738/1‑149 (MQ=255)
                                                                                          cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCt                                                       <  3:227257/149‑1 (MQ=255)
                                                                                          cTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCt                                                       >  7:91598/1‑149 (MQ=255)
                                                                                           tCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCtt                                                      >  5:337052/1‑149 (MQ=255)
                                                                                                  cGGTATTACCAGACGCGCCAGGGGCGAACGGTCAGCCTCCCAGGCCCGCCATGGAGGCATCCATTTCCCCAATGCCCTGCTCTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGAc                                                 <  1:286196/149‑1 (MQ=255)
                                                                                                   ggTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGt                                              >  8:61134/1‑149 (MQ=255)
                                                                                                           cAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCa                                      <  7:111228/149‑1 (MQ=255)
                                                                                                            aGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCaa                                     <  1:353648/149‑1 (MQ=255)
                                                                                                              aCGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtc                                   <  7:61134/149‑1 (MQ=255)
                                                                                                                gcgcCAGGGGCGAAGGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTGCCCAATGCCCTG‑‑CTGCACCGCGACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCTAtctc                                 <  5:183784/149‑1 (MQ=255)
                                                                                                                gcgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtctc                                 <  8:91598/149‑1 (MQ=255)
                                                                                                                gcgcCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAAtctc                                 <  3:267307/149‑1 (MQ=255)
                                                                                                                    cAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCCCCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCg                              >  7:90308/1‑148 (MQ=255)
                                                                                                                     aGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTa                            >  6:243885/1‑149 (MQ=255)
                                                                                                                     aGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTa                            >  6:279682/1‑149 (MQ=255)
                                                                                                                                 tggatgCCTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGc                                                         <  5:346032/105‑1 (MQ=38)
                                                                                                                                 tggatgCCTCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGc                                                         >  6:346032/6‑110 (MQ=38)
                                                                                                                                  ggCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGc                 <  6:306738/149‑1 (MQ=255)
                                                                                                                                    cAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCGGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCAc               <  8:54792/149‑1 (MQ=255)
                                                                                                                                        cTCCCTGGCGGGCCTGGGAGGCTGCCCGTTCGCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCAtt           >  6:168232/1‑149 (MQ=255)
                                                                                                                                         tCCCTGGCGGGCCTGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg          >  3:186156/1‑149 (MQ=255)
                                                                                                                                         tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTc                                             <  6:201142/114‑1 (MQ=255)
                                                                                                                                         tCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAAAGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGAGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTg          >  5:259131/1‑149 (MQ=255)
                                                                                                                                          cccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGc         >  8:38973/1‑149 (MQ=255)
                                                                                                                                           ccAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGcc        <  5:240283/149‑1 (MQ=255)
                                                                                                                                              ggCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCAcc     <  5:279682/149‑1 (MQ=255)
                                                                                                                                               gCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCt    >  5:38655/1‑149 (MQ=255)
                                                                                                                                                cccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTgaga                                                  >  2:326518/1‑102 (MQ=255)
                                                                                                                                                cccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTgaga                                                  <  1:326518/102‑1 (MQ=255)
                                                                                                                                                 ccGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTgg  <  3:136874/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGCTTCGCTTAGTACCTGCCGGGTTTCGATTAGCGCAGGCAAGTCGATGCCGGTGTCGTAGCCCATGCTTTGCAGCATGAACACCAGGTCCTCGGTCACGGTATTACCAGACGCGCCAGGGGCGAACGGGCAGCCTCCCAGGCCCGCCAGGGAGGCATCCAGTTCCCCAATGCCCTG‑‑CTGCACCGCCACCAGGCTGTTGGCCAAGGCCAGGCCCCGGGTGTCATGGAAGTGCGCGGCCTTGAGACGGTCGCCTGCAATCTCCCGTACGGCCTGAAGCACCATTGCCACCTGG  >  NC_002947/3844785‑3845076

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: