Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,587,767 A→G 7.4% C32R (TGC→CGC)  PP_4063 ← long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,587,7670AG 7.4% 195.0 / 4.5 67C32R (TGC→CGC) PP_4063long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base G (3/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCGGGAATACTCTG  >  NC_002947/4587619‑4587916
                                                                                                                                                    |                                                                                                                                                     
gcgcAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCaa                                                                                                                                                      >  2:955779/1‑150 (MQ=255)
    aGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGc                                                                                                                                                  <  1:1036056/150‑1 (MQ=255)
              ccAGATGCCAACACGGTCGCCAGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGt                                                                                                                                        >  1:150659/1‑150 (MQ=255)
                 gATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGc                                                                                                                                     >  1:943072/1‑150 (MQ=255)
                      cAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCGGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGa                                                                                                                                >  1:981776/1‑150 (MQ=255)
                       aaCACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGaa                                                                                                                               >  2:409237/1‑150 (MQ=255)
                          acGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGc                                                                                                                            >  2:824342/1‑150 (MQ=255)
                             gTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTg                                                                                                                         <  2:847324/150‑1 (MQ=255)
                                gCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGc                                                                                                                                                  >  2:936942/1‑122 (MQ=255)
                                gCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTCTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGc                                                                                                                                                  <  1:936942/122‑1 (MQ=255)
                                    gtgtTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATg                                                                                                                  <  2:1194842/150‑1 (MQ=255)
                                    gtgtTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATg                                                                                                                  <  1:195757/150‑1 (MQ=255)
                                       ttCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTc                                                                                                               <  1:770158/150‑1 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGc                                                                                                                                              <  1:1169623/116‑1 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGc                                                                                                                                              >  2:1169123/1‑116 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGc                                                                                                                                              >  2:1169623/1‑116 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGc                                                                                                                                              <  1:1169123/116‑1 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                                            >  1:418783/1‑150 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                                            <  2:291387/150‑1 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                                            <  2:943072/150‑1 (MQ=255)
                                          aCGCCCAGTGCGATCAGGGCGCGGGCGTATATTTAGCCCGGGGCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTgg                                                                                                            <  2:286053/150‑1 (MQ=255)
                                           cGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTggg                                                                                                           >  1:1097512/1‑150 (MQ=255)
                                                  tGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCaa                                                                                                    <  2:150659/150‑1 (MQ=255)
                                                   gCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAg                                                                                                   <  2:418783/150‑1 (MQ=255)
                                                      aTCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGggt                                                                                                <  2:150298/150‑1 (MQ=255)
                                                       tCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGggtc                                                                                               >  1:972362/1‑150 (MQ=255)
                                                         aGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGggtcgg                                                                                             >  2:842794/1‑150 (MQ=255)
                                                              gcgGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCATg                                                                                        <  1:34036/150‑1 (MQ=255)
                                                                  gCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAg                                                                                                                              <  1:834042/108‑1 (MQ=255)
                                                                  gCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAg                                                                                                                              >  2:834042/1‑108 (MQ=255)
                                                                       tatTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCa                                                                                                                                    <  2:64853/97‑1 (MQ=255)
                                                                       tatTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCa                                                                                                                                    >  1:64853/1‑97 (MQ=255)
                                                                        atTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCgg                                                                              >  2:1189618/1‑150 (MQ=255)
                                                                        atTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCgg                                                                              <  2:972362/150‑1 (MQ=255)
                                                                          ttCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGt                                                                            >  1:661134/1‑150 (MQ=255)
                                                                          ttCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGAAACAGCCGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGt                                                                            <  1:122926/150‑1 (MQ=255)
                                                                              aCCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACgg                                                                                                                                         >  2:712426/1‑85 (MQ=255)
                                                                              aCCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACgg                                                                                                                                         <  1:712426/85‑1 (MQ=255)
                                                                               cttGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGc                                                                       <  1:131836/148‑1 (MQ=255)
                                                                                  gTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCg                                                                    <  1:955779/150‑1 (MQ=255)
                                                                                     cGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCGCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTTAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGtt                                                                 <  1:1152268/150‑1 (MQ=255)
                                                                                        ccAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTTGCGACAGCgag                                                                                                                                               >  2:1128364/1‑67 (MQ=255)
                                                                                        ccAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTTGCGACAGCgag                                                                                                                                               <  1:1128364/69‑3 (MQ=255)
                                                                                                    aGCTGTACCGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCTGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTcc                                                  <  2:121097/150‑1 (MQ=255)
                                                                                                       tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGTGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTg                                                                                        >  2:573827/1‑109 (MQ=255)
                                                                                                       tGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGTGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTg                                                                                        <  1:573827/109‑1 (MQ=255)
                                                                                                            cGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTtgt                                          <  2:661134/150‑1 (MQ=255)
                                                                                                                  cccTGGTGTCCCGCCACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCcttctt                                            >  2:1031458/1‑142 (MQ=255)
                                                                                                                  cccTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGgttgt                                    >  2:450356/1‑150 (MQ=255)
                                                                                                                       gtggCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTgg                               >  1:1060909/1‑150 (MQ=255)
                                                                                                                            gggACACCAGGGCCTCGCTGTCGCAACAGCGGGCCAGGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACg                          <  2:285384/150‑1 (MQ=255)
                                                                                                                                 acCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACgg                         <  1:835437/146‑1 (MQ=255)
                                                                                                                                 acCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACgg                         >  2:835437/1‑146 (MQ=255)
                                                                                                                                  ccAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTtg                                           >  2:542400/1‑127 (MQ=255)
                                                                                                                                  ccAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTtg                                           <  1:542400/127‑1 (MQ=255)
                                                                                                                                  ccAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGaa                    <  1:51429/150‑1 (MQ=255)
                                                                                                                                  ccAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGaa                    <  1:842794/150‑1 (MQ=255)
                                                                                                                                   cAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTcc                                                  >  1:838090/1‑119 (MQ=255)
                                                                                                                                   cAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTcc                                                  <  2:838090/119‑1 (MQ=255)
                                                                                                                                        ccTCGCTGTCGCAAAAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGaagcaagc              >  2:676663/1‑150 (MQ=255)
                                                                                                                                         cccGCTGTTGCGACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGgttgt                                    >  2:526306/3‑127 (MQ=255)
                                                                                                                                         cccGCTGTTGCGACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGgttgt                                    <  1:526306/125‑1 (MQ=255)
                                                                                                                                         cccGCTGTTGCGACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCg             >  2:711416/3‑150 (MQ=255)
                                                                                                                                          tCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCgg            <  1:824342/150‑1 (MQ=255)
                                                                                                                                          tCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCgg            >  1:349035/1‑150 (MQ=255)
                                                                                                                                               gTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCGGGAATa       <  2:511004/150‑1 (MQ=255)
                                                                                                                                                tCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTg                                >  2:220460/1‑124 (MQ=255)
                                                                                                                                                tCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTg                                <  1:220460/124‑1 (MQ=255)
                                                                                                                                                    aaCAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGACGCAAGCGGGAATACTCTg  <  1:1189618/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCAGTTGGGCGACCAGATGCCAACACGGTCGCCGGTGTTCACGCCCAGTGCGATCAGGGCGCGGGCGTATATTTCGACCTGTTCGGCCAACTGCCGCCAGCTGTAACGCAGGCCCTGGTGGCGGGACACCAGGGCCTCGCTGTCGCAACAGCGGGCCACGGTGGCATCGAAGGCCTGGCCGATGGTCTGGGTCAGCAAGGGTCGGTCCTGGCGACCGCGGGTATAGCTCGGTTGACTCATGGGTGTCCCTTCTTGTGGTTGTTTTGGGCACGGCTGAAGCAAGCGGGAATACTCTG  >  NC_002947/4587619‑4587916

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