Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 20.0% A337T (GCC→ACC)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,1310GA20.0% 143.4 / 34.3 79A337T (GCC→ACC) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base A (9/7);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAA  >  NC_000913/1195895‑1196364
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
tCGCCGGTGGCCCGTTGCTGAAAGTTAAAACCCCGAACGCGTGCAAAGAGATCGTCATTAAAGACGGGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGACCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCGGTTTGAAGcc                                                                                                                                                                                                                                          <  1:733414/239‑1 (MQ=255)
 cggcgggggccctgtgCTGAAAGTTAAAAACCCCACCACTGGCAAGGAGATCGGCAATAAAGACGGTATTGCTGAGGCAGTCCTGGAACAGATCCTGCTGCGTCCGGCTGAATATGATGGTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcc                                                                                                                                                                                                                                          <  1:1113124/223‑1 (MQ=255)
    cGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                                                                                                                                                            >  2:1091021/1‑233 (MQ=255)
     ggTGGCCCTTGGTTGAAAGTTAAAAACCCGAAACCTGGAAAAGAGTACGTCATTAAAGGCGTGAGTGCGGATGCATGCCTGCAACAGTTCCTGCTGCGTCCGGCTGAATATGATGTTATTGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGGTTGAAAccaccca                                                                                                                                                                                                                                     <  1:991934/239‑1 (MQ=255)
      gttgctcgtgtggtaaaGTTAAAGACCCGAACACTGGCAAAGAGTGCGTCATTAAAGCCGTGATTGCGGATGCATTCCGTCAACAGATCCGGCTGCGACCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCGGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGccaccca                                                                                                                                                                                                                                     <  1:418618/224‑1 (MQ=255)
      gTGGCCCGGGTTTTAAAATTAAGAACCCGTGCACTGGCACAGAGGTCGGTCTTAAAGACGTTATTGCTGAGTCCTTCCCGCAACAAATCCCTCTGCGTCCGGCTGAATGTTATGGTACCGCATGTAGGAACCTGAACCGTGACTACATTTCTGGCGCGCTGGCAGCCCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcctccca                                                                                                                                                                                                                                     <  1:7337/238‑1 (MQ=255)
         gCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCCGGTTGGCGGTATCGGGATCGCCCCTGGTGCAACCATCGGTGACGAATGCGCCCTGTTTGAGGCCACCCACgt                                                                                                                                                                                                                                  >  2:793090/1‑237 (MQ=255)
         gCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGGATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGt                                                                                                                                                                                                                                 >  2:374450/1‑239 (MQ=255)
          cccggtgatGAAAGTTAAAACCCCGAACCCTGGCAAAGAGATCGTCATGTAGGAGGTTATGTCGGATGCATGCCTGGATACGGTCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGGACCTGAACGGTGACTACATTTCTTACGGCCTGGCAGCGCAGGTTGGCGGTAGCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                                                                                                                                                                <  1:389229/231‑1 (MQ=255)
           ccGTGGCTGAAAGTTAAAAACCCGAACCCTGGCAAAGGGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                                                                                                                                                                <  2:1298479/238‑1 (MQ=255)
              tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGAAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGGACTg                                                                                                                                                                                                                             >  2:87230/1‑238 (MQ=255)
              tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGGATCGCCCCTGTGGCAACCATCGGGGACGAATGCGCCCTGTTTGAAGcgtgtctcttataaac                                                                                                                                                                                                                            >  2:169601/1‑224 (MQ=255)
              tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGCGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTg                                                                                                                                                                                                                             >  2:347148/1‑238 (MQ=255)
                gCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTgcg                                                                                                                                                                                                                           <  2:1281387/238‑1 (MQ=255)
                     aatttAAAAACCCGAACACTGGCAAAGGACTCGTCATTAAATGCGTTATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTTGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGaa                                                                                                                                                                                                                      <  2:646220/235‑1 (MQ=255)
                         ttaaacCCCGAAACCTGGGAAAGGGATCGTTCTTAAAAGCGTTATTGCGGATGCGTTCCTGCAACAGATCCTGGTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATg                                                                                                                                                                                                                 <  1:385789/237‑1 (MQ=255)
                              aCCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGt                                                                                                                                                                                                            <  2:263141/239‑1 (MQ=255)
                                  gAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAgg                                                                                                                                                                                                        >  2:704470/1‑239 (MQ=255)
                                    acccTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGAc                                                                                                                                                                                                      <  2:42345/236‑1 (MQ=255)
                                              agagATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGGGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGCCGGTATCGGGATCGCCCCTCGTGCAAACATCTGTGACTAGAGCGCCCGGTTTGTAACCCCCCCCGGGACCGAGCCCAAAAAAGCCAGACAGGCAAAAGaaaaac                                                                                                                                                                                            >  1:111905/1‑237 (MQ=255)
                                                agaTCGGCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCt                                                                                                                                                                                          <  2:939565/239‑1 (MQ=255)
                                                 gtttcttaTTAAAGACGTGATTTGTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCt                                                                                                                                                                                          <  1:1077554/231‑1 (MQ=255)
                                                 gaTCGTCATTAAAGTCTTGATTGCTGGTGCATTCCGGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACTGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATcc                                                                                                                                                                                           <  1:34260/237‑1 (MQ=255)
                                                   tgtgcATTAAAGGTGTGGTTTCGGATGCTTTCCCGAAAACAATCCGGCTGCGCTCGGCTGAATATGATGTTATCCCCTGTATGAAACTGAGCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                        <  1:659780/234‑1 (MQ=255)
                                                   tCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                        <  2:1272238/238‑1 (MQ=255)
                                                    cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGGCAAAGTAAATCCTGGc                                                                                                                                                                                       >  2:34260/1‑238 (MQ=255)
                                                    cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATAGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTCTGAAGCACCCCCCGGTACTTGCCCGAAATAAGCCGGTCCGGCAAAAGaaaaaacgggag                                                                                                                                                                                      >  1:841921/1‑231 (MQ=255)
                                                                 gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGCCGAATGCGCCCTGTTTGCAGCCACCCACGGTATTGCGCCGAAATATGCCGGTCAGGCCAAAGTAAATCCCGGCTCTATTATtctcc                                                                                                                                                                          >  1:771150/1‑237 (MQ=255)
                                                                 gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtctct                                                                                                                                                                          >  2:348225/1‑238 (MQ=255)
                                                                     ttGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                      <  2:936317/238‑1 (MQ=255)
                                                                                 tCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAtgctg                                                                                                                                                         <  2:716686/239‑1 (MQ=255)
                                                                                       aaCAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                                                                                                   <  2:887099/239‑1 (MQ=255)
                                                                                            aTCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGAGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTgcg                                                                                                                                                                                                                           >  1:905832/1‑162 (MQ=255)
                                                                                            aTCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGAGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTgcg                                                                                                                                                                                                                           <  2:905832/162‑1 (MQ=255)
                                                                                              ccTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGCAATATGCCGGTCCGGCCAAAGTAAATCCTGGCCCTTTTATTTTCCCCGCTAAGACGATTCTCGCCCACACGGGtg                                                                                                                                            >  1:477447/1‑238 (MQ=255)
                                                                                              ccTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGtt                                                                                                                                            >  2:65404/1‑239 (MQ=255)
                                                                                                       tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTTACTACATTTCTGACGACCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGCATGCGACCTGTTTTGAGCCACCCTCGGTACTGCGCCGCAATATACCGGTCAGGACAAAGAACACTCTGCCCCTACTATTTTATCTCCGCCGATGATTCGGCGCAACCTGGGGTGGTGCGAAg                                                                                                                                   >  1:866389/1‑239 (MQ=255)
                                                                                                           ttaggactctgttttATCGGCTGGTTGAACCTGAACGGGGACTCCAGTTTGTACTCCCTGGCAGGCCATGTTGGCGGTATCGGTGTCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGc                                                                                                                               <  1:269594/226‑1 (MQ=255)
                                                                                                           ctcaaaTATGTTGTTATGCCCTGGATGAACGTTGAGGGTGTCTACATTTGTGACGCCTTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACAGCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGc                                                                                                                               <  1:889942/235‑1 (MQ=255)
                                                                                                                      tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcc                                                                                                                                                                                                                                          >  2:264575/1‑121 (MQ=255)
                                                                                                                      tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTg                                                                                                                    >  2:905931/1‑239 (MQ=255)
                                                                                                                      tgTTATCGCCTGTATGAACCTGAACGGGGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcc                                                                                                                                                                                                                                          <  1:264575/121‑1 (MQ=255)
                                                                                                                         tttCGCCTGTATGGACCGTAAAGGGGAGTTCAGTTCTGACGGCCGGGCAGCGTAGGTTGGAGGTAGCGGGATCGCCCCGTGGGAAAAGCGCGGTGACGAATGCGCCCTGTTTGAAGCCACGTGCGGTACTGCGCCGAAATAGGCCGGGCAGGAGAAAGTAATTCCTGGCGCTATTATTCTCTCCGCTGAGATGATGATGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGtt                                                                                                                  <  1:501560/236‑1 (MQ=255)
                                                                                                                         tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                        >  2:656366/1‑168 (MQ=255)
                                                                                                                         tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                                 >  2:269056/1‑239 (MQ=255)
                                                                                                                         tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                                                                 <  2:1097141/239‑1 (MQ=255)
                                                                                                                         tATCGCCTGTATGAACCTGAACGGGGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                        <  1:656366/168‑1 (MQ=255)
                                                                                                                             ccttGTTTGACCCTGAAGGTGTCTCCATTTAATGAGCCCGGGGAGCGCGGGTTGGGGGTTCCGTACCCGCCACTGGTGCGAACCATGGTGGCGAATGTGACCTGTTTGAAGCAACCCACGGTACTGCGCGGAAATATGCCGGTCAGGACAAAGTAAATCCTGGATCTATTTTTCTCTCTGCTGAGATGATGCTGCGCGACATGGGTTTGACCGAAGCGGCTGCCTTAATTGTTAAAggg                                                                                                            <  1:1295696/236‑2 (MQ=255)
                                                                                                                               cTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                                                                            <  2:477447/238‑1 (MQ=255)
                                                                                                                                tGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGCTGACTTAAtt                                                                                                                     >  2:269594/1‑228 (MQ=255)
                                                                                                                                  ttttaaCCTGAATGGTTACTAAGTTTGCGACGCCCTGGCAGGGCAGGTGGGCGGGCTCGGTATCGCCCCTGGTGAAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTTCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                                                                        <  1:704470/235‑1 (MQ=255)
                                                                                                                                  tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAACGGCTGACTTAATTGTTAAAGGTATgg                                                                                                        >  2:493024/1‑239 (MQ=255)
                                                                                                                                  tATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                                                                        >  1:746433/1‑239 (MQ=255)
                                                                                                                                       aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCCCTATTATTCTCTCCGCTGAGATGATGCTGCCACGCATGGGTTGGCACGAAGCGGCCGACTTCATTGTTAAAGggaaggaaaga                                                                                                   >  1:521419/1‑236 (MQ=255)
                                                                                                                                         cTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGATGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTAAAGGTATGGAAGgcgc                                                                                                 >  2:1036142/1‑239 (MQ=255)
                                                                                                                                                          ttCTGACGCCCTGGCAGCGGAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCGCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                               <  2:591183/239‑1 (MQ=255)
                                                                                                                                                          ttCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                               <  2:111905/239‑1 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                            >  2:376571/1‑237 (MQ=255)
                                                                                                                                                                       ggCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCAACATGGGTTGGACCGAAACGGCTGACTTAATTGTTTAAGGTATGGAAGGCGCAATCAAAGCGAAAACCGGAACCTACGaat                                                                   >  1:211421/1‑237 (MQ=255)
                                                                                                                                                                             gcAGTTTGGCGTTTTCTTTTTCCCCCCTGTTGCAAACCTCGGGTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGc                                                             <  1:493024/239‑1 (MQ=255)
                                                                                                                                                                              cAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATAACTTAGAGCg                                                            >  2:920146/1‑239 (MQ=255)
                                                                                                                                                                                        tgttttggttgcgacccTGGTGCAAACACCGGGCCCGAAGGCTCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGAGGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                   <  1:347148/225‑1 (MQ=37)
                                                                                                                                                                                         tatcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTTTTAACGGGATGGAAGGGCACAACAACGCGAAAAACGGAACATATAACTTCGAGAGTcgggaggggg                                                 >  1:435922/1‑230 (MQ=255)
                                                                                                                                                                                           tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGGTGgc                                               >  2:236692/1‑239 (MQ=255)
                                                                                                                                                                                           tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATGGTTAAAGGCATGGAAGGGCCAATCAACGCGAAAACCGTAACCTATGACTTCTAGCGTCTGatggatggg                                               >  2:976985/1‑238 (MQ=255)
                                                                                                                                                                                                    gCCCCTGGTTCAAATATCGGTGACGAAGGGGCCATGTTTGAAGCCCCCACGGGTCGGCGCCGAAATAGGCCGGTCAGGAAAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAActg                                      <  1:750988/239‑1 (MQ=37)
                                                                                                                                                                                                          ggggCAAACCTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGaa                                 <  2:771150/236‑1 (MQ=255)
                                                                                                                                                                                                          ggTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCCCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGaa                                 <  2:872444/239‑1 (MQ=255)
                                                                                                                                                                                                              cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGCAAATCCTGGCTCTATTATTCTCCCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTATTTGTTAAAGGGATGGAAAGCGAAATCAACCCGGAAACCGAAACCTTTGACTCCTTGAGTCTAATAGATGGGCCAAAAACGCTGaaaatt                             >  1:601428/1‑236 (MQ=255)
                                                                                                                                                                                                               aaaCTTCGGTGCAGAATGGCCCCTGTTTGAAGCCACCAACGGTACTTGGCCGAAATAAGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCCCTCCGCTGAGAAGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAGCTAGGACTGCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                            <  1:348225/239‑1 (MQ=25)
                                                                                                                                                                                                                aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCCCTATTATTCTCTCCGCTGAGATGAGGCTGCGCCACATGGGTTGGTCCGAAGAGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACCCGAGAGCAGTACCCTAGATCATCACGGGTCTCATGGGTGGCGCAAAAAGTCTGAAATGTTc                           >  1:364539/1‑239 (MQ=255)
                                                                                                                                                                                                                aaCATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGGTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGAAAACCCTAACCCCTGACTACGGGCGTCTGTAGCATGGCGCTAAACGGTTAAAATGTTc                           >  1:62282/1‑239 (MQ=255)
                                                                                                                                                                                                                 aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTACGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                          <  2:211421/239‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTGACGAATGCGCCCTGTTTGCAACCACCCACGGTACTGCGCCGTAATATGCCGGCCAGGGCAAAGTAAATCCTGGGTCTATTCTTCTCTACGCTGCGATGATATTGGGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGGTAAAGGTAGGGAAGGGGCAATCAACGCAAAAACCGAACCCAATGAAGTCCAGCGTCTGACGAATGGCACaaaaaagcttaacattccagaat                      >  1:1243574/1‑220 (MQ=255)
                                                                                                                                                                                                                     cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGa                                                                                                                                     <  2:707422/129‑1 (MQ=255)
                                                                                                                                                                                                                     cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGa                                                                                                                                     >  1:707422/1‑129 (MQ=255)
                                                                                                                                                                                                                          aCGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGCTGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCCAAGCGAAAAACGTAACCTATGACTTCGAGCGACTGAGGGATGGGGCTAAAATGCTGAAAGGTTCAGAGTTTggg                 >  1:649206/1‑238 (MQ=255)
                                                                                                                                                                                                                               tGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTCAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGCCCGAAACGGCTGGATTAATTGTTACAGGTCTGGAAGGCGCAATCAACGCAAAACCCGTAACCTATGACTTCGGGCGTCTGATGGATCTCGCTCAACTCTTGAATTGTTCCGCGTTTGGTGAcgc            >  1:63320/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gcCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg           <  2:746433/237‑1 (MQ=255)
                                                                                                                                                                                                                                         ttGAAGCCCCCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCGATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa  <  2:841921/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
TCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAA  >  NC_000913/1195895‑1196364

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