Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,305 C→T 79.2% R395C (CGT→TGT)  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,3050CT79.2% 128.8 / 26.1 68R395C (CGT→TGT) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (22/28);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAG  >  NC_000913/1196078‑1196541
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
gtatcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                                     <  1:953988/238‑1 (MQ=255)
        atcgCCCCTGTTGCAAACAGCGGTGACGAAGTCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTa                                                                                                                                                                                                                             <  1:1013266/238‑1 (MQ=255)
                tGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGa                                                                                                                                                                                                                    >  1:394764/1‑239 (MQ=255)
                 tggggcaaCTTCGGGTACGAATGCGCCTTTTTTAAAGCCCCCCACGGTACTGCGACGAAATATGCCGGTAAGGACAAAGTAAATCCTGGTCCTATTATTCCCTCCGGTGGGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGaa                                                                                                                                                                                                                   <  2:394764/233‑1 (MQ=255)
                 ggTGCAAACATGGGTTACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGaa                                                                                                                                                                                                                   <  2:179037/239‑1 (MQ=255)
                 ggTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaa                                                                                                                                                                                                                   <  1:574421/239‑1 (MQ=255)
                       aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGAAACCTATGACTTCAAGGGTCTGATGAATAGCCCTAACATGTGGAAATGTTc                                                                                                                                                                                                             >  1:1087598/1‑239 (MQ=38)
                             ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt                                                                                                                                                                                                       <  1:297871/239‑1 (MQ=255)
                                    aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCATGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCCTAACCCATGACTTCGAGAGTCTGATGGATGGCGCAAAACAGTGGAAATGTTAAGAGTTGGGGGag                                                                                                                                                                                                >  2:549289/1‑238 (MQ=37)
                                                 tagaGCCCCCCACTGTACTGCGCCGAAATACTCCGGGCAGGACAAAGTAAATCCTGGCTTTATTATTCGCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAATGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCAACGaaa                                                                                                                                                                                   <  2:633200/236‑1 (MQ=255)
                                                 tGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                   <  1:624386/239‑1 (MQ=255)
                                                      ccccccaaGGTACAGCCGCCAAAATTGCCGGTTAGGACAAATTAAACCCTGGCACTATTATTCTCTCCGCTGAGATGATGCTGCGCCCAATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAACCAACGCGAAAACCGTACCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                              <  2:640257/236‑1 (MQ=255)
                                                         cccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGATTGGGGCCGCGATCATCGAAAACATGTaa                                                                                                                                                                           >  1:817489/1‑239 (MQ=255)
                                                                 aCTGCGCCGAAATAAGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTa                                                                                                                                                                   <  1:379223/238‑1 (MQ=255)
                                                                   tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                                                                                                                                                                                                                           <  2:40344/181‑1 (MQ=255)
                                                                   tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGAGGGATGGCGCTaaa                                                                                                                                                                                                                           >  1:40344/1‑181 (MQ=255)
                                                                        ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAAATGTAAta                                                                                                                                                                         >  1:516886/1‑226 (MQ=255)
                                                                            atttttccGGTCAGGAAAAAGTAAATCCTGGCTCTATTATTCTCTCCGTTGAGATGATGCAGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAt                                                                                                                                                        <  2:210599/233‑1 (MQ=255)
                                                                            attttGCGGGGCAGGACAAAGTAAATCCGTGCTGTATTATTCTCTCCGTTGAGAGGATGCTGCGGCACATGGGTTTGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGGGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAt                                                                                                                                                        <  2:1087598/235‑1 (MQ=255)
                                                                                tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                   <  1:865942/209‑1 (MQ=255)
                                                                                tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                                                                                                                                   >  2:865942/1‑209 (MQ=255)
                                                                                     gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatggt                                                                                                                                                                                                                                     >  2:117162/1‑153 (MQ=255)
                                                                                     gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                                     <  1:117162/154‑1 (MQ=255)
                                                                                       cAGGCCAAAGTAAATCCTGGTTTTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACgg                                                                                                                                              <  2:462113/238‑1 (MQ=255)
                                                                                              aaGTAAATCCTGGGTTTGTTATTTTCTCCGCTGAGATGATGGTGCGCCACATGGGTTGGACCGAAGCTGCTGAATTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                                                                                                                      <  2:784851/239‑1 (MQ=255)
                                                                                               aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGGAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                                     >  1:288560/1‑144 (MQ=255)
                                                                                               aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGGAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgga                                                                                                                                                                                                                                     <  2:288560/144‑1 (MQ=255)
                                                                                                    aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCCTCGAAAAAATGTAATGCCGGAGTTTGTTAAATTTATTAACGGAAgggaaaacgccc                                                                                                                                >  2:283471/1‑228 (MQ=21)
                                                                                                      tcTGGCTCTATTATTCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                                                              <  1:32926/238‑1 (MQ=255)
                                                                                                         ggCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCTGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTATCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTCAAATGTTAAGAGCTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttg                                                                                                                           >  1:56267/1‑239 (MQ=255)
                                                                                                                ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                                                                                                                                              <  1:1099739/182‑1 (MQ=255)
                                                                                                                    tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTACTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAg                                                                                                                <  1:4274/239‑1 (MQ=255)
                                                                                                                     ttctcCCGTGGGATGATGCGTCGCCAAAGGGGTTGGACCGAAGGGGCTGACTTAATGGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCGATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTGTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                               <  2:1073750/238‑1 (MQ=255)
                                                                                                                     catgccGCTTAGATGATGCTAGCCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAAACCGTTACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTGTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                                                               <  1:598112/235‑1 (MQ=255)
                                                                                                                        tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAg                                                                                                                >  1:1054569/1‑235 (MQ=255)
                                                                                                                        tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgct                                                                                                             >  1:1077275/1‑238 (MQ=255)
                                                                                                                        tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgct                                                                                                             >  1:930341/1‑238 (MQ=255)
                                                                                                                        tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCACCCGTTGTTTTTTGTTAg                                                                                                                >  2:330452/1‑235 (MQ=255)
                                                                                                                        tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgct                                                                                                             >  1:429174/1‑238 (MQ=255)
                                                                                                                          cGCTGGGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATAGTTAAAGGTATGGAAGGCGCAATCAACGCGACACCCGTAACCTATGACTTCAAGGTTATGCTGGAGAGCGCTAAACTGCTGAAATGTTCAGGCTTTAGTTACACGCTCATCGAAAACAGTCATGACCTTAGTTTTGTAACTTTATTAACGGGAGCGAAACGATCCCGTTGTTTTTTGTTAGgcggta                                                                                                          >  2:548834/1‑235 (MQ=255)
                                                                                                                          cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCACAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAACTTTATTAACGGGAGCGAAAAGCTCCCTTTGTTTTTTTTTAGgctgtt                                                                                                          >  2:422876/1‑237 (MQ=255)
                                                                                                                                gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                                                     <  1:836829/238‑1 (MQ=255)
                                                                                                                                    atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTa                                                                                                                                                                            <  1:868789/164‑1 (MQ=255)
                                                                                                                                    atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTa                                                                                                                                                                            >  2:868789/1‑164 (MQ=255)
                                                                                                                                                          ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                                                                          >  1:523462/1‑239 (MQ=255)
                                                                                                                                                             aaGAGGCTGACTTAATTGTTAAAGGTATGGAAGGCTCAAGCAACGCGAAAAAGGTAACCTCTGACTTCGAGTGGCTGCTGGAAGAAGCTAAACGGCTGAAATGTACAGAAGTTGGTGACGCATGCATCGAAACCCGGAAAGGACGTAGGTTGTTAAGATTATTAACGGGGCGGTAACGCACCCGTCGCTTTTTGGTGAGAGGCTACAGGGCAACAAAATTTTATAAAAAAAAGtagtta                                                                       >  2:1024880/1‑234 (MQ=255)
                                                                                                                                                                  gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGCCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                                                   >  1:554081/1‑238 (MQ=255)
                                                                                                                                                                     gACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTGCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                                                               <  1:83715/239‑1 (MQ=255)
                                                                                                                                                                       cTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGCAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTCGTTTTTTGTTAGGCGGCTAACGGTTATAAAAATTTTATCAAAAAAAGTTACAAAAACCCCCCg                                                             >  2:321136/1‑239 (MQ=37)
                                                                                                                                                                          aaTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTTTCAAAACCCCTCGGTa                                                          >  1:327714/1‑239 (MQ=255)
                                                                                                                                                                                tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGCAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAATAGTTATCAAAACACCTCGGTAGTTTTg                                                    >  1:585896/1‑239 (MQ=255)
                                                                                                                                                                                     gTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATtta                                                                                                                                                     <  1:662076/137‑1 (MQ=255)
                                                                                                                                                                                            aGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTgg                                         <  2:327714/238‑1 (MQ=255)
                                                                                                                                                                                                gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGt                                    <  2:1054569/239‑1 (MQ=255)
                                                                                                                                                                                                       accggAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGGTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtgg                             <  1:432921/235‑1 (MQ=255)
                                                                                                                                                                                                       aCGCGAAAACCGTAACCTATTACTTCGATTGGCTGATGGATGGCGTTAAACTGCTGAAATTTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAggtgg                             <  1:1006768/239‑1 (MQ=255)
                                                                                                                                                                                                                       cTATGGCTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGt             <  1:254793/239‑1 (MQ=255)
                                                                                                                                                                                                                         aTGACTTCGATCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGtct           <  1:825888/239‑1 (MQ=255)
                                                                                                                                                                                                                         aTGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTCGTTTTGGGGTAGGCGTGCCGGTCAGGTGGTAGTTCTACTACTAGtc            >  1:71925/1‑238 (MQ=255)
                                                                                                                                                                                                                         aTGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGtct           <  1:72691/239‑1 (MQ=255)
                                                                                                                                                                                                                          tggcTTCGGGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGtctc          <  1:121490/236‑1 (MQ=255)
                                                                                                                                                                                                                          tGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTTAAATGTTCAGATTTTGGTGACGCGATCATCGAAAACATGCAATGCCGTAGTTTGTTAAATTTATTACCGACAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATGAACATTTTATCAAAACAACATATCAAACCACCCCGGTAGGTTTGGGTTAGGCCGACCGGCAAGGTGTAAGTTataatacaaatccc          >  2:997275/1‑232 (MQ=255)
                                                                                                                                                                                                                             cTTCGAGGGTCGGATGGATGGCGCTAAACTGGTGAAATGTTCAGAGTTTGGTGACGCTATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCca       <  1:581160/239‑1 (MQ=255)
                                                                                                                                                                                                                                cGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgt                                                                                                                             <  1:225442/118‑1 (MQ=255)
                                                                                                                                                                                                                                cGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGCGGTAGTTCTACTACTAGTCTccaccat    >  1:44742/1‑234 (MQ=255)
                                                                                                                                                                                                                                cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTCGTTTTGGGGTAGGCTGGCCGGTCAGGTGGAAGTTTACCTACTAGTCCCCcacaa    >  2:121490/1‑238 (MQ=255)
                                                                                                                                                                                                                                agAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTACCGGGAGCGTAACGCTCCCgt                                                                                                                             >  2:225442/2‑118 (MQ=255)
                                                                                                                                                                                                                                  aGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTTTCAAAAAAAGTTATCAAAACCCCTCGGGAGTTTTGGGGTAGGCTGGCCGGTCAGGGGGTAGTTTTACCACTAGTTCCCCAAATAg  >  2:72691/1‑239 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
GTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGGTCAGGTGGTAGTTCTACTACTAGTCTCCCACATAG  >  NC_000913/1196078‑1196541

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