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,9800AT97.3% 225.7 / ‑3.5 74T530T (ACA→ACTileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (33/38);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
*NC_00091323,9810TC94.6% 224.4 / ‑3.6 74L531L (TTG→CTC) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (33/37);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23753‑24215
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                            
aTGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCgg                                                                                                                                                                                                                       >  1:253088/1‑251 (MQ=255)
  ggTTGCTAACCGTCCGGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTGGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCGCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCGCGACGTATGGTTTGACTCCGGAt                                                                                                                                                                                                                     <  2:277396/251‑1 (MQ=255)
       cTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAACGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTATCGGCATCCAGGCGTGGTGGGATCTCGCTGCGAAAGAGATCCTCGGCGACCAAGCTGACCAGTAGGTGAAAGTGCCGGATA‑‑CTCTCGACGGATGGTTTGACGCCGGATCtccc                                                                                                                                                                                                                >  2:138466/1‑248 (MQ=255)
         aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTAccc                                                                                                                                                                                                               <  1:102422/250‑1 (MQ=255)
         aaCCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTAcc                                                                                                                                                                                                                <  1:67659/249‑1 (MQ=255)
          acattCCTGACTTGTGTAGTTCCCGTCAGCGAATAGGGGACGTATCGATGTCACTGATCGTGAACAAAGAAACGGAAGAGCTGCAACCGCGTCCCCACGAACTGATGGTAGAAGTGGCAAAACTCGTGGCAGTTAACGGTACCCAGGCGTGGTGGGAGATCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTAAAAGTGCCGGATATTCAT‑‑GACGTATGGTTTGACTCCGTATCTACCCat                                                                                                                                                                                                             <  2:2411/247‑2 (MQ=255)
               ccTGGCTGGGGTATCGCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACtctt                                                                                                                                                                                                         <  2:101891/250‑1 (MQ=255)
               ccTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACtctt                                                                                                                                                                                                         <  1:110932/250‑1 (MQ=255)
                cTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGGGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACtcttct                                                                                                                                                                                                       >  2:2811/1‑251 (MQ=255)
                 tGACTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTAAGGTTTGACTCCGGATCTACCCACTCTTCtg                                                                                                                                                                                                      >  2:432619/1‑251 (MQ=255)
                    cTGGTGTATCTCCCGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACCGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGCAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACCCTTCtgttg                                                                                                                                                                                                   >  2:232415/1‑251 (MQ=255)
                               cccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTCCCCACTCTTCTGTTGTTGACGTGCGt                                                                                                                                                                                        <  1:480933/251‑1 (MQ=255)
                                ccGTCAGCGCACCTGGGGCGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTc                                                                                                                                                                                       >  1:471822/1‑251 (MQ=255)
                                                 gcgtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCAc                                                                                                                                                                      >  1:384567/3‑251 (MQ=255)
                                                   gtACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTTACCAGTATGTGAAATTGCCGGATA‑‑CTTTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                                                    <  1:317329/251‑1 (MQ=255)
                                                     aCCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGc                                                                                                                                                                    >  1:115482/1‑249 (MQ=255)
                                                        gATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGTTCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGACTTTGCCGGTCACGCAGCgg                                                                                                                                                               >  2:132675/1‑251 (MQ=255)
                                                         aTGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGa                                                                                                                                                              <  2:476720/251‑1 (MQ=255)
                                                             cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACAt                                                                                                                                                          <  1:31894/250‑1 (MQ=255)
                                                             cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                                         <  1:158470/251‑1 (MQ=255)
                                                             cACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACGCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATg                                                                                                                                                         <  2:47776/251‑1 (MQ=255)
                                                                  ttCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTCCCGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTACCGCAGCggcctagtatcc                                                                                                                                                     >  2:437855/1‑241 (MQ=255)
                                                                        cacaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGAAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGCGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGAATCCGGATCTACCCACTCTTCTGTTGTTGCCGTGGGTCCGGAATTTGCCGGTCACGAAACGGGAATGGATCTGGAAgg                                                                                                                                              >  2:234815/1‑251 (MQ=255)
                                                                         acaAAGACACGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGt                                                                                                                                             <  1:193222/251‑1 (MQ=255)
                                                                               acacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAc                                                                                                                                       <  2:215461/251‑1 (MQ=255)
                                                                               acacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAc                                                                                                                                       <  1:247324/251‑1 (MQ=255)
                                                                                cacGGAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGAcc                                                                                                                                      <  1:197507/251‑1 (MQ=255)
                                                                                   ggAAGAGCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACCTGTATCTGGAAGGTTCTGCCCagc                                                                                                                                   >  2:359122/1‑249 (MQ=255)
                                                                                      agagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCAcc                                                                                                                                <  1:299602/251‑1 (MQ=255)
                                                                                       gagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCGAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                                               <  2:464702/251‑1 (MQ=255)
                                                                                       gagCTGCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcg                                                                                                                               <  1:233390/251‑1 (MQ=255)
                                                                                         gatGCATCCGCGTACCCTCGAACTGTAGGAAGAAGTGGCAAAACGCGTTTAAGTTAACTGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtg                                                                                                                             <  1:287138/249‑3 (MQ=255)
                                                                                            gCATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGAAGTATGGTTTGGCTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTATGGCCAGCACcgtggtt                                                                                                                          >  2:86077/1‑246 (MQ=255)
                                                                                             cATCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTGGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACcgtggttg                                                                                                                         <  2:384567/251‑7 (MQ=255)
                                                                                              aTCCGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGAGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTACGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCCATGTATCTGGAAGGTTCTGACCAGCACcgtggttgg                                                                                                                        >  1:224604/1‑244 (MQ=255)
                                                                                                ccGCGTACCCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGtggttggt                                                                                                                       <  1:205248/250‑1 (MQ=255)
                                                                                                       cccTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCt                                                                                                               <  1:363694/251‑1 (MQ=255)
                                                                                                           cgAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTccc                                                                                                           >  2:333697/2‑251 (MQ=255)
                                                                                                            gAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                                          <  2:150542/251‑1 (MQ=255)
                                                                                                             aaCTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCt                                                                                                          >  2:456419/1‑250 (MQ=255)
                                                                                                              acggATGGAAGAAGGGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGCAGCTGACCAGTATGTGGAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGGTCTACCCACGCTTCTGTTGTTGTCGCGCGTCCGGCATTTGCCGGGAACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCCTCCCTat                                                                                                        >  2:161011/4‑250 (MQ=255)
                                                                                                                 gATGGAAGAAGTGGCAAAACGCGTTTGAGTTAACGGCCCCCAGGCGTGGTTGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCCCCGTTGTTGGTTCATGTCTTCCCTAAt                                                                                                       <  2:297164/249‑1 (MQ=255)
                                                                                                                  aTGGAAGAAGTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGACTTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGAt                                                                                                    >  2:74095/1‑251 (MQ=255)
                                                                                                                     gaagaaGTGGTACAACGCGTTGAAGTCAACGGCATCAAGGCGTGATGAGATCTCGATGCGCAAGAGATCCTCTGCGATGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCGGTTGTTGACTTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCATCACCGTGGTTGGTTCATGTCGTCCCTaagtatctc                                                                                                 <  2:202036/251‑8 (MQ=255)
                                                                                                                       agaagTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCa                                                                                               <  2:77394/251‑1 (MQ=255)
                                                                                                                         aagTGGCAAAACGCGTTGAAGTTAACGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCAcc                                                                                             <  1:198477/251‑1 (MQ=255)
                                                                                                                            tGGCAAAACGCGTTGAAGTTAACGGCATCCATGAGTGGTGGGATCTCGATGCAAAAGAGATCCTCGGCGACGAATCTGACAAGTATGTGAACGTGACGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGCACTGTATCTGGACGGTTCTGCACAGCCCCGTGGTTGGTTCATGCCTTCCCTAATAATCTCCACcgcg                                                                                          >  2:268165/1‑251 (MQ=255)
                                                                                                                                                 aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCagc                                                                                                                                   <  2:78396/187‑3 (MQ=255)
                                                                                                                                                 aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCagc                                                                                                                                   >  1:78396/4‑188 (MQ=255)
                                                                                                                                                 aacggCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGt                                                                       >  1:346015/4‑250 (MQ=255)
                                                                                                                                                     gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGt                                                                  >  1:125609/1‑251 (MQ=255)
                                                                                                                                                      cATCCAGGCGTGGTGGGAGCCCGATGCGAAAGATATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATC‑‑CTCTCGACGGATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTc                                                                 <  2:224604/251‑1 (MQ=255)
                                                                                                                                                         ccAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGTCGAAGCTGACCAGTATGTGAAAGTGCCTGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCGCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAg                                                               <  2:40357/250‑1 (MQ=255)
                                                                                                                                                           aGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGt                                                             >  1:413836/1‑250 (MQ=255)
                                                                                                                                                                    gggATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCAc                                                   <  2:276684/251‑1 (MQ=255)
                                                                                                                                                                     ggATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACg                                                  <  1:311493/251‑1 (MQ=255)
                                                                                                                                                                          gcGATGCGAAAGAGGGCCTCGGCGGTGAAGCTGACCAGTATGTTAAAGTGCCGGATA‑‑CTCTCGTCGTATTGTGTGACTCCGGATCTTCCCACTCTTCGTTTGTGGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGGTCATGGCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCTGTATCGTCAGGTACTGACCCACGGCtt                                              <  2:431823/249‑1 (MQ=255)
                                                                                                                                                                           cGATGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCttt                                             <  2:242327/250‑1 (MQ=255)
                                                                                                                                                                             aTGCGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTAcc                                          <  2:95405/251‑1 (MQ=255)
                                                                                                                                                                                cGAAAGAGATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTg                                       >  1:125602/1‑251 (MQ=255)
                                                                                                                                                                                    agagATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATg                                   >  1:241828/1‑251 (MQ=255)
                                                                                                                                                                                    agagATCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGAACTGACCCACGGCTTTACCGTGGATg                                   >  1:205447/1‑251 (MQ=255)
                                                                                                                                                                                      agaTCCTCGGCGACGCAGCTGACCAGTATGTGCAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCGACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGTCATGTCTCTGGAAGGTTCTGCCCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCCTTACCGTGGATGGt                                 >  2:18457/1‑251 (MQ=255)
                                                                                                                                                                                      agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt                                 >  1:25444/1‑251 (MQ=255)
                                                                                                                                                                                      agaTCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGt                                 >  2:40061/1‑251 (MQ=255)
                                                                                                                                                                                       gaTCCTCGGCGACGATGCTGACCAGTATGTGAAATTGCCGGATA‑‑CTCTCGCCGTAGGGTTTGACTCCGGATCTACCCACTCTGCTGTTGTTGACGTGCGTCCGGAAATTGCCGGTTACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGGCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTATCGTGGATGGTc                                <  2:282799/251‑1 (MQ=255)
                                                                                                                                                                                         tCCTCGGCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTc                                >  2:425456/1‑249 (MQ=255)
                                                                                                                                                                                               gCGACGAAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGc                        <  1:408489/251‑1 (MQ=255)
                                                                                                                                                                                                agacgaAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGAGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATAGTCAGGTACGGACCCACGGCTTTACCGGGGATGGTCacgggcgc                        >  2:5231/2‑243 (MQ=255)
                                                                                                                                                                                                   cgaAGCTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCa                       >  2:306390/1‑248 (MQ=255)
                                                                                                                                                                                                        cTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATg                                                                                                                  <  2:320078/152‑1 (MQ=255)
                                                                                                                                                                                                        cTGACCAGTATGTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGTTCATg                                                                                                                  >  1:320078/1‑152 (MQ=255)
                                                                                                                                                                                                                   gTGAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGGATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGTTGGATCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCgg    >  1:286386/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAAAGTGCCGGATA‑‑CTCTCGACGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAGCACCGTGGGTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAAGGCAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCa  <  1:432619/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                            
ATGGTTGCTAACCGTCCTGACTGGTGTATCTCCCGTCAGCGCACCTGGGGTGTACCGATGTCACTGTTCGTGCACAAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACA‑‑CATTGGATGTATGGTTTGACTCCGGATCTACCCACTCTTCTGTTGTTGACGTGCGTCCGGAATTTGCCGGTCACGCAGCGGACATGTATCTGGAAGGTTCTGACCAACACCGCGGCTGGTTCATGTCTTCCCTAATGATCTCCACCGCGATGAAGGGTAAAGCGCCGTATCGTCAGGTACTGACCCACGGCTTTACCGTGGATGGTCAGGGCCGCAAGATGTCTAAATCCATCGGCA  >  NC_000913/23753‑24215

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