Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,131 G→A 28.5% 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.5% 104.3 / 48.4 71A337T (GCC→ACC) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base A (10/10);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCG  >  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)
                                    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)
                                                                               ctttCCTTCAACACTTCCTGTTTCTTCCGGTTGAATCTGATGTTATCGCCGGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGatgatg                                                                                                                                                         <  2:482936/237‑1 (MQ=255)
                                                                                      gCAACAGATCCTGCTGCTTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGAGCCCTGTTTGAAACCACCCACGGTACTCCGCCGAAAAATGCCGGTCAGGACAAACTAAATCCTGCCCCTTTTCTTCTCTCCCCTGGGATGAggcgtgacc                                                                                                                                                  >  2:544271/1‑233 (MQ=255)
                                                                                        aaCAGTTCCTGCTGCGGCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                                                                                                <  1:544271/239‑1 (MQ=255)
                                                                                         aCAGATCCTGCTGGGTCCGGCTTAATATGGTGTTTTCGCCTGTATGAACCTTAACGGTGACTTCATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                                                                                                <  2:411974/238‑1 (MQ=255)
                                                                                              tCCTGCGGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcac                                                                                                                                                <  1:236860/233‑1 (MQ=255)
                                                                                                                   atgatgTTATCCCCTGTATGAACCGGAACGGTGACTTCATTGCTGGCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACAGCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                                                                     <  1:418350/239‑1 (MQ=255)
                                                                                                                    tggatgtaccGCCTGTATGAACCTGAACGGTGACTCCATTTCTTACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                                                                    <  1:848444/230‑1 (MQ=255)
                                                                                                                       tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTTGTGCAAACATCGGTGACGCATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAACCTGGCTGTTTTATTCTCTCCCCTGGAGGAAGGCAGCGCCACATGGGTTGAACCGAACCGGCGGACTTaaagt                                                                                                                 >  2:2933/1‑236 (MQ=255)
                                                                                                                          tATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTCTTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAACGGCTGACTTAATTGTTa                                                                                                              >  1:33432/1‑239 (MQ=255)
                                                                                                                           aTCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCGGACTTAATTGTTaa                                                                                                             >  1:673300/1‑239 (MQ=255)
                                                                                                                               ccTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                                                                         >  1:321565/1‑239 (MQ=255)
                                                                                                                                        aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    >  2:794376/1‑172 (MQ=255)
                                                                                                                                        aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGc                                                                                                                                                                    <  1:794376/172‑1 (MQ=255)
                                                                                                                                             aaaggTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTATCCGCTGAGATGATGCTGCGCCACCTGGGTTGGACCGAAGGGGCTGCCTTAATTGTTAAAGGAATGGAAGGCGCAAt                                                                                           >  1:451561/4‑239 (MQ=255)
                                                                                                                                             aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGCAATATGCCGGCCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCCGGGCCACATCGGTTGGACCGAAGCGGCCGACTTAATTGTTAAAGGGATGGAAGGCGAAAt                                                                                           >  2:9376/1‑239 (MQ=255)
                                                                                                                                             aaCGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAAt                                                                                           >  1:99893/1‑239 (MQ=255)
                                                                                                                                                  tGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGGTGCGCCACATGGGTTGGACCGAAGAGGCTGACTTAATTGTTAAAGGGATGGAAGGGGGAAACaaag                                                                                      >  2:910558/1‑237 (MQ=255)
                                                                                                                                                  tGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTAACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGCCTTAATTTTTAAAGGAAAGGAAGGGGCAAACaaag                                                                                      >  2:357750/1‑237 (MQ=255)
                                                                                                                                                     ctccATGTCTGGCGCCCTGGGAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                                   <  2:99893/236‑1 (MQ=255)
                                                                                                                                                     cTACCTTTCCGCCGCCCTGACAGCTCAGTGTGGCCGTTTCGGTTACGCCCCGGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTTCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGa                                                                                   <  2:451561/239‑1 (MQ=255)
                                                                                                                                                      tccATTTCTGACGCCCTGGCAGCGCGGGTTGGCGTTTTCGGTATCGCCCCTGGGGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                                                                  <  2:673919/237‑1 (MQ=255)
                                                                                                                                                          gttttagcGCCCTGGCAGCGCAGGTTGGTGTTTCTGGTTTCGCCCCTGGTGCAAAACTTGGTGACGAATGCGCCTTGTTTGAAGCCACCCACGGTACTGCGCCGAAATAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                             <  2:380273/232‑1 (MQ=255)
                                                                                                                                                           tgcTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                                             <  1:420529/237‑1 (MQ=255)
                                                                                                                                                              tGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtct                                                                                                                                                                          <  1:430019/144‑1 (MQ=255)
                                                                                                                                                              tGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCCGGCTCTATTATtct                                                                                                                                                                          >  2:430019/1‑144 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAACTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAGGGAAGGCGCAATAAACCCGAAAACCGTAAc                                                                         >  1:359320/1‑238 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      >  2:169778/1‑130 (MQ=255)
                                                                                                                                                                aCGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTgg                                                                                                                                                                                      <  1:169778/130‑1 (MQ=255)
                                                                                                                                                                                            cggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGgcg                                            <  1:548390/239‑1 (MQ=255)
                                                                                                                                                                                                atcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTa                                         <  1:910558/238‑1 (MQ=255)
                                                                                                                                                                                                 tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGCAGCGGCTGACTTAATTGTTCAAGGTATGGAAGGCGCACTCAACGCGCAAACCGACACCAATGACTTCGGGGGTCTGATGGAAGGGGCTcaa                                       >  2:831038/1‑239 (MQ=255)
                                                                                                                                                                                                 tcgCCCCGGGTGCAAACATCGGTGAAGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                                       <  1:357750/239‑1 (MQ=255)
                                                                                                                                                                                                   gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAAc                                      >  1:913729/1‑238 (MQ=255)
                                                                                                                                                                                                   gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGACCGTCTGATGACTGGCGCTaaccg                                     >  1:35998/1‑236 (MQ=255)
                                                                                                                                                                                                         ggTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAACCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTGCGCCACATGGGTCGGCCCGAAGCGGCTGCCTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAACCCGTCACCTATGACTTCGAGGGGCGGAAGGAAGGGGCTAAAAAGCTGaa                               >  2:434045/1‑239 (MQ=255)
                                                                                                                                                                                                              aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGtt                          <  1:448116/239‑1 (MQ=255)
                                                                                                                                                                                                                aCATCGGTGACGAATGCGCCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa                        <  1:726084/239‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGCAGGCTGGCGCTAAACTGCTGAAAAGGTCAGAGt                    >  2:320376/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGGAAACCGTAACCTATGACTTCGAGTGGCTGATGGATGGGGCCAAACTACTAAAATGTTCAGGGtt                   >  2:448116/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTATCTCCGCCGAGCTGATCCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGGAATCAACGCGAAAAACGTAAACTATAACTTCGAGAGTGTATAGTATGGCGCTAAAATGCTAAAAtttttaaag                     >  2:18692/1‑232 (MQ=255)
                                                                                                                                                                                                                            aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGTGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGGTGGAGGGCGCTAAACTGCTGAAATTTTCAAAGTTTggggcc            >  2:15854/1‑235 (MQ=255)
                                                                                                                                                                                                                              tGCGCCCTGTTTGAAACCAACCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGAATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTAGACCGTCTGATGGATGCCGCTAAACTGATGAAATGTTCAGAGTTTGGtacggc          >  1:728157/1‑234 (MQ=255)
                                                                                                                                                                                                                               gggcCCTGTTTGAAACCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg         <  1:143690/237‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGCAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTAAGAGTTTGGTGACGCGatca     >  1:164224/1‑239 (MQ=255)
                                                                                                                                                                                                                                      gTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCg  <  1:15854/239‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
ATCGACGGTGGCCCGTGGCTGAAAGTTAAAAACCCGAACACTGGCAAAGAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCG  >  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 < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: