Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,064 C→G 8.8% L305L (CTG→CTC 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,0640CG 8.8% 586.6 / 36.1 182L305L (CTG→CTCetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (79/87);  new base G (7/9);  total (86/96)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 187 total aligned reads displayed.

GTGGTGACACCAACAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGC  >  NC_002947/4747916‑4748213
                                                                                                                                                    |                                                                                                                                                     
gtggtgACACCAACAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCa                                                                                                                                                      <  2:2400837/150‑1 (MQ=255)
          cAACAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGacc                                                                                                                                            >  2:4383730/1‑148 (MQ=255)
          cAACAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAAc                                                                                                                                            >  1:800680/1‑150 (MQ=255)
            aCAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCg                                                                                                                                          >  1:496254/1‑150 (MQ=255)
                  ccACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTcc                                                                                                                                                       <  1:2626494/131‑1 (MQ=255)
                  ccACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTcc                                                                                                                                                       <  1:5352171/131‑1 (MQ=255)
                  ccACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTcc                                                                                                                                                       >  2:2626494/1‑131 (MQ=255)
                  ccACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGa                                                                                                                                    >  1:1435009/1‑150 (MQ=255)
                     cacaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCACAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACa                                                                                                                                 <  1:3790924/150‑1 (MQ=255)
                      acaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAg                                                                                                                                <  2:4442761/150‑1 (MQ=255)
                        aTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGc                                                                                                                                         <  2:5377370/139‑1 (MQ=255)
                          gCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                            >  2:1958153/1‑150 (MQ=255)
                           ccAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                            <  2:3837533/149‑1 (MQ=255)
                           ccAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                            >  1:3837533/1‑149 (MQ=255)
                               cACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGAc                                                                                                                       <  1:4789998/150‑1 (MQ=255)
                                acaatGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGAcc                                                                                                                      >  2:3352511/5‑150 (MQ=255)
                                     gCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGcc                                                                                                                 <  1:2729701/150‑1 (MQ=255)
                                        cGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTa                                                                                                              <  2:1987464/150‑1 (MQ=255)
                                          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)
                                                   aCGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAACCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCt                                                                                                   <  1:976199/150‑1 (MQ=255)
                                                      tCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                <  1:4711274/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)
                                                                 cccTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGct                                                                                     <  1:3485079/150‑1 (MQ=255)
                                                                 cccTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGct                                                                                     >  1:3074149/1‑150 (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)
                                                                             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)
                                                                               aaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAAc                                                                                                                                            >  2:3892374/1‑81 (MQ=255)
                                                                               aaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAAc                                                                                                                                            <  1:3892374/81‑1 (MQ=255)
                                                                               aaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgtt                                                                       >  1:3616970/1‑150 (MQ=255)
                                                                                 aaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGa                                                                     >  1:2018949/1‑150 (MQ=255)
                                                                                   aCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                   <  2:3874176/150‑1 (MQ=255)
                                                                                     ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                 >  2:2139218/1‑150 (MQ=255)
                                                                                     ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                 >  2:5214400/1‑150 (MQ=255)
                                                                                     ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                 >  1:4204478/1‑150 (MQ=255)
                                                                                     ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGAAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGgc                                                                                        >  2:3691247/1‑127 (MQ=255)
                                                                                      cGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGcc                                                                <  2:4988922/150‑1 (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)
                                                                                                 aCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTGTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGTGGCGATCaacaa                                                     <  2:3769342/150‑5 (MQ=255)
                                                                                                  cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                    <  2:2499312/150‑1 (MQ=255)
                                                                                                    aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                  <  1:244121/150‑1 (MQ=255)
                                                                                                    aGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                  <  1:3009118/150‑1 (MQ=255)
                                                                                                         cGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCtt                                             >  1:2210442/1‑150 (MQ=255)
                                                                                                              tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                        <  2:358682/150‑1 (MQ=255)
                                                                                                              tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                        <  1:2654415/150‑1 (MQ=255)
                                                                                                                 ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                  <  1:1215883/137‑1 (MQ=255)
                                                                                                                 ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                     <  2:5267311/150‑1 (MQ=255)
                                                                                                                     aaGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCa                                 >  2:3643660/1‑150 (MQ=255)
                                                                                                                      aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                >  2:4747843/1‑150 (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:1296438/150‑1 (MQ=255)
                                                                                                                            gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                          <  1:3358951/150‑1 (MQ=255)
                                                                                                                             cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                  <  2:5054580/125‑1 (MQ=255)
                                                                                                                             cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                  >  1:5054586/1‑125 (MQ=255)
                                                                                                                               ttataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGAt                       >  1:537043/6‑150 (MQ=255)
                                                                                                                                tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                      >  2:117786/5‑150 (MQ=255)
                                                                                                                                tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                      <  2:2722572/146‑1 (MQ=255)
                                                                                                                                tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                      >  1:5213295/5‑150 (MQ=255)
                                                                                                                                 ataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcg                     >  2:3528389/4‑150 (MQ=255)
                                                                                                                                   aaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcgca                   >  1:3320286/2‑149 (MQ=255)
                                                                                                                                        gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCATGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga              <  2:5324436/150‑1 (MQ=255)
                                                                                                                                        gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACt                                                       >  1:5399184/1‑109 (MQ=255)
                                                                                                                                        gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga              <  2:1034958/150‑1 (MQ=255)
                                                                                                                                         aCCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGgag             <  2:3113241/150‑1 (MQ=255)
                                                                                                                                         aCCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACGTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGgag             >  2:5362959/1‑150 (MQ=255)
                                                                                                                                          gcactttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCgaga            <  1:3398296/145‑1 (MQ=255)
                                                                                                                                          gcactttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCgaga            <  1:3528389/145‑1 (MQ=255)
                                                                                                                                           cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt           >  2:2757719/1‑150 (MQ=255)
                                                                                                                                           cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt           >  1:2943323/1‑150 (MQ=255)
                                                                                                                                           cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt           >  1:4526428/1‑150 (MQ=255)
                                                                                                                                            actttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtg          <  1:3972162/147‑2 (MQ=255)
                                                                                                                                            aGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATa          <  1:284793/150‑1 (MQ=255)
                                                                                                                                             gCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc         <  2:5261914/150‑1 (MQ=255)
                                                                                                                                             gCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc         <  2:2885788/150‑1 (MQ=255)
                                                                                                                                               ttttCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCg       >  2:70870/1‑150 (MQ=255)
                                                                                                                                                  tCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACg    <  2:3288491/150‑1 (MQ=255)
                                                                                                                                                  tCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACg    <  2:562284/150‑1 (MQ=255)
                                                                                                                                                    cAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGc  >  1:1567511/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTGACACCAACAGCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGC  >  NC_002947/4747916‑4748213

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: