Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,481 C→G 6.0% intergenic (+36/+51) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,4810CG 6.0% 951.6 / 34.3 276intergenic (+36/+51)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base C (121/136);  new base G (7/10);  total (129/147)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 283 total aligned reads displayed.

TGATACCGCTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCA  >  NC_000913/3959211‑3959758
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                     
tGAGACCGATTAGTACGGTAACGATCTGTGCTCTTACGCTGGTTTGCAGTTGCGGAAACCGTCTATGGCAGAGCGGCAACCGGGCGCCCTGGGTAAAGCTATTGCGGAAGGGGCTGTAGATAACGGGCAACTGCGTGAGGTGAACGAAGCGATTCGTAGCCGTTCGATGTAGCAGGTAGGTAAGAAACTGCGCGGCTATATGGCAGGTGTGAAACGTATTGATGTGGCGGGTTAAGTGCGCACTGATGCCCTCACCCGGACGCTCTCCCACaggg                                                                                                                                                                                                                                                                                   <  4:191550/275‑1 (MQ=255)
 gATACCGCTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCCCCCACaggg                                                                                                                                                                                                                                                                                   >  2:464693/1‑274 (MQ=255)
        cTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaccaggacggggg                                                                                                                                                                                                                                                                            >  2:404325/1‑263 (MQ=255)
                              ctctTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                                                                                                                                                             >  2:29099/1‑251 (MQ=255)
                                           gtgCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                                                                                                                                >  2:194145/1‑251 (MQ=255)
                                                                                                                           cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTTAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                                                                                                                                                      <  8:1392624/149‑1 (MQ=255)
                                                                                                                           cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                                                                                                                                                                      <  5:1176085/149‑1 (MQ=255)
                                                                                                                            gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGAAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACag                                                                                                                                                                                                                                                                                     >  7:711211/1‑149 (MQ=255)
                                                                                                                              gcaaatggggagtgtgAACGAAGCGGTTGCCGGCGTTGCGAGTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGAAATGAAACTACTTGCGGCTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCCCCAACAGGGAGAGGTAGAAAACACTCAAGGCCTTCTCCTTGAGAAGGCGTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                                                                                                                      <  4:295726/262‑1 (MQ=255)
                                                                                                                              gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggg                                                                                                                                                                                                                                                                                   <  5:1120802/149‑1 (MQ=255)
                                                                                                                               cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                                                                                                                                                                  >  8:1115941/1‑149 (MQ=255)
                                                                                                                                  cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                                                                                                                                                                               <  6:402728/149‑1 (MQ=255)
                                                                                                                                  cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                                                                                                                                                                               <  8:711211/149‑1 (MQ=255)
                                                                                                                                   tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                                                                                                                                                              <  7:1115941/149‑1 (MQ=255)
                                                                                                                                         aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGGCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagggtcggg                                                                                                                                                                                                                                                                        <  8:1512281/149‑6 (MQ=255)
                                                                                                                                          tgtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                                                                                                                                                                                       <  7:263326/149‑1 (MQ=255)
                                                                                                                                             gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAAAAAACCACTTTAt                                                                                                                                      >  2:30396/1‑275 (MQ=255)
                                                                                                                                                 gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                                                                                                                                                             <  7:237853/136‑1 (MQ=255)
                                                                                                                                                 gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                                                                                                                                >  6:1272012/1‑149 (MQ=255)
                                                                                                                                                 gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                                                                                                                                >  6:439158/1‑149 (MQ=255)
                                                                                                                                                  aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                                                                                                                               <  6:746593/149‑1 (MQ=255)
                                                                                                                                                  aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGATTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                                                                                                                               <  8:159409/149‑1 (MQ=255)
                                                                                                                                                       aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                                                                                                                                                              <  6:121120/129‑1 (MQ=255)
                                                                                                                                                       aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCtt                                                                                                                                                                                                                                                          <  7:746995/149‑1 (MQ=255)
                                                                                                                                                        ttCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                                                                                                                                                                                         >  8:339708/1‑149 (MQ=255)
                                                                                                                                                          cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                                                                                                                                                       <  7:339708/149‑1 (MQ=255)
                                                                                                                                                            cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAACACTCAAGGCCTTCTCCTg                                                                                                                                                                                                                                                    <  7:616649/149‑1 (MQ=255)
                                                                                                                                                                  tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                                                                                                                                     <  6:675766/143‑1 (MQ=255)
                                                                                                                                                                     gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                                                                                                                            >  6:1259578/1‑149 (MQ=255)
                                                                                                                                                                      aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAGCACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                                                                                                                                                           <  5:439158/149‑1 (MQ=255)
                                                                                                                                                                       ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                                                                                                                          <  5:986391/149‑1 (MQ=255)
                                                                                                                                                                        tgtgCAGGTAGGTAAGAACCGTCGCGGCTATACGACAGATTTCAACCGTTTGGCTGTTGCGGGTTAAGTGGGCGCTGATGCCCTCACCCCGACCCTCTCCCCACGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg                                                                                                            <  4:404321/271‑1 (MQ=255)
                                                                                                                                                                            cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                                                                                                                                     <  5:1638683/149‑1 (MQ=255)
                                                                                                                                                                              ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgt                                                                                                                                                                                                                                   <  5:418668/149‑6 (MQ=255)
                                                                                                                                                                              ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACACGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                                                                                                                                                                                                                                   <  6:510363/149‑1 (MQ=255)
                                                                                                                                                                               gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                                                                                                                                  <  8:132827/149‑1 (MQ=255)
                                                                                                                                                                               gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                                                                                                                                                              >  6:874240/1‑105 (MQ=255)
                                                                                                                                                                               gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                                                                                                                                  >  8:1029876/1‑149 (MQ=255)
                                                                                                                                                                                   gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                                                                                                                                                                                         <  2:707371/154‑1 (MQ=255)
                                                                                                                                                                                    taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                                                                                                                             >  6:1555686/1‑149 (MQ=255)
                                                                                                                                                                                    taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                                                                                                                             >  7:306707/1‑149 (MQ=255)
                                                                                                                                                                                      aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                                                                                                                           <  8:356601/149‑1 (MQ=255)
                                                                                                                                                                                              cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                                                                                                                                      <  1:35852/146‑1 (MQ=255)
                                                                                                                                                                                               gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTACGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                                                                                                                                  <  7:149530/149‑1 (MQ=255)
                                                                                                                                                                                                      tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                                                                                                                                                                                                                                    <  6:1461020/76‑1 (MQ=255)
                                                                                                                                                                                                      tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                                                                                                                                                                                                           >  7:485212/1‑149 (MQ=255)
                                                                                                                                                                                                        tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                                                                                                                                          <  8:554636/116‑1 (MQ=255)
                                                                                                                                                                                                        tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                                                                                                                         <  6:1704383/149‑1 (MQ=255)
                                                                                                                                                                                                         gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                        <  5:1315802/149‑1 (MQ=255)
                                                                                                                                                                                                           cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                                                                                                                                      <  5:1259578/149‑1 (MQ=255)
                                                                                                                                                                                                             gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                                                                                                                                                                                       <  3:64979/82‑1 (MQ=255)
                                                                                                                                                                                                              atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATggg                                                                     >  2:572887/1‑275 (MQ=255)
                                                                                                                                                                                                               tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                                                                                                                               <  5:683386/88‑1 (MQ=255)
                                                                                                                                                                                                                atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                                                                                                                                 <  8:942002/149‑1 (MQ=255)
                                                                                                                                                                                                                atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                                                                                                                                 <  5:1198149/149‑1 (MQ=255)
                                                                                                                                                                                                                  gAAACGTATTGCTGTTGGGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                                                                                                                                   <  4:734332/177‑1 (MQ=255)
                                                                                                                                                                                                                   aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                                                                                                                                                           <  7:1410605/88‑1 (MQ=255)
                                                                                                                                                                                                                   aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                                                                                                                              >  5:878847/1‑149 (MQ=255)
                                                                                                                                                                                                                    aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                                                                                                                                                              <  8:647092/68‑1 (MQ=255)
                                                                                                                                                                                                                    aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                                                                                                                             <  7:434759/149‑1 (MQ=255)
                                                                                                                                                                                                                     aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                                                                                                                                                              >  7:1287643/1‑67 (MQ=255)
                                                                                                                                                                                                                      cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                                                                                                                                     >  6:1527733/1‑139 (MQ=255)
                                                                                                                                                                                                                      cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                                                                                                                                     <  5:1527733/139‑1 (MQ=255)
                                                                                                                                                                                                                       gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                          <  8:485212/149‑1 (MQ=255)
                                                                                                                                                                                                                                  gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                                                                                                                                               <  7:1198380/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                                                                                                                                                     <  7:555715/43‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                               >  6:1415691/1‑81 (MQ=255)
                                                                                                                                                                                                                                      gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                               <  5:1415691/81‑1 (MQ=255)
                                                                                                                                                                                                                                      gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGata                                                                                                                                                                                                    <  2:737723/124‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                                                                                                                                                      >  6:114250/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                                                                                                                                                                  >  5:614456/1‑149 (MQ=255)
                                                                                                                                                                                                                                                cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCAc                                                                                         >  2:248484/1‑221 (MQ=255)
                                                                                                                                                                                                                                                 gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAg                                                                                                                                                                >  5:344101/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGccc                             <  1:165426/279‑1 (MQ=255)
                                                                                                                                                                                                                                                    gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCAccgcc                                                                                     <  4:248481/221‑1 (MQ=255)
                                                                                                                                                                                                                                                      tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAGTGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                                                                                                                                                           >  7:40192/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aaccccgtccccGGCCCTCAGCCACGGGGAAATGGGGAAAAAGATAAAGGCGTTTGCCTGGGAAGGGCCTTGCGATGAGTTGCGGTACAAAAACGTAATGGTGGGATATCCGAACTGGGAATGCAGCGGGACGGTCGGCACCGGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCGGCTGAAACTAGCGGAATTCACCTAACTCACGGCCGCGTTTGCCTGATAGGGAAAAGGGGTATATCTACGCCAGCTTGCCGAAATCGGccc                             <  4:560122/267‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                                                                                                                        >  5:1474485/1‑132 (MQ=255)
                                                                                                                                                                                                                                                           tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                                                                                                                          >  8:1243598/1‑97 (MQ=255)
                                                                                                                                                                                                                                                           tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                                                                                                                      >  6:338536/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                                                                                                                  >  7:1051/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                                                                                                                             <  6:202076/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                 >  7:1690310/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGCTGCGGTGCAGCACCTTAATGATGTGATATCCGAACTGAGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgagac                                                                                                                                                 >  5:1230400/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   aCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATg                                                                       >  3:87386/1‑220 (MQ=255)
                                                                                                                                                                                                                                                                    cccttccccACAGGTAGAGGAGAAAATACTCAATGCCTTCTCCTGTGGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGGGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTGATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTGGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCt              <  2:49386/269‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                                                                                                                            >  8:809063/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCt              >  2:629334/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                       tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCg                                                                                                                                                                       <  3:287284/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                                                                                                                              >  5:948994/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                       tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                                                                                                                          >  6:972543/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                                                                                                                          >  8:1491245/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg                                                                                                                            <  2:118687/163‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                                                                                                                                         <  8:1461665/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                                                                                                                             <  1:83163/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt                                                                                                                                        >  5:195447/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCga           >  4:67810/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa                                                                                                                                       <  8:224083/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                                                                                                                      >  5:67761/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTTCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCTAGAAGATccc  >  1:268697/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                     
TGATACCGCTGAGTACGGTAACTATCTGTTCTCTTACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCA  >  NC_000913/3959211‑3959758

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