Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,106 G→A 15.1% D291D (GAC→GAT etfA ← electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,1060GA15.1% 460.4 / 62.3 188D291D (GAC→GATetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (81/78);  new base A (13/15);  total (94/94)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAC  >  NC_002947/4747958‑4748255
                                                                                                                                                    |                                                                                                                                                     
ttCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGt                                                                                                                                                      <  2:4923800/150‑1 (MQ=255)
gtCCGCGCAACGTCATCCCTCACCCCTCACGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGt                                                                                                                                                      <  2:1604950/149‑1 (MQ=255)
  ccGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCg                                                                                                                                                    <  1:3185003/150‑1 (MQ=255)
    gcgcAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGc                                                                                                                                                  <  2:1463359/150‑1 (MQ=255)
        aaCGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCAcc                                                                                                                                              >  2:2861684/1‑150 (MQ=255)
            tCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                                                          <  1:4711339/150‑1 (MQ=255)
                cccTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGAt                                                                                                                                      >  1:1032413/1‑150 (MQ=255)
                cccTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGAt                                                                                                                                      >  1:4601495/1‑150 (MQ=255)
                  cTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCg                                                                                                                                    >  2:2208929/1‑150 (MQ=255)
                   tCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCgg                                                                                                                                   <  2:1110791/150‑1 (MQ=255)
                   tCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCgg                                                                                                                                   >  1:2543295/1‑150 (MQ=255)
                   tCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCGTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCgg                                                                                                                                   >  1:4995495/1‑150 (MQ=255)
                       cccTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGct                                                                                                                               <  1:3485079/150‑1 (MQ=255)
                         cTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttc                                                                                                                             >  1:3974310/1‑150 (MQ=255)
                          tCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttct                                                                                                                            >  2:202351/1‑150 (MQ=255)
                          tCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttct                                                                                                                            <  1:4795369/150‑1 (MQ=255)
                               ccATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTc                                                                                                                       <  1:3635784/150‑1 (MQ=255)
                                 aTAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCt                                                                                                                     >  2:5306879/1‑150 (MQ=255)
                                 aTAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCGGGAAGATCGGCGCTTCTGCGTCCt                                                                                                                     >  2:4594213/1‑150 (MQ=255)
                                   aaaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttg                                                                                                                   >  2:3883854/1‑150 (MQ=255)
                                   aaaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttg                                                                                                                   >  1:2155149/1‑150 (MQ=255)
                                    aaaaaaCCCGACCCTCAAACCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                                                  <  2:1714830/150‑1 (MQ=255)
                                    aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                                                  >  2:1296438/1‑150 (MQ=255)
                                    aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                                                  <  1:1640574/150‑1 (MQ=255)
                                    aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                                                  >  2:1640574/1‑150 (MQ=255)
                                     aaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgtt                                                                                                                 >  1:3616970/1‑150 (MQ=255)
                                       aaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGa                                                                                                               >  1:2018949/1‑150 (MQ=255)
                                         aCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                                             <  2:3812588/150‑1 (MQ=255)
                                           ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGgc                                                                                                                                  <  1:3691247/127‑1 (MQ=255)
                                           ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                                                           >  2:2139218/1‑150 (MQ=255)
                                           ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                                                           >  2:5214400/1‑150 (MQ=255)
                                           ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                                                           >  2:1071215/1‑150 (MQ=255)
                                            cGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGcc                                                                                                          <  2:4988922/150‑1 (MQ=255)
                                                  tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                    >  1:300917/1‑150 (MQ=255)
                                                  tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                                                    >  2:2287221/1‑150 (MQ=255)
                                                    aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                                                  >  1:4361257/1‑150 (MQ=255)
                                                    aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                                                  >  1:1352046/1‑150 (MQ=255)
                                                    aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                                                  >  1:1010043/1‑150 (MQ=255)
                                                    aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                                                  >  1:1181367/1‑150 (MQ=255)
                                                     aaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACt                                                                                                 >  2:2556498/1‑150 (MQ=255)
                                                       aCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTGTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGTGGCGATCaacaa                                                                                               <  2:3769342/150‑5 (MQ=255)
                                                        cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                                                              <  1:2990039/150‑1 (MQ=255)
                                                        cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                                                              <  2:2499312/150‑1 (MQ=255)
                                                        cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                                                              <  2:280704/150‑1 (MQ=255)
                                                         gAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGTGGCGATCaacaagg                                                                                             <  1:3250574/150‑7 (MQ=255)
                                                          aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            <  1:3352511/150‑1 (MQ=255)
                                                          aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            <  1:244121/150‑1 (MQ=255)
                                                          aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTGTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            <  2:2896503/150‑1 (MQ=255)
                                                               cGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCtt                                                                                       >  1:2210442/1‑150 (MQ=255)
                                                                gggTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCttt                                                                                      <  1:1495748/150‑1 (MQ=255)
                                                                    tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                                                                  <  2:358794/150‑1 (MQ=255)
                                                                    tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                                                                  <  1:2654415/150‑1 (MQ=255)
                                                                       ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            <  1:1215883/137‑1 (MQ=255)
                                                                       ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                                               <  2:1003954/150‑1 (MQ=255)
                                                                        tataAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCgg                                                                              >  1:3373091/1‑150 (MQ=255)
                                                                            aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                                                          >  2:3358951/1‑150 (MQ=255)
                                                                            aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                                                          >  2:3670292/1‑150 (MQ=255)
                                                                            aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                                                          <  2:3695092/150‑1 (MQ=255)
                                                                            aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCACTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                                                          <  2:918401/150‑1 (MQ=255)
                                                                             gTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                                                        >  1:629906/1‑119 (MQ=255)
                                                                             gTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTGTGGAGTCTTTCATGCCGGCCAgg                                                                         >  2:1508541/1‑150 (MQ=255)
                                                                              tGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                                                        <  1:5023470/118‑1 (MQ=255)
                                                                              tGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGt                                                                        >  2:1589933/1‑150 (MQ=255)
                                                                              tGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGt                                                                        >  2:5126028/1‑150 (MQ=255)
                                                                                  gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                                                                    >  1:3502394/1‑150 (MQ=255)
                                                                                   cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            >  1:5054586/1‑125 (MQ=255)
                                                                                   cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                                            <  2:5054580/125‑1 (MQ=255)
                                                                                     ttataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGAt                                                                 >  1:537043/6‑150 (MQ=255)
                                                                                      tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                                                                <  2:2722572/146‑1 (MQ=255)
                                                                                      tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                                                                >  1:5213295/5‑150 (MQ=255)
                                                                                      tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                                                                >  2:117786/5‑150 (MQ=255)
                                                                                       ataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcg                                                               >  2:3528389/4‑150 (MQ=255)
                                                                                          aTCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCgcgc                                                            <  2:4867954/150‑1 (MQ=255)
                                                                                              gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACt                                                                                                 <  2:5399184/109‑1 (MQ=255)
                                                                                               aCCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGgag                                                       <  2:3113241/150‑1 (MQ=255)
                                                                                               aCCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACGTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGgag                                                       >  2:5362959/1‑150 (MQ=255)
                                                                                                 cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt                                                     >  1:4526428/1‑150 (MQ=255)
                                                                                                  actttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtg                                                    <  1:3972193/147‑2 (MQ=255)
                                                                                                   gCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                                                   <  2:5261914/150‑1 (MQ=255)
                                                                                                   ctttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtgc                                                   <  1:4231028/148‑3 (MQ=255)
                                                                                                    cTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAcc                                                  >  2:1756416/1‑150 (MQ=255)
                                                                                                        tCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACg                                              <  2:3288491/150‑1 (MQ=255)
                                                                                                                   aaCCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                                             >  1:929120/1‑76 (MQ=255)
                                                                                                                   aaCCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                                             <  2:929120/76‑1 (MQ=255)
                                                                                                                        cTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGc                                                                             >  1:2671992/1‑103 (MQ=255)
                                                                                                                              aCAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCaa                        >  2:1535569/1‑150 (MQ=255)
                                                                                                                                agatCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGAt                                         <  1:2605059/128‑1 (MQ=255)
                                                                                                                                  atCACCGACGAGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAAccac                    <  2:578302/149‑4 (MQ=255)
                                                                                                                                   tCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATc                   >  1:768154/1‑150 (MQ=255)
                                                                                                                                   tCAGCAACCAACCCATAATCGGCAACCTGGAAGATCGGCGCTTCCTCATCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcgca                                                             >  1:3936815/1‑107 (MQ=255)
                                                                                                                                       cgacgaGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGc               <  1:4189723/150‑1 (MQ=255)
                                                                                                                                       cGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGc               <  1:4897292/150‑1 (MQ=255)
                                                                                                                                       cGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGc               <  1:2326034/150‑1 (MQ=255)
                                                                                                                                         acgaGGCCATACTCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCg             <  1:1265141/147‑1 (MQ=255)
                                                                                                                                         acgaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGccc             >  2:2617620/4‑149 (MQ=255)
                                                                                                                                           gaGGCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGt           <  1:528676/149‑1 (MQ=255)
                                                                                                                                           cAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGt           >  2:3995212/1‑150 (MQ=255)
                                                                                                                                           cAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGt           >  2:67348/1‑150 (MQ=255)
                                                                                                                                                cGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGc      <  2:800680/150‑1 (MQ=255)
                                                                                                                                                    gTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAc  <  1:2861684/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAC  >  NC_002947/4747958‑4748255

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