Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,079 T→G 9.5% E300D (GAA→GAC) ‡ 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,0790TG 9.5% 535.4 / 36.4 177E300D (GAA→GAC) ‡etfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (76/84);  new base G (8/9);  total (84/93)
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.94e-01

Only 100 of 181 total aligned reads displayed.

GCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGG  >  NC_002947/4747931‑4748228
                                                                                                                                                    |                                                                                                                                                     
gcgcCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAACCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCtt                                                                                                                                                      <  1:3533759/150‑1 (MQ=255)
    cacacaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGaa                                                                                                                                                  <  1:3178289/150‑1 (MQ=255)
      cacaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCACAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACa                                                                                                                                                <  1:3790924/150‑1 (MQ=255)
       acaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAg                                                                                                                                               <  2:4442761/150‑1 (MQ=255)
       acaTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAg                                                                                                                                               <  1:1403718/150‑1 (MQ=255)
         aTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGt                                                                                                                                             >  2:2667839/1‑150 (MQ=255)
         aTGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGt                                                                                                                                             >  1:2923265/1‑150 (MQ=255)
           gCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                                           >  2:1958153/1‑150 (MQ=255)
            ccAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                                           >  1:3837533/1‑149 (MQ=255)
            ccAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCa                                                                                                                                           <  2:3837732/149‑1 (MQ=255)
              aGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCg                                                                                                                                        <  2:5248661/150‑1 (MQ=255)
                cACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGAc                                                                                                                                      <  1:4789998/150‑1 (MQ=255)
                    cTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAgg                                                                                                                                  <  1:895666/150‑1 (MQ=255)
                      gCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGcc                                                                                                                                <  1:2729701/150‑1 (MQ=255)
                         cGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTa                                                                                                                             <  2:1987464/150‑1 (MQ=255)
                           ttCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGt                                                                                                                           <  2:4923800/150‑1 (MQ=255)
                                 gcAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCa                                                                                                                     <  1:456853/150‑1 (MQ=255)
                                   aaCGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCAcc                                                                                                                   >  2:2861684/1‑150 (MQ=255)
                                      gTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTgg                                                                                                                >  1:4466778/1‑150 (MQ=255)
                                       tCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                               <  1:4711274/150‑1 (MQ=255)
                                       tCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGa                                                                                                               <  2:1621977/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                                                                                                        >  1:2543295/1‑150 (MQ=255)
                                                  cccTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCGTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGct                                                                                                    >  1:2436521/1‑150 (MQ=255)
                                                   ccTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGctt                                                                                                   >  2:2158416/1‑150 (MQ=255)
                                                    cTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGcttc                                                                                                  <  2:3331918/150‑1 (MQ=255)
                                                          ccATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTc                                                                                            <  1:3635784/150‑1 (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:1640574/1‑150 (MQ=255)
                                                               aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                       >  2:1296438/1‑150 (MQ=255)
                                                               aaaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgt                                                                                       <  1:1640574/150‑1 (MQ=255)
                                                                aaaaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCttgtt                                                                                      >  1:3616970/1‑150 (MQ=255)
                                                                  aaaCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGa                                                                                    >  1:2018949/1‑150 (MQ=255)
                                                                    aCCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                  <  2:3812588/150‑1 (MQ=255)
                                                                      ccGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGc                                                                                >  2:1071215/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)
                                                                             tCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATc                                                                         >  1:300917/1‑150 (MQ=255)
                                                                               aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                       >  1:1181367/1‑150 (MQ=255)
                                                                               aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                       >  1:1352046/1‑150 (MQ=255)
                                                                               aaaaCGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCAc                                                                       >  1:1010043/1‑150 (MQ=255)
                                                                                   cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                                   <  2:2499312/150‑1 (MQ=255)
                                                                                   cGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTg                                                                   <  1:2990039/150‑1 (MQ=255)
                                                                                    gAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGTGGCGATCaacaagg                                                                  <  1:3250574/150‑7 (MQ=255)
                                                                                           gggTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCttt                                                           <  1:1495748/150‑1 (MQ=255)
                                                                                            ggTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTc                                                          <  2:158387/150‑1 (MQ=255)
                                                                                              ttCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCAt                                                        >  2:1156818/1‑150 (MQ=255)
                                                                                               tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                                       <  2:358794/150‑1 (MQ=255)
                                                                                               tCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATg                                                       <  1:2654415/150‑1 (MQ=255)
                                                                                                 tttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGcc                                                     >  1:5240751/1‑150 (MQ=255)
                                                                                                  ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                 >  2:1215883/1‑137 (MQ=255)
                                                                                                  ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                 <  1:1215883/137‑1 (MQ=255)
                                                                                                  ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                    <  2:1003954/150‑1 (MQ=255)
                                                                                                  ttATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCg                                                    <  2:5267311/150‑1 (MQ=255)
                                                                                                       aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                               >  2:3670292/1‑150 (MQ=255)
                                                                                                       aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                               >  2:3358951/1‑150 (MQ=255)
                                                                                                       aGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCACTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAg                                               <  2:918521/150‑1 (MQ=255)
                                                                                                        gTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                             >  1:629906/1‑119 (MQ=255)
                                                                                                        gTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                             <  2:629906/119‑1 (MQ=255)
                                                                                                         tGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                             <  1:5023470/118‑1 (MQ=255)
                                                                                                           aaGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGc                                           <  1:2158416/150‑1 (MQ=255)
                                                                                                            aGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCt                                          <  2:2062025/150‑1 (MQ=255)
                                                                                                             gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                                         <  1:1296438/150‑1 (MQ=255)
                                                                                                             gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                                         <  1:3358951/150‑1 (MQ=255)
                                                                                                             gCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTg                                         >  1:3502394/1‑150 (MQ=255)
                                                                                                              cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                 >  1:5054586/1‑125 (MQ=255)
                                                                                                              cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                 <  2:5054586/125‑1 (MQ=255)
                                                                                                              cAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGa                                                                 <  2:5054580/125‑1 (MQ=255)
                                                                                                                ttataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGAt                                      >  1:537043/6‑150 (MQ=255)
                                                                                                                 tataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATc                                     <  2:2722572/146‑1 (MQ=255)
                                                                                                                  ataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcg                                    >  2:3528389/4‑150 (MQ=255)
                                                                                                                  ataaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcg                                    >  2:1464871/4‑150 (MQ=255)
                                                                                                                    aaTCAGAGCACTTTTTCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATcgca                                  >  1:3320286/2‑149 (MQ=255)
                                                                                                                     aTCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCgcgc                                 <  2:4867954/150‑1 (MQ=255)
                                                                                                                         gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCATGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga                             <  2:159147/150‑1 (MQ=255)
                                                                                                                         gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCATGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga                             <  2:5324436/150‑1 (MQ=255)
                                                                                                                         gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACt                                                                      <  2:5399184/109‑1 (MQ=255)
                                                                                                                         gACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGga                             <  2:885213/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)
                                                                                                                            cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                             <  1:5162484/99‑1 (MQ=255)
                                                                                                                            cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCAc                                                                             >  2:5162484/1‑99 (MQ=255)
                                                                                                                            cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt                          >  2:2757719/1‑150 (MQ=255)
                                                                                                                            cAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCATAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGAt                          >  1:4103063/1‑150 (MQ=255)
                                                                                                                             actttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtg                         <  1:3972193/147‑2 (MQ=255)
                                                                                                                             actttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtg                         <  1:3972162/147‑2 (MQ=255)
                                                                                                                              gCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATAc                        <  2:2885788/150‑1 (MQ=255)
                                                                                                                              ctttttCGAGCTCCGGCAACACGGTGAAGAGATCGCCGACGAGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGAtgc                        <  1:4231028/148‑3 (MQ=255)
                                                                                                                                   tCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACg                   <  2:3288491/150‑1 (MQ=255)
                                                                                                                                              aaCCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATc                                                                                  >  1:929120/1‑76 (MQ=255)
                                                                                                                                                   cTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGc                                                  >  1:2671992/1‑103 (MQ=255)
                                                                                                                                                    ttCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTgg  >  2:3710698/1‑150 (MQ=255)
                                                                                                                                                    ttCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTgg  >  1:4947202/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCCACACATGCCAGCACCCTGCCCGTTCCGCGCAACGTCATCCCTCACCCCTCAAGCCATAAAAAAACCCGACCCTCAAAACGAGGGCCGGGTTCTTTATAAGTGAAGCAGGCGGATCAGACCAGCTTTTCCAGCTCCGGAACCGCTTCGAACAGGTCAGCGACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGG  >  NC_002947/4747931‑4748228

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: