Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 20.4% 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.4% 131.1 / 34.1 73A337T (GCC→ACC) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base A (8/7);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

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