Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 28.2% 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,1310GA28.2% 108.9 / 49.1 75A337T (GCC→ACC) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (11/10);  total (38/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 = 1.00e+00

ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCG  >  NC_000913/1195894‑1196362
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
aTCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCCTTCCTGCAACAGATCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAACAATCGGGGACGAATGCGCCCCGTTTGAAg                                                                                                                                                                                                                                          >  1:210263/1‑238 (MQ=255)
 ttgccGGTGGCCCGTGGCTGAAATTTAAAAACCCGAACCCTGGCAAAGAGATCGGCATTAAAGGCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcc                                                                                                                                                                                                                                        <  2:709672/235‑1 (MQ=255)
  cGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGcc                                                                                                                                                                                                                                        <  1:806775/238‑1 (MQ=255)
      ggTGGCCCGGGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGccaccca                                                                                                                                                                                                                                   <  1:210523/239‑1 (MQ=255)
          gCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGGGCCGAATGCGCCCTGTTTGAAGCCACCCACgga                                                                                                                                                                                                                               >  1:380273/1‑238 (MQ=255)
          gCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACACCGGGGACGAATGGGCCCGGTTTGAAACCACCCACgga                                                                                                                                                                                                                               >  2:848444/1‑238 (MQ=255)
            ccGGGGCTGAAAGTTAAAAACCCGAACCCTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                                                                                                                                                              <  1:243815/238‑1 (MQ=255)
               tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGCCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGTCTGAATATGATGTTATCGCCTGTATGAACCTGAACGTGGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCTCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGCAGCCACCCACGGTACTgc                                                                                                                                                                                                                          >  1:669351/1‑239 (MQ=255)
               tGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATCGGCCCCGTTTGAAGCCACCCACGGGACTgc                                                                                                                                                                                                                          >  2:210523/1‑239 (MQ=255)
                  cTGAAAGTTAAAAACCCGACCTCTGGCAAAGGGTTCGTCATTAAAGGCGGTATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTgcgc                                                                                                                                                                                                                        <  2:3645/238‑1 (MQ=255)
                          tAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGGCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATg                                                                                                                                                                                                               <  1:661840/239‑1 (MQ=255)
                          tAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATg                                                                                                                                                                                                               <  1:228390/239‑1 (MQ=255)
                               accccaacctcggCAAAGAGACGTTCATTAAAGACGGTGTTGCTGATGCATTCCTGCAACGGATCCTGGTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTTAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGt                                                                                                                                                                                                          <  2:669351/228‑1 (MQ=255)
                                 ccGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGTGCAGGTTGGCGGTATCGGTATCACCCCTGGTGCAACCCTCGGTGACCACTGACCCCTGTTTTAAGCCACCCCCGGTAccccccccacaaatcccgtgcg                                                                                                                                                                                                        >  2:473124/1‑218 (MQ=255)
                                    aaCACTGGCAAAGAGAGCGGTATTAAAGGCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGGATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGa                                                                                                                                                                                                     <  1:275963/239‑1 (MQ=255)
                                         tGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGCTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAACCATCGGTGACCAATGCGCCCTGTTTGAAGCACCCCCCGGGACTCCGCCGAAAAATCCCGGCCCGGGCAAAg                                                                                                                                                                                                >  2:389193/1‑239 (MQ=255)
                                                  gaTCGTCAGTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATcc                                                                                                                                                                                         <  1:764216/237‑1 (MQ=255)
                                                   aTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGAACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAACCTg                                                                                                                                                                                       >  1:369636/1‑238 (MQ=255)
                                                    tCGTCTTTTAAGGCGTGATTGGTGAATCATTCCTGCAACGGGTCCGGCTGCGTCCGGCTGAATATGATGTTATCTCCGGTAAGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  2:210263/238‑1 (MQ=255)
                                                    tCGTCATTAAAGACTTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:389193/238‑1 (MQ=255)
                                                     cGTCATTAAGGGCGGTATGGCTGATGCAGTCCTGCAACGGGTTCTGCTGCGTCCGGCTGAATAGGAGGTTATCGCCTGTATGAACATGAACGGTGACTACATTTCGGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATGGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  2:565175/237‑1 (MQ=255)
                                                           tAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTCCATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAACCATCGGTGACGAATGCGCCCCGTTTGAACCCACCCCCGGTCCTGCGCCGAAAAACGCCGGTCCGGACAAAGTAAAACCTGGccccctta                                                                                                                                                                              >  2:418350/1‑232 (MQ=255)
                                                                  gTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAATAAATCCTGGCTCATTTTTtctcc                                                                                                                                                                        >  2:143690/1‑237 (MQ=255)
                                                                     atttgttgtgaaTTCCTGCCACAGACCCGGCTCTGTCCGGCTGAATATGAGTTTATCGCCTTTATGAACCTGAACGGTGACTACCTTTCTGACGCCCTGGCAGCTCAGGGTGGCGGTATCGGTATCGCCCCTGGTTCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    <  2:138336/228‑1 (MQ=255)
                                                                     aTTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGAACTGCGCCGAAATAGGCCGGTCAGGACAAAGTAAATCCTTGCCCCAATATTCTCTcccc                                                                                                                                                                    >  1:193747/1‑237 (MQ=255)
                                                                     aTTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGCAATATGCCGGTCAGGACAAAGTAAAACCCGGCCCTATTATTCCCCCCGc                                                                                                                                                                    >  2:726084/1‑239 (MQ=255)
                                                                      ttGCTGATGCATTCCTTCAACATGTCCTGCGGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCt                                                                                                                                                                   <  2:306451/239‑1 (MQ=255)
                                                                            aTGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:672801/214‑1 (MQ=255)
                                                                               ctttCCTTCAACACTTCCTGTTTCTTCCGGTTGAATCTGATGTTATCGCCGGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGatg‑atg                                                                                                                                                         <  2:482936/237‑1 (MQ=255)
                                                                                      gCAACAGATCCTGCTGCTTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGAGCCCTGTTTGAAACCACCCACGGTACTCCGCCGAAAAATGCCGGTCAGGACAAACTAAATCCTGCCCCTTTTCTTCTCTCCCCTGGGATG‑Aggcgtgacc                                                                                                                                                  >  2:544271/1‑233 (MQ=255)
                                                                                      gCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACCTCGGGGCCGAATGCGCCCTGCTTTAAGCCACCCACCGTACTGCCCCGAAATATGCCGGTCAGGCCAAAGAAAATCCTTGCTCTATTATtcctgtcccttttaaccacctcacg                                                                                                                                                    >  2:672801/1‑215 (MQ=255)
                                                                                        aaCAGTTCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCcac                                                                                                                                                <  1:544271/239‑1 (MQ=255)
                                                                                         aCAGATCCTGCTGGGTCCGGCTTAATATGGTGTTTTCGCCTGTATGAACCTTAACGGTGACTTCATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCcac                                                                                                                                                <  2:411974/238‑1 (MQ=255)
                                                                                              tCCTGCGGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCcac                                                                                                                                                <  1:236860/233‑1 (MQ=255)
                                                                                                        tCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAAGGCGCCCTGTTTGAAGCCACCCTCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTCAATCCTGGCTCTATTATTCTCTCCCCTGAGATG‑GTGCAGGGCcaacgggggtgggccag                                                                                                                                  >  2:236860/1‑222 (MQ=255)
                                                                                                                   atgatgTTATCCCCTGTATGAACCGGAACGGTGACTTCATTGCTGGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACAGCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                                                                     <  1:418350/239‑1 (MQ=255)
                                                                                                                    tggatgtaccGCCTGTATGAACCTGAACGGTGACTCCATTTCTTACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                                    <  1:848444/230‑1 (MQ=255)
                                                                                                                       tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTTGTGCAAACATCGGTGACGCATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAACCTGGCTGTTTTATTCTCTCCCCTG‑GAGGAAGGCAGCGCCACATGGGTTGAACCGAACCGGCGGACTTaaagt                                                                                                                 >  2:2933/1‑236 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTCTTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAACGGCTGACTTAATTGTTa                                                                                                              >  1:33432/1‑239 (MQ=255)
                                                                                                                           aTCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCGGACTTAATTGTTaa                                                                                                             >  1:673300/1‑239 (MQ=255)
                                                                                                                               ccTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                                                                         >  1:321565/1‑239 (MQ=255)
                                                                                                                                        aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    <  1:794376/172‑1 (MQ=255)
                                                                                                                                        aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    >  2:794376/1‑172 (MQ=255)
                                                                                                                                             aaaggTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTATCCGCTGAGATG‑ATGCTGCGCCACCTGGGTTGGACCGAAGGGGCTGCCTTAATTGTTAAAGGAATGGAAGGCGCAAt                                                                                           >  1:451561/4‑239 (MQ=255)
                                                                                                                                             aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGCAATATGCCGGCCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCCGGGCCACATCGGTTGGACCGAAGCGGCCGACTTAATTGTTAAAGGGATGGAAGGCGAAAt                                                                                           >  2:9376/1‑239 (MQ=255)
                                                                                                                                             aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAAt                                                                                           >  1:99893/1‑239 (MQ=255)
                                                                                                                                                  tGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGGTGCGCCACATGGGTTGGACCGAAGAGGCTGACTTAATTGTTAAAGGGATGGAAGGGGGAAACaaag                                                                                      >  2:910558/1‑237 (MQ=255)
                                                                                                                                                  tGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTAACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTTTTAAAGGAAAGGAAGGGGCAAACaaag                                                                                      >  2:357750/1‑237 (MQ=255)
                                                                                                                                                     ctccATGTCTGGCGCCCTGGGAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                                   <  2:99893/236‑1 (MQ=255)
                                                                                                                                                     cTACCTTTCCGCCGCCCTGACAGCTCAGTGTGGCCGTTTCGGTTACGCCCCGGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTTCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                                   <  2:451561/239‑1 (MQ=255)
                                                                                                                                                      tccATTTCTGACGCCCTGGCAGCGCGGGTTGGCGTTTTCGGTATCGCCCCTGGGGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                                                                  <  2:673919/237‑1 (MQ=255)
                                                                                                                                                          gttttagcGCCCTGGCAGCGCAGGTTGGTGTTTCTGGTTTCGCCCCTGGTGCAAAACTTGGTGACGAATGCGCCTTGTTTGAAGCCACCCACGGTACTGCGCCGAAATAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTT‑AAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                             <  2:380273/232‑1 (MQ=255)
                                                                                                                                                           tgcTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                             <  1:420529/237‑1 (MQ=255)
                                                                                                                                                              tGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtct                                                                                                                                                                          <  1:430019/144‑1 (MQ=255)
                                                                                                                                                              tGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCCGGCTCTATTATtct                                                                                                                                                                          >  2:430019/1‑144 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAACTCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAGGGAAGGCGCAATAAACCCGAAAACCGTAAc                                                                         >  1:359320/1‑238 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      >  2:169778/1‑130 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:169778/130‑1 (MQ=255)
                                                                                                                                                                            cgcAGGGTGGCGGTATCGGTATCGCCCCTGGTGCAAACCTCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAAGATGCCGGTCAGGACAAAGTAAATCCTGGCACAATTATTCTCTCCGACGAGCTC‑ATGCTGCGCCACAGGGGTTGGACCGAAGCGGCTGTCCTACTTGTTAAAGGAATGGAAGACACAATCAACCCAGAAACCCTAAccccagacatggaa                                                            >  2:623744/1‑227 (MQ=255)
                                                                                                                                                                                            cggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                            <  1:548390/239‑1 (MQ=255)
                                                                                                                                                                                                atcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTa                                         <  1:910558/238‑1 (MQ=255)
                                                                                                                                                                                                 tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGCAGCGGCTGACTTAATTGTTCAAGGTATGGAAGGCGCACTCAACGCGCAAACCGACACCAATGACTTCGGGGGTCTGATGGAAGGGGCTcaa                                       >  2:831038/1‑239 (MQ=255)
                                                                                                                                                                                                 tcgCCCCGGGTGCAAACATCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                                       <  1:357750/239‑1 (MQ=255)
                                                                                                                                                                                                   gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAAc                                      >  1:913729/1‑238 (MQ=255)
                                                                                                                                                                                                   gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGACCGTCTGATGACTGGCGCTaaccg                                     >  1:35998/1‑236 (MQ=255)
                                                                                                                                                                                                         ggTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAACCCTGGCTCTATTATTCTCTCCGCTGCGATG‑ATGCTGCGCCACATGGGTCGGCCCGAAGCGGCTGCCTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAACCCGTCACCTATGACTTCGAGGGGCGGAAGGAAGGGGCTAAAAAGCTGaa                               >  2:434045/1‑239 (MQ=255)
                                                                                                                                                                                                              aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                          <  1:448116/239‑1 (MQ=255)
                                                                                                                                                                                                                aCATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                        <  1:726084/239‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGCAGGCTGGCGCTAAACTGCTGAAAAGGTCAGAGt                    >  2:320376/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGGAAACCGTAACCTATGACTTCGAGTGGCTGATGGATGGGGCCAAACTACTAAAATGTTCAGGGtt                   >  2:448116/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTATCTCCGCCGAGCTG‑ATCCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGGAATCAACGCGAAAAACGTAAACTATAACTTCGAGAGTGTATAGTATGGCGCTAAAATGCTAAAAtttttaaag                     >  2:18692/1‑232 (MQ=255)
                                                                                                                                                                                                                            aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGTGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGGTGGAGGGCGCTAAACTGCTGAAATTTTCAAAGTTTggggcc            >  2:15854/1‑235 (MQ=255)
                                                                                                                                                                                                                              tGCGCCCTGTTTGAAACCAACCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTAGACCGTCTGATGGATGCCGCTAAACTGATGAAATGTTCAGAGTTTGGtacggc          >  1:728157/1‑234 (MQ=255)
                                                                                                                                                                                                                               gggcCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg         <  1:143690/237‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGCAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTAAGAGTTTGGTGACGCGatca     >  1:164224/1‑239 (MQ=255)
                                                                                                                                                                                                                                      gTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCg  <  1:15854/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATG‑ATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCG  >  NC_000913/1195894‑1196362

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