Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,318,705 C→T 5.3% R302C (CGC→TGC)  PP_2039 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,318,7050CT 5.3% 226.2 / 2.5 78R302C (CGC→TGC) PP_2039acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

CCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGT  >  NC_002947/2318560‑2318854
                                                                                                                                                 |                                                                                                                                                     
ccccACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACgcag                                                                                                                                                   >  1:8387/1‑150 (MQ=255)
   cACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCAt                                                                                                                                                >  2:187699/1‑150 (MQ=255)
    aCCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATg                                                                                                                                               <  2:8387/150‑1 (MQ=255)
    aCCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATg                                                                                                                                               <  2:319317/150‑1 (MQ=255)
    aCCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATg                                                                                                                                               <  1:2234300/150‑1 (MQ=255)
       ttGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCGTGCAGCATgcag                                                                                                                                            <  1:1103923/150‑4 (MQ=255)
       ttGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTc                                                                                                                                            >  1:2397064/1‑150 (MQ=255)
              ggCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCa                                                                                                                                     <  2:2397064/150‑1 (MQ=255)
                 tGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACgg                                                                                                                                  >  1:1535111/1‑150 (MQ=255)
                   gccggccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCa                                                                                                                                >  1:2299886/1‑150 (MQ=255)
                      ggccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCcgcg                                                                                                                             >  1:959960/1‑150 (MQ=255)
                        ccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCAc                                                                                                                                    >  2:280972/1‑141 (MQ=255)
                        ccTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCAc                                                                                                                                    <  1:280972/141‑1 (MQ=255)
                          tGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGTTGCGTGCAGCATgc                                                                                                                                              >  1:2164808/1‑128 (MQ=255)
                          tGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGTTGCGTGCAGCATgc                                                                                                                                              <  2:2164808/129‑2 (MQ=255)
                           ggCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGaa                                                                                                                        <  2:1535111/150‑1 (MQ=255)
                            gCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       <  2:959960/150‑1 (MQ=255)
                            gCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       <  2:89982/150‑1 (MQ=255)
                                tGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACg                                                                                                                   >  1:2296535/1‑150 (MQ=255)
                                    ccTCTCGATCGACTATGCCGGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAATCGATCAACCATATGATCGCCGACATGGCCACCGATCTGCATGCTGCACGCAGCATGTTCTACTCCGCGGCAGCGGTGAATGACGCcgg                                                                                                               >  1:859073/1‑150 (MQ=255)
                                      tctcGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCcggcg                                                                                                             >  1:1905212/1‑150 (MQ=255)
                                        tcgatcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCggcggc                                                                                                           >  2:1297272/1‑150 (MQ=255)
                                          gatcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGTACgcagc                                                                                                                                                  <  2:1735963/109‑1 (MQ=255)
                                          gatcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGTACgcagc                                                                                                                                                  >  1:1735963/1‑109 (MQ=255)
                                          gatcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                         <  1:1341942/134‑1 (MQ=255)
                                          gatcgaCTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                         >  2:1341942/1‑134 (MQ=255)
                                                cTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCg                                                                                                   >  2:1181729/1‑150 (MQ=255)
                                                 tATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGc                                                                                                  <  2:2296535/150‑1 (MQ=255)
                                                 tATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCAGTTTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGc                                                                                                  <  1:1833084/150‑1 (MQ=255)
                                                      ccGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCa                                                                                             >  1:1332693/1‑150 (MQ=255)
                                                       cGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAg                                                                                            <  2:888546/150‑1 (MQ=255)
                                                       cGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCAAGAAGACCCCCTGTAGAATCTCGGACATTGGTACCGTGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGGCGCCGGCGGCAATATCCGCACCCAg                                                                                            >  1:45183/1‑150 (MQ=255)
                                                              gcAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGa                                                                                     <  2:2313471/150‑1 (MQ=255)
                                                                aaGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGAtg                                                                                   >  2:1303919/1‑150 (MQ=255)
                                                                     gTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCaa                                                                              <  2:1076609/150‑1 (MQ=255)
                                                                      ttCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       >  1:336205/1‑108 (MQ=255)
                                                                      ttCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       <  2:336205/108‑1 (MQ=255)
                                                                      ttCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAg                                                                             >  1:792995/1‑150 (MQ=255)
                                                                      ttCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTCCATGCTGCAAGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAg                                                                             >  1:1962517/1‑150 (MQ=255)
                                                                      ttCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATAGGCACAGGGTGGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAg                                                                             >  2:1532304/1‑150 (MQ=255)
                                                                           ccAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGtt                                                                        <  2:780810/150‑1 (MQ=255)
                                                                            cAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTc                                                                       <  1:1167382/150‑1 (MQ=255)
                                                                            cAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTc                                                                       <  2:1905212/150‑1 (MQ=255)
                                                                                      ccATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGgaga                                                             >  1:1513885/1‑150 (MQ=255)
                                                                                        aTGTTCCAGTCGATCAACCATATGATCGCCGACATTGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACg                                                           <  2:1027584/150‑1 (MQ=255)
                                                                                               aGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCaa                                                                              >  2:2322948/1‑124 (MQ=255)
                                                                                               aGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCaa                                                                              <  1:2322948/124‑1 (MQ=255)
                                                                                                gTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       <  1:922037/82‑1 (MQ=255)
                                                                                                gTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAAt                                                                                                                       >  2:922037/1‑82 (MQ=255)
                                                                                                 tCGATCAACCATATGATCGCCGACCCGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTcc                                                    <  1:1220/148‑1 (MQ=255)
                                                                                                 tCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTcc                                                    >  2:1220/1‑148 (MQ=255)
                                                                                                    aTCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATc                                               <  1:187699/150‑1 (MQ=255)
                                                                                                       aaCCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGcc                                            <  2:1513885/150‑1 (MQ=255)
                                                                                                         ccATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGa                                          >  1:1096000/1‑150 (MQ=255)
                                                                                                            tatGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATcg                                       >  1:2163040/1‑150 (MQ=255)
                                                                                                                aTCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGGACCCGGGCCCCGATGTGTAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGcc                                   <  1:1472771/150‑1 (MQ=255)
                                                                                                                    ccGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTAc                                                                                                                                         >  2:1014915/1‑44 (MQ=255)
                                                                                                                    ccGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTAc                                                                                                                                         <  1:1014915/44‑1 (MQ=255)
                                                                                                                    ccGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGc                               >  2:335818/1‑150 (MQ=255)
                                                                                                                     cggcATGGGCACCGAGCTGCATGCTGCAGGCAGCATGGTATACGCCACGGCAGCGGTGAATGACGCGGGGGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCa                              <  2:279586/147‑1 (MQ=255)
                                                                                                                     cGACATGGCCACCGAGCTGCATGCTGCACGCAGCATTGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCTTTTCGTAGACGTCCTTCCGCATCACCGATCGCGCCGTGCa                              <  2:2299886/150‑1 (MQ=255)
                                                                                                                       aCATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGa                            <  1:335818/150‑1 (MQ=255)
                                                                                                                       aCATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGa                            <  1:1671732/150‑1 (MQ=255)
                                                                                                                           ggCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTc                                                                  <  1:1506648/108‑1 (MQ=255)
                                                                                                                           ggCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTc                                                                  >  2:1506648/1‑108 (MQ=255)
                                                                                                                           ggCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTc                                                     >  2:1689869/1‑121 (MQ=255)
                                                                                                                           ggCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTc                                                     <  1:1689869/121‑1 (MQ=255)
                                                                                                                                ccGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggcg                   <  1:683084/150‑1 (MQ=255)
                                                                                                                                 cGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACAcggcgg                  >  1:396359/1‑150 (MQ=255)
                                                                                                                                      tGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGa                                                                                     >  2:250654/1‑78 (MQ=255)
                                                                                                                                      tGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGa                                                                                     <  1:250654/78‑1 (MQ=255)
                                                                                                                                       gCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCgg            <  1:1303919/150‑1 (MQ=255)
                                                                                                                                        cATGCTGCGTGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCggg           >  2:1259410/1‑150 (MQ=255)
                                                                                                                                         aTGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGACAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGgc          >  2:282385/1‑150 (MQ=255)
                                                                                                                                           gctgcACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                         >  1:667199/1‑37 (MQ=255)
                                                                                                                                           gctgcACGCAGCATGGTCTACGCCACGGCAGCGGTGa                                                                                                                         <  2:667199/37‑1 (MQ=255)
                                                                                                                                               cACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCtgctgc    >  2:1896982/1‑150 (MQ=255)
                                                                                                                                                aCGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCGGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCg   >  1:2106605/1‑150 (MQ=255)
                                                                                                                                                 cGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGt  <  2:1332693/150‑1 (MQ=255)
                                                                                                                                                 cGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGt  >  1:483981/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCCCACCTTGCTGGGGCTGGCCGGCCTGGCACTGAACCTCTCGATCGACTATGCCCGCACGCGCAAGCAGTTCGGCCAGCCTATCGCCATGTTCCAGTCGATCAACCATATGATCGCCGACATGGCCACCGAGCTGCATGCTGCACGCAGCATGGTCTACGCCACGGCAGCGGTGAATGACGCCGGCGGCAATATCCGCACCCAGGCCCCGATGTGCAAGCTGTTCGTTTCGGAGACGGCCTTCCGCATCGCCGATCGCGCCGTGCAGATACACGGCGGTGCCGGGCTGCTGCGT  >  NC_002947/2318560‑2318854

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