Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,556 C→G 8.6% L505L (CTG→CTC PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,5560CG 8.6% 251.0 / 13.3 81L505L (CTG→CTCPP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base G (4/3);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

AGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCT  >  NC_002947/4248415‑4248704
                                                                                                                                              |                                                                                                                                                    
aGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGcc                                                                                                                                              >  2:276258/1‑150 (MQ=255)
aGCCTCCAGGGCGCGGGCATCGCGCAGGCTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGcc                                                                                                                                              >  2:416844/1‑150 (MQ=255)
  cTTCCAGGGCGAGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCg                                                                                                                                            >  2:235989/1‑150 (MQ=255)
    tCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCa                                                                                                                                          <  2:111845/150‑1 (MQ=255)
     ccAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAg                                                                                                                                         <  2:1344996/150‑1 (MQ=255)
      cAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAgg                                                                                                                                        <  1:276258/150‑1 (MQ=255)
       aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                       >  1:1032475/1‑150 (MQ=255)
       aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                       >  2:1242902/1‑150 (MQ=255)
       aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                       <  2:532819/150‑1 (MQ=255)
        gggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTc                                                                                                                                                      <  1:836022/134‑1 (MQ=255)
        gggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTc                                                                                                                                                      >  2:836022/1‑134 (MQ=255)
                 cATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGAc                                                                                                                             <  1:416844/150‑1 (MQ=255)
                     gcgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGCATCAACCGCGGGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCTTCACGCCCGCAAACACCGGCAGACCGaa                                                                                                                         >  1:247440/1‑150 (MQ=255)
                      cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                        <  2:3309/150‑1 (MQ=255)
                      cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                        <  2:801315/150‑1 (MQ=255)
                         aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                        <  1:783047/147‑1 (MQ=255)
                         aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                        >  2:783047/1‑147 (MQ=255)
                           gTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGgcgcg                                                                                                                   >  1:90911/1‑150 (MQ=255)
                            ttCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGTAGACCGAAGgcgcgc                                                                                                                  >  1:414264/1‑150 (MQ=255)
                                ccGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGGCGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGc                                                                                                                                                   <  1:672517/112‑1 (MQ=255)
                                aCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGc                                                                                                                                                   >  2:672517/1‑113 (MQ=255)
                                 cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGc                                                                                                             <  2:601091/150‑1 (MQ=255)
                                          aCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                    <  2:90911/150‑1 (MQ=255)
                                           ccAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                   <  1:60044/150‑1 (MQ=255)
                                                ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGGGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCa                                                                                              >  1:1314479/1‑150 (MQ=255)
                                                ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGTAGGCGCGCAGGGCCGGCTTGGATTGCCa                                                                                              >  2:434693/1‑150 (MQ=255)
                                                ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCa                                                                                              <  1:719542/150‑1 (MQ=255)
                                                   ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTc                                                                                           >  1:142387/1‑150 (MQ=255)
                                                    gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGc                                                                                                >  1:411747/1‑144 (MQ=255)
                                                    gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGc                                                                                                <  2:411747/144‑1 (MQ=255)
                                                       ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                       <  1:1099007/150‑1 (MQ=255)
                                                       ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                       <  1:1099095/150‑1 (MQ=255)
                                                       ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                       >  1:1295780/1‑150 (MQ=255)
                                                         gtgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                     <  2:446985/150‑1 (MQ=255)
                                                              gCAGCTGCGCGCTGACCGCCTTCAGTACAACGGGGTCACCGACGCGGGCGGCATATGACACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACAGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTc                                                                                <  1:945338/150‑1 (MQ=255)
                                                                            aCCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGCGCCTCCTTCCATCGCCATTCGTCCAACACTCCGCTTTCGCCCGTg                                                                  >  2:231894/1‑150 (MQ=255)
                                                                             ccGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGa                                                                                                                              >  2:699207/1‑89 (MQ=255)
                                                                             ccGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGa                                                                                                                              <  1:699207/89‑1 (MQ=255)
                                                                              cGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAg                                                                <  2:1295780/150‑1 (MQ=255)
                                                                              cGCCTTCAGTACCAAGGGGTAACCGACGCAGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAg                                                                >  1:147338/1‑150 (MQ=255)
                                                                               gCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAgg                                                               >  1:177105/1‑150 (MQ=255)
                                                                                ccTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAggg                                                              >  1:269064/1‑150 (MQ=255)
                                                                                   tCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCt                                                           >  1:1139833/1‑150 (MQ=255)
                                                                                      gTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCg                                                        <  1:243455/150‑1 (MQ=255)
                                                                                       tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCAt                                                                                             <  2:764902/112‑1 (MQ=255)
                                                                                       tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCAt                                                                                             >  1:764902/1‑112 (MQ=255)
                                                                                       tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                     <  1:1277904/120‑1 (MQ=255)
                                                                                       tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                     >  2:1277904/1‑120 (MQ=255)
                                                                                       tACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCgg                                                       <  1:1173360/150‑1 (MQ=255)
                                                                                          cAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                    >  1:195568/1‑150 (MQ=255)
                                                                                              gggTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGc                                                <  2:1032475/150‑1 (MQ=255)
                                                                                                     ccGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGAGGGCACAGCGGTAcc                                          >  2:988205/1‑150 (MQ=255)
                                                                                                      cGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAccc                                         <  2:195568/150‑1 (MQ=255)
                                                                                                       gACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCg                                        >  2:536698/1‑150 (MQ=255)
                                                                                                          gcgggcggCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGcag                                     <  2:385399/150‑2 (MQ=255)
                                                                                                             ggcggcGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGc                                  >  1:332572/1‑150 (MQ=255)
                                                                                                                   gTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGcc                            <  2:1139833/150‑1 (MQ=255)
                                                                                                                   ctCTGCCACGGCGCGTTCGGGCGTCCTCAGCAGGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCCCACAGCCTCGcc                            <  2:147338/149‑1 (MQ=255)
                                                                                                                      tGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTg                         <  2:269064/150‑1 (MQ=255)
                                                                                                                      tGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTg                         >  2:1393567/1‑150 (MQ=255)
                                                                                                                      tGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTg                         >  2:509570/1‑150 (MQ=255)
                                                                                                                        ccACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGt                       <  1:584500/150‑1 (MQ=255)
                                                                                                                         cACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTc                      >  2:715080/1‑150 (MQ=255)
                                                                                                                           cGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCa                    >  1:341831/1‑150 (MQ=255)
                                                                                                                            ggCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCAt                   >  1:1322577/1‑150 (MQ=255)
                                                                                                                             gccagtacGGGCGTGCTCTGCACGCCGGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATc                  >  2:1124285/8‑150 (MQ=255)
                                                                                                                               gcgTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTTAAACAGTCGGCTTGCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCaa                <  1:1242902/150‑1 (MQ=255)
                                                                                                                                atctCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGAAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAAc               <  2:970329/147‑1 (MQ=255)
                                                                                                                                 ctgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                        >  1:1075677/4‑140 (MQ=255)
                                                                                                                                 ctgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                        <  2:1075677/137‑1 (MQ=255)
                                                                                                                                  tgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                   >  2:150480/3‑64 (MQ=25)
                                                                                                                                  tgcGGGCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                   <  1:150480/62‑1 (MQ=25)
                                                                                                                                     gggCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                        >  2:65052/1‑136 (MQ=255)
                                                                                                                                     gggCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTgg                        <  1:65052/136‑1 (MQ=255)
                                                                                                                                     gggCGTGCTGAGTACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTg          >  1:965883/1‑150 (MQ=255)
                                                                                                                                     gggCGTACTCAGCACGCCCGCAGGCGCAGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAAccagat          >  1:970329/1‑146 (MQ=255)
                                                                                                                                     gggCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTg          >  1:242353/1‑150 (MQ=255)
                                                                                                                                      ggCGTACTCAGCACGCCCGCAGGCACCGGCAGAACGCTCTCTCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTgg         <  1:434693/150‑1 (MQ=255)
                                                                                                                                          tACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGt     >  1:315175/1‑150 (MQ=255)
                                                                                                                                            cTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTc   >  2:576424/1‑150 (MQ=255)
                                                                                                                                            cTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTc   >  2:638098/1‑150 (MQ=255)
                                                                                                                                             tCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCt  >  1:654681/1‑150 (MQ=255)
                                                                                                                                             tCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCAATCGGTTCt  >  1:60305/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTCT  >  NC_002947/4248415‑4248704

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: