Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,717,206 G→C 5.4% A626A (GCG→GCC cti → esterified fatty acid cis/trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,717,2060GC 5.4% 295.5 / 6.5 92A626A (GCG→GCCctiesterified fatty acid cis/trans isomerase
Reads supporting (aligned to +/- strand):  ref base G (38/49);  new base C (2/3);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GACTGAAGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTACC  >  NC_002947/2717058‑2717345
                                                                                                                                                    |                                                                                                                                           
gACTGAAGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCcggcgg                                                                                                                                            >  1:1890362/1‑150 (MQ=255)
gACTGAAGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCcggcgg                                                                                                                                            >  2:1854411/1‑150 (MQ=255)
   tGAAGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCg                                                                                                                                         >  2:507000/1‑150 (MQ=255)
      aGCTGCTGCAGCGCACCGGGAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAAAATTCCGCAATGCCGATCAGGCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGc                                                                                                                                      <  2:711324/150‑1 (MQ=255)
      aGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCATGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGc                                                                                                                                      <  1:1531958/150‑1 (MQ=255)
       gctgctGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCt                                                                                                                                     >  2:1650728/1‑150 (MQ=255)
          gctgcAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAAAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGaa                                                                                                                                  >  1:1772511/1‑150 (MQ=255)
           ctgcAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCAGCCCGGCCGCCGGGcgtgaca                                                                                                                                <  2:1839449/150‑7 (MQ=255)
               aGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGa                                                                                                                             <  2:810001/150‑1 (MQ=255)
                gcgcACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGAt                                                                                                                            <  1:1866251/150‑1 (MQ=255)
                  gcACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCa                                                                                                                          >  1:129936/1‑150 (MQ=255)
                    aCCGGCAGCTGTAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        <  1:574476/150‑1 (MQ=255)
                    aCCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        <  2:1009247/150‑1 (MQ=255)
                    aCCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        <  1:1463343/150‑1 (MQ=255)
                      cGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTGAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAccc                                                                                                                      >  1:604730/1‑149 (MQ=255)
                      cGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCa                                                                                                                      <  1:1650728/150‑1 (MQ=255)
                        gCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        >  2:585378/1‑146 (MQ=255)
                        gCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        <  1:585378/146‑1 (MQ=255)
                         cAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGtt                                                                                                                   <  1:1854411/150‑1 (MQ=255)
                         cAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCGGAAGGTGATCAACCAGtt                                                                                                                   >  2:308082/1‑150 (MQ=255)
                              tGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCg                                                                                                              <  1:807825/150‑1 (MQ=255)
                                aaTGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGa                                                                                                                             >  2:1900255/1‑133 (MQ=255)
                                aaTGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGa                                                                                                                             <  1:1900255/133‑1 (MQ=255)
                                aaTGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAg                                                                                                            >  2:1208155/1‑150 (MQ=255)
                                    ccGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGa                                                                                                        <  1:406198/150‑1 (MQ=255)
                                          cGGCCCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGTCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAtgc                                                                                                  <  1:1375049/150‑1 (MQ=255)
                                          cGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        <  1:716490/128‑1 (MQ=255)
                                          cGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAAc                                                                                                                        >  2:716490/1‑128 (MQ=255)
                                          cGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAtgc                                                                                                  <  1:1843588/150‑1 (MQ=255)
                                            gACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAtgctg                                                                                                >  2:175074/1‑150 (MQ=255)
                                             aCCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAtgctgc                                                                                               <  1:1534182/150‑1 (MQ=255)
                                             aCCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAtgctgc                                                                                               <  2:1585646/150‑1 (MQ=255)
                                              cccGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCg                                                                                              <  2:984070/150‑1 (MQ=255)
                                              cccGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCCGCCGGGcgtga                                                                                                                                  <  2:960972/114‑5 (MQ=255)
                                              cccGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCCGCCGGGcgtga                                                                                                                                  >  1:960972/1‑110 (MQ=255)
                                               ccGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGt                                                                                             >  2:1701809/1‑150 (MQ=255)
                                                 gATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCa                                                                                                                      >  2:1726225/1‑123 (MQ=255)
                                                 gATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCa                                                                                                                      <  1:1726225/123‑1 (MQ=255)
                                                     aaCCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGaa                                                                                       <  2:1067859/150‑1 (MQ=255)
                                                      aCCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAg                                                                                      <  1:540345/150‑1 (MQ=255)
                                                          cTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCa                                                                                  >  1:1833994/1‑150 (MQ=255)
                                                             ccAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGa                                                                                                                             >  2:173275/1‑104 (MQ=255)
                                                             ccAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGa                                                                                                                             <  1:173275/104‑1 (MQ=255)
                                                             ccAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAt                                                                                                    >  1:289097/1‑129 (MQ=255)
                                                             ccAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCAt                                                                                                    <  2:289097/129‑1 (MQ=255)
                                                             ccAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGa                                                                               >  2:1113436/1‑150 (MQ=255)
                                                                gggCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACgg                                                                            >  1:975149/1‑150 (MQ=255)
                                                                 ggCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGc                                                                           >  2:1012246/1‑150 (MQ=255)
                                                                 ggCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGc                                                                           >  2:1687277/1‑150 (MQ=255)
                                                                  gcgcGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGcc                                                                          >  1:1392073/1‑150 (MQ=255)
                                                                   cgcgTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCa                                                                         <  1:1701809/150‑1 (MQ=255)
                                                                       ttCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACTATGCTGCGTATTGAAGGGCTGGACGGCGAGCGt                                                                     >  2:837141/1‑150 (MQ=255)
                                                                       ttCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGt                                                                     <  2:146686/150‑1 (MQ=255)
                                                                          tGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAg                                                                  <  2:522785/150‑1 (MQ=255)
                                                                          tGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAg                                                                  <  2:1782954/150‑1 (MQ=255)
                                                                          agcccACGGCCACAGATGAGCGAAGAATTCCGGAATGCCTAGCAGTCGCTCAGCCGACTGGTGTCACCCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCGGAGGCGCCCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAg                                                                  <  2:604730/146‑1 (MQ=255)
                                                                            cTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGgt                                                                >  1:932013/1‑150 (MQ=255)
                                                                                 ggCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACa                                                           <  2:932013/150‑1 (MQ=255)
                                                                                   ccGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGc                                                         <  1:1469718/150‑1 (MQ=255)
                                                                                    cGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGcc                                                        >  1:1046478/1‑150 (MQ=255)
                                                                                    cGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTATACAGcc                                                        <  2:866018/150‑1 (MQ=255)
                                                                                        gATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCctgct                                                    <  1:154312/150‑1 (MQ=255)
                                                                                            aGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTgcgc                                                <  1:1208155/150‑1 (MQ=255)
                                                                                               gaagaaTTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAAc                                             >  2:1167903/1‑150 (MQ=255)
                                                                                                aagaaTTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAAcc                                            >  2:1786637/1‑150 (MQ=255)
                                                                                                    aTTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACcgcgc                                        <  1:1687277/150‑1 (MQ=255)
                                                                                                     ttCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCACCcgtgcg                                       <  2:1044444/150‑1 (MQ=255)
                                                                                                     ttCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACcgcgcg                                       >  2:1343876/1‑150 (MQ=255)
                                                                                                     ttCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACcgcgcg                                       >  2:864041/1‑150 (MQ=255)
                                                                                                     ttCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACcgcgcg                                       >  1:868965/1‑150 (MQ=255)
                                                                                                      tCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCgcgcgc                                      >  1:1654559/1‑150 (MQ=255)
                                                                                                       ccGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGca                                     <  2:1833994/150‑1 (MQ=255)
                                                                                                        cGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGcac                                    >  2:568347/1‑150 (MQ=255)
                                                                                                         gCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAg                                                                  >  2:1034590/1‑119 (MQ=255)
                                                                                                         gCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCGGCGTATCGAAGGGCAGGACGGCCAGCGTCAg                                                                  <  1:1034590/119‑1 (MQ=255)
                                                                                                           aaTGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGc                                 >  1:1604575/1‑150 (MQ=255)
                                                                                                             tGCCCAGCAGGCGCTCAGCCGACTTGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCaa                               <  1:837141/150‑1 (MQ=255)
                                                                                                                       tCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGcaca                                   >  1:803867/1‑136 (MQ=255)
                                                                                                                       tCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGcaca                                   <  2:803867/136‑1 (MQ=255)
                                                                                                                         gCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCct                   >  1:1339696/1‑150 (MQ=255)
                                                                                                                          cTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGCCCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCctg                  <  1:721000/150‑1 (MQ=255)
                                                                                                                             aGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCctgctg               <  1:175074/150‑1 (MQ=255)
                                                                                                                              gCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTgg              <  1:1167903/150‑1 (MQ=255)
                                                                                                                              gCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTgg              <  2:1046478/150‑1 (MQ=255)
                                                                                                                                   cTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAg         <  1:864041/150‑1 (MQ=255)
                                                                                                                                      gtgtCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAAGACGGCCAACGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCg      >  2:1611278/1‑150 (MQ=255)
                                                                                                                                         tCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATTGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTAc   >  2:893611/1‑150 (MQ=255)
                                                                                                                                         tCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTAc   <  1:1796312/150‑1 (MQ=255)
                                                                                                                                         tCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTAc   <  2:1623630/150‑1 (MQ=255)
                                                                                                                                          cACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTAcc  <  1:1163962/150‑1 (MQ=255)
                                                                                                                                           aCGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTgg                          >  2:581706/1‑125 (MQ=255)
                                                                                                                                           aCGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTgg                          <  1:581706/125‑1 (MQ=255)
                                                                                                                                             gCCCGGCCGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACa                                                           >  1:1357665/1‑90 (MQ=255)
                                                                                                                                             gCCCGGCCGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACa                                                           <  2:1357665/90‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GACTGAAGCTGCTGCAGCGCACCGGCAGCCTGAATGCCGCACCGGACCCGATCAACCGCTGCCAGGGCGCGTTCTGCTCACGGCCGCAGATGAGCGAAGAATTCCGCAATGCCGAGCAGTCGCTCAGCCGACTGGTGTCACGCCCGGCGGCCGGGCTGAAGGTGATCAACCAGTTGCCCGAGGCGACCATGCTGCGTATCGAAGGGCAGGACGGCCAGCGTCAGGTGTACAGCCTGCTGCGCAACCGCGCGCACAGCAACGTGGCGTTCCTGCTGGGTGAGGCGTACC  >  NC_002947/2717058‑2717345

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: