Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,553 A→G 9.9% S506S (AGT→AGC)  PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,5530AG 9.9% 209.4 / 10.3 79S506S (AGT→AGC) PP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base G (5/3);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

TGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGT  >  NC_002947/4248404‑4248701
                                                                                                                                                      |                                                                                                                                                    
tGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTa                                                                                                                                                      <  2:1284490/150‑1 (MQ=255)
           aGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGcc                                                                                                                                           >  2:276258/1‑150 (MQ=255)
           aGCCTCCAGGGCGCGGGCATCGCGCAGGCTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGcc                                                                                                                                           >  2:416844/1‑150 (MQ=255)
             cTTCCAGGGCGAGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCg                                                                                                                                         >  2:235989/1‑150 (MQ=255)
               tCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCa                                                                                                                                       <  2:111845/150‑1 (MQ=255)
                ccAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAg                                                                                                                                      <  2:1344996/150‑1 (MQ=255)
                 cAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAgg                                                                                                                                     <  1:276258/150‑1 (MQ=255)
                  aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                    >  1:1032475/1‑150 (MQ=255)
                  aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                    >  2:1242902/1‑150 (MQ=255)
                  aGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                    <  2:532819/150‑1 (MQ=255)
                   gggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTc                                                                                                                                                   <  1:836022/134‑1 (MQ=255)
                   gggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTc                                                                                                                                                   >  2:836022/1‑134 (MQ=255)
                            cATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGAc                                                                                                                          <  1:416844/150‑1 (MQ=255)
                                gcgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGCATCAACCG‑CGGGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCTTCACGCCCGCAAACACCGGCAGACCGaa                                                                                                                      >  1:247440/1‑150 (MQ=255)
                                 cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                     <  2:3309/150‑1 (MQ=255)
                                 cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                     <  2:801315/150‑1 (MQ=255)
                                    aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                     >  2:783047/1‑147 (MQ=255)
                                    aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAg                                                                                                                     <  1:783047/147‑1 (MQ=255)
                                      gTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGgcgcg                                                                                                                >  1:90911/1‑150 (MQ=255)
                                       ttCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGTAGACCGAAGgcgcgc                                                                                                               >  1:414264/1‑150 (MQ=255)
                                           ccGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGGCGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGc                                                                                                                                                <  1:672517/112‑1 (MQ=255)
                                           aCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGc                                                                                                                                                >  2:672517/1‑113 (MQ=255)
                                            cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGc                                                                                                          <  2:601091/150‑1 (MQ=255)
                                                     aCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGa                                                                                                 <  2:90911/150‑1 (MQ=255)
                                                      ccAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                <  1:60044/150‑1 (MQ=255)
                                                           ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGGGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCa                                                                                           >  1:1314479/1‑150 (MQ=255)
                                                           ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGTAGGCGCGCAGGGCCGGCTTGGATTGCCa                                                                                           >  2:434693/1‑150 (MQ=255)
                                                           ttCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCa                                                                                           <  1:719542/150‑1 (MQ=255)
                                                              ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTc                                                                                        >  1:142387/1‑150 (MQ=255)
                                                               gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGc                                                                                             >  1:411747/1‑144 (MQ=255)
                                                               gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGc                                                                                             <  2:411747/144‑1 (MQ=255)
                                                                  ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                    <  1:1099095/150‑1 (MQ=255)
                                                                  ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                    <  1:1099007/150‑1 (MQ=255)
                                                                  ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                    >  1:1295780/1‑150 (MQ=255)
                                                                    gtgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                  <  2:446985/150‑1 (MQ=255)
                                                                         gCAGCTGCGCGCTGACCGCCTTCAGTACAACGGGGT‑CACCGACGCGGGCGGCATATGACACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACAGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTc                                                                             <  1:945338/150‑1 (MQ=255)
                                                                                       aCCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGCGCCTCCTTCCATCGCCATTCGTCCAACACTCCGCTTTCGCCCGTg                                                               >  2:231894/1‑150 (MQ=255)
                                                                                        ccGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGa                                                                                                                           <  1:699207/89‑1 (MQ=255)
                                                                                        ccGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGa                                                                                                                           >  2:699207/1‑89 (MQ=255)
                                                                                         cGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAg                                                             <  2:1295780/150‑1 (MQ=255)
                                                                                         cGCCTTCAGTACCAAGGGGT‑AACCGACGCAGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAg                                                             >  1:147338/1‑150 (MQ=255)
                                                                                          gCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAgg                                                            >  1:177105/1‑150 (MQ=255)
                                                                                           ccTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAggg                                                           >  1:269064/1‑150 (MQ=255)
                                                                                              tCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCt                                                        >  1:1139833/1‑150 (MQ=255)
                                                                                                 gTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCg                                                     <  1:243455/150‑1 (MQ=255)
                                                                                                  tACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCAt                                                                                          <  2:764902/112‑1 (MQ=255)
                                                                                                  tACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCAt                                                                                          >  1:764902/1‑112 (MQ=255)
                                                                                                  tACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                  >  2:1277904/1‑120 (MQ=255)
                                                                                                  tACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                  <  1:1277904/120‑1 (MQ=255)
                                                                                                  tACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCgg                                                    <  1:1173360/150‑1 (MQ=255)
                                                                                                     cAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                 >  1:195568/1‑150 (MQ=255)
                                                                                                         gggT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGc                                             <  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)
                                                                                                                                                      |                                                                                                                                                    
TGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGT‑AACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGCGTACTCAGCACGCCCGCAGGCACCGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGT  >  NC_002947/4248404‑4248701

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