Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 23,980 2 bp→TC coding (1590‑1591/2817 nt) ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091323,9800AT94.4% 211.0 / 3.9 73T530T (ACA→ACTileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (27/40);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01
*NC_00091323,9810TC93.2% 195.2 / 3.7 73L531L (TTG→CTC) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (27/41);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23748‑24215
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                           
gggcGATGGTTGCGAACCGGCCTGACTGGGGTAACTCCAGAAATCGCACCAGGGGCGGACCGATGTCACTGTACGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATTCAGTCGTGGTGGGATCTCGATACGAAAGAGAGCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTggc                                                                                                                                                                                                                           <  2:249766/248‑3 (MQ=255)
     aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCgg                                                                                                                                                                                                                      <  1:470153/251‑1 (MQ=255)
     aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCgg                                                                                                                                                                                                                      <  2:435776/251‑1 (MQ=255)
             tAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGGATGTGAAAGTGCCGGACACTCTCGCCGAATGGTTTGCCTCCGGATCTcccc                                                                                                                                                                                                              >  2:154678/1‑251 (MQ=255)
              aaCCGGCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTGCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACGGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAGCGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACGCTCGACGTATGGTTTGACTCCGGATCTACCCa                                                                                                                                                                                                             <  2:411794/251‑1 (MQ=255)
               aCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCAc                                                                                                                                                                                                            >  1:238582/1‑251 (MQ=255)
                         cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGCTGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGACAGTGCCGGATACTCTCGCCGTATGGTTTGACTCCGGATCTACCCCCTCTTCtgttg                                                                                                                                                                                                  >  2:369942/1‑251 (MQ=255)
                         cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGCCCAGTATGTGAAAGTGCCGGATCCTCTCGACGTATGGTTTGACTACGGATCTACCCACTCTTCtgttg                                                                                                                                                                                                  >  2:420372/1‑251 (MQ=255)
                           ggTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCGTCtgttg                                                                                                                                                                                                  <  1:226231/249‑1 (MQ=255)
                               tATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACg                                                                                                                                                                                            >  1:306020/1‑251 (MQ=255)
                                  ctcCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGc                                                                                                                                                                                         >  1:46217/1‑251 (MQ=255)
                                    cccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGt                                                                                                                                                                                       >  1:174440/1‑251 (MQ=255)
                                       gTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCg                                                                                                                                                                                    <  2:83428/251‑1 (MQ=255)
                                          aGCGCAACTGGCGCGTACCAGTGTCCCTGTTCGTGCACAAAGGCACGGCAGAGCTGGATTCGAGTACCCTCGAACTGAGGGAAGAAGGTGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCGGTTGTTGACGTGCGTCCGGaa                                                                                                                                                                                 <  2:118366/251‑1 (MQ=255)
                                           gcgcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTTCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAt                                                                                                                                                                                <  1:355593/251‑1 (MQ=255)
                                             gcACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAttt                                                                                                                                                                              <  1:231656/251‑1 (MQ=255)
                                              ccccTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACAAGGAAGAGCTGCAGCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGAGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTg                                                                                                                                                                             <  2:225673/249‑1 (MQ=255)
                                              cACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGACGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTACGGAATTTg                                                                                                                                                                             >  2:415489/1‑251 (MQ=255)
                                                ccTGGGGCGTACCGATGGCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGcc                                                                                                                                                                           <  1:360500/251‑1 (MQ=255)
                                                       cgtacctaggtcactgtgccggCACAAAGACACGGAAGAGCTGCATCCTCGGACCCGCGAACTGATGGAAGAAGTGGCAAAGCGCGTTGAAGTTATCGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACTAATCTGACCAGTAGGTGAAAGTGCCGGATACTCTCGACGTATGGTTGGATTCCGGATCTTCCCACTTTTCTTTTGTTGACGTGCGTCCGGAATTTGCCGGTCACg                                                                                                                                                                    <  2:463813/230‑1 (MQ=255)
                                                        gtACCGATGTCACTGCTCGGGCTCAAGGACTCTGTAGTGTGGCATCCGCGTACACCCGAACTGATGGGAGAAGTGGCAAAGCGCGTTGAAGTTAACTGCATCCAGGCGGGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGAGCCTGATACTCTCTTCGTAGGGTGCGACTCCGGAGCTACCCACTCGTCTGTTGGTGACGTGCGTCCGGAATTTGCCggctacgc                                                                                                                                                                   <  2:458354/251‑7 (MQ=255)
                                                            cGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCg                                                                                                                                                               >  1:119022/1‑251 (MQ=255)
                                                                gTCACTGTTAGAGCACAAAGAAACGGAAGAGATGCCTCCGCGTACCCTCGACCTGCTGGAAGAAGTGGCAAACCGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGTCGTCTGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGCCGTGAGTCCGGAATTTGCCGGTCACGCAGCGGAcc                                                                                                                                                           >  2:500538/1‑250 (MQ=255)
                                                                 tCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAt                                                                                                                                                          >  1:213132/1‑251 (MQ=255)
                                                                  cACTGTTCGTGCACAAATGCACCGAAGATCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAt                                                                                                                                                          <  1:30699/250‑1 (MQ=255)
                                                                  cACTGTTCGTGCACAAAGACTCGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGAGCTCGATGCGAAAGAGATCCTCGGCGACGAAGTTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                                         <  2:177918/251‑1 (MQ=255)
                                                                      ggtCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCGCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCGTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTa                                                                                                                                                       <  2:137850/247‑1 (MQ=255)
                                                                        tCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTg                                                                                                                                                   >  1:289218/1‑251 (MQ=255)
                                                                               caAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGCCCGGAATTTGCCGGTAACGCAGCGGACATGTATCTGGAAGGtt                                                                                                                                            <  1:441694/251‑1 (MQ=255)
                                                                                aaaGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCCTCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGCAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAACGGCCATGTATCTGGAAGGTTc                                                                                                                                           >  2:131690/1‑251 (MQ=255)
                                                                                           agagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGTAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGTCGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGCATTTGCCGGTCACGCAGCGGCCATGTATCTGGTAGGTTCTGAcccgcacc                                                                                                                                >  2:150548/1‑245 (MQ=255)
                                                                                            gagCTGCATCCGCGTACCCTCGAACTGATGTAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCTGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                                               <  1:23998/251‑1 (MQ=255)
                                                                                                 gCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggtt                                                                                                                          <  1:262738/251‑6 (MQ=255)
                                                                                                   aTCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggttgg                                                                                                                        <  1:385456/251‑8 (MQ=255)
                                                                                                          tACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCCGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATgg                                                                                                                 >  2:11691/1‑250 (MQ=255)
                                                                                                           aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGt                                                                                                                 <  1:371679/250‑1 (MQ=255)
                                                                                                           aCCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                >  1:7771/1‑251 (MQ=255)
                                                                                                            cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTc                                                                                                                <  1:502111/250‑1 (MQ=255)
                                                                                                            cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCt                                                                                                               >  1:366720/1‑251 (MQ=255)
                                                                                                              cgcgAACTGATGGAAGAAGGGGCCAAACGCGTTGCAGGTTACGGCATCCAGGCGTGGTGGGATCTATATGCGAAAGAGATCCTCGGCGAAGAAGCTGACAAGTATGTGAAATTGCCGGATACTCTCGACGTATGGTTTGCCTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCTGCACCGTGGTTGGTTCATGTCtgc                                                                                                             <  2:373539/248‑3 (MQ=255)
                                                                                                                cgAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTccc                                                                                                           >  1:8109/2‑251 (MQ=255)
                                                                                                                 gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTAccc                                                                                                                                                                                                              <  1:470990/151‑1 (MQ=255)
                                                                                                                 gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTAcac                                                                                                                                                                                                              >  2:470990/1‑149 (MQ=255)
                                                                                                                            agaagTGGCAAAACGCGTTGCAGTTAACGGCATCCAGGCGTAGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTCGACGTATGGATTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACAGTGGTTGGTTGATGTATTCCCTAATGATCTccc                                                                                               >  2:333189/1‑250 (MQ=255)
                                                                                                                             gaagTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCAc                                                                                              <  1:73381/251‑1 (MQ=255)
                                                                                                                                      aaaCGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGaa                                                                                     <  1:170297/251‑1 (MQ=255)
                                                                                                                                          cgcggtgaaattaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAATATATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCGTCTTTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTGCATGTCTTCCCTAATGATCTCCACCGCGATGaaagg                                                                                  <  2:420144/236‑4 (MQ=255)
                                                                                                                                             ttGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGt                                                                                                                                             <  2:144464/188‑1 (MQ=255)
                                                                                                                                             ttGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGt                                                                                                                                             >  1:144464/1‑188 (MQ=255)
                                                                                                                                             gttgaagttaacggCATCCAGGCGGGGTGGGATCTCGATGAGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGGATCTGGAAGGTTCTGACCCGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGCAGAAgggca                                                                                >  2:405826/13‑248 (MQ=255)
                                                                                                                                                   gttaacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGcc                                                                         <  1:375991/245‑1 (MQ=255)
                                                                                                                                                      aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGt                                                                       >  1:326704/4‑250 (MQ=255)
                                                                                                                                                         ggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCg                                                                   <  1:426267/251‑1 (MQ=255)
                                                                                                                                                          gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGt                                                                  <  1:286426/251‑1 (MQ=255)
                                                                                                                                                                   cGTGGTGGGATCGCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTGCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTg                                                         <  1:293340/251‑1 (MQ=255)
                                                                                                                                                                         gggagcccGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCAc                                                   <  2:158392/244‑1 (MQ=255)
                                                                                                                                                                           gATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCAACCGCGATGAAAGGAAAAGCGCCGTATCGTCAGGTACTTACCAACgg                                                 >  2:368565/1‑251 (MQ=255)
                                                                                                                                                                                 gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTc                                                                                                                                           <  1:201754/155‑1 (MQ=255)
                                                                                                                                                                                 gATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCtttcc                                           >  2:51790/1‑249 (MQ=255)
                                                                                                                                                                                 gATGCGAAAGAGATCCTCGGCGAAGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTc                                                                                                                                           >  2:201754/1‑155 (MQ=255)
                                                                                                                                                                                       aaaGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGa                                     <  1:370176/251‑1 (MQ=255)
                                                                                                                                                                                        aaGAGATCCTCGGCGACGAAGCTGCCGAGTATGTGAAAGTGCCGGATTCTCTCGACGTATGGTTTGACTCCGGACCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATTTATCTGGAAGGTTCTGACCAGCGCCGTGGTTGGTTCATGGCTTTCTTGATGATCGCCTCCGTGATGAAAGGCAAGGCGCCGTAGCGTCAGGATCTGATCCAGGGCTCTACCGTGGAt                                    <  2:211957/251‑1 (MQ=255)
                                                                                                                                                                                         agagATCCTCGGCGACGAAGCTGACCAGTAAGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCTTCAGGTACTGATCCACGGCTTTACCGTGGATg                                   <  2:289218/251‑1 (MQ=255)
                                                                                                                                                                                           agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt                                 <  1:481228/251‑1 (MQ=255)
                                                                                                                                                                                                      gacgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGAGGATGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGAACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                       >  2:384546/1‑250 (MQ=255)
                                                                                                                                                                                                       acgaAGCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                       >  1:259883/1‑249 (MQ=255)
                                                                                                                                                                                                          aaGCTGACCAGTATGTGAAAGTGCCGGATACTCGCGACGGATGGTTTGACTCCGGATCTACCGACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCTATGAAAGGCAAAGCGCCGTATCGGCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATg                  <  2:56897/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTGACCAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTc                <  1:427841/251‑1 (MQ=255)
                                                                                                                                                                                                               gaccAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAAGTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTTACCCACGGTTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTaa             <  1:191976/248‑1 (MQ=255)
                                                                                                                                                                                                                  cAGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACGCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTTTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCTGCAAGATGTCTAAATc          <  1:409289/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGTATGTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATcc         <  2:358757/251‑1 (MQ=255)
                                                                                                                                                                                                                        gTGAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAGGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAGCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGTCCCACGGCTTTACCGTGGATGGTTAGGGCCGCAAGATGTCTAAATCCATCgg    <  1:444009/251‑1 (MQ=255)
                                                                                                                                                                                                                          gAAAGTGCCGGATACTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGcc  >  2:253612/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                           
AGTCGATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23748‑24215

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