Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,177 C→T E134K (GAG→AAG)  thiF ← sulfur carrier protein ThiS adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,1770CT98.3% 178.6 / NA 60E134K (GAG→AAG) thiFsulfur carrier protein ThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/33);  minor base G (0/1);  total (26/34)
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.95e-01

CAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGT  >  NC_000913/4192931‑4193417
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
cAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGAGCCACGCAGGCGGCATTAATCTTCTgg                                                                                                                                                                                                                                              >  1:387974/1‑251 (MQ=255)
     gtCTCTATACCGCTTAATAACTTAAGGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCAACCCCGCCCGCCGTGCGGCAGCTGCGTTCTGGCGCCTGTTTATCTGGTCACAGGCATCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGCCGGCATTAATCTTCTGTCGAGt                                                                                                                                                                                                                                         <  1:269731/251‑1 (MQ=255)
      tctcTATAACGCTTACTAACTTAATGGCTTCCAGTGCCTGCCAAGTGCCCATAACCCCGAACACCGGGACAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGATGAGGGGCACGCAGGCGGACTTAATCTTCTGGccgggc                                                                                                                                                                                                                                        >  2:234252/1‑246 (MQ=255)
       ctctATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGGGTCg                                                                                                                                                                                                                                       >  2:154615/1‑251 (MQ=255)
         cttgaCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCATCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTc                                                                                                                                                                                                                                        <  2:157468/244‑1 (MQ=255)
          tataCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGcc                                                                                                                                                                                                                                     >  1:68041/1‑250 (MQ=255)
           ataCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCAACCACCGGGCCAACCACCACCCCCGTGCGAAAGTTGCGTTCTAGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGCACGCGCTGGCGGTGATAAGCGGCGTGTTGAGCGCCACGCAGGACGCATTAATCTTCCGGCGAGTCGccct                                                                                                                                                                                                                                   >  1:31175/1‑249 (MQ=255)
                gCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGCCCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTGATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAGCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGt                                                                                                                                                                                                                              <  1:182393/251‑1 (MQ=255)
                 cTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGAGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGGGTCGCCATATTGTc                                                                                                                                                                                                                             >  1:38983/1‑251 (MQ=255)
                     ataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCAAACCCCCGACCACCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTa                                                                                                                                                                                                                         <  2:49927/251‑1 (MQ=255)
                                             gCAAAGTGCCCATACCCCCGCCCACCGGGCCAACCCCGCCCGCCGGGCGGCAGTTGCGTTCTGGCTCCGGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCAGTAATCTTCTGGCGAGTCGCCATATTTTCGGTACAGTCGAGCACCACATCGGccc                                                                                                                                                                                                   <  2:192760/249‑1 (MQ=255)
                                               aaaGGGCCCATCACCCCGCCCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCGCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTGCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCTCCACATCGGCCCGTGc                                                                                                                                                                                               <  2:90439/251‑1 (MQ=255)
                                                            ccccGCCCCCCGGCCCACCCACCCCCGCCGTGCGGCGGTTGCGTTCTGGCTCCGGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTGGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCTCCATATTGTCGGTACAGTCGAGCTCCACATCGGCCCGTCCAACCGCATCtttt                                                                                                                                                                                  <  2:169566/251‑1 (MQ=255)
                                                             cccGACCACCGGGCCAACCACGCCCGCCGAGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGCCCGCCAAATCCGACCGCGCTGGCGGAGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGCGCAACCCATCGGCCCGTGACACCGCACCttttt                                                                                                                                                                                 >  2:262741/1‑250 (MQ=255)
                                                                  cccccGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGcc                                                                                                                                                                            <  2:322501/248‑1 (MQ=255)
                                                                   caccGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCt                                                                                                                                                                           <  1:308184/251‑1 (MQ=255)
                                                                        ggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGAGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGACACAGCATCTTTTAACGCCTCAccc                                                                                                                                                                      >  1:87763/1‑251 (MQ=255)
                                                                           ccacctcgcccgccGTGCGGCAGTTCGGGTCTGGCGCCTGGTTATCTGGCCAAAGGTAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTTCCATCAACTTACCGGCTAATCCGGCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCTCGCAGGCGGCAGTAATCTTCTGGCTAGGCGCCATATTGTCTGTACAGTCGATCACAGCATCGGCCCGTGCAACCGCAGCTTTTAACGCTTCACCCGtt                                                                                                                                                                   <  2:383380/245‑1 (MQ=255)
                                                                            aaCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACGGGCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTa                                                                                                                                                                  <  2:123605/251‑1 (MQ=255)
                                                                                 cgcccgccGTGCGGCATTGGCGTTCTGGCTCCTGTTTAACTGGCCACAGGCCGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCTGTGATAAGCGGCGTGTTGAGTGCCTCGCAGGCGTCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGATCACCACATCGGCCCGTGCAGCCGCATCTTTTAACGCCTCACTCGTGAACCgt                                                                                                                                                             <  2:375561/251‑1 (MQ=255)
                                                                                   cccgccGTGCGGCAGTTGCGTTCTGGGTCGTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCCCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGTCGATTCGCCATATTGTCTGTACAGTCGAGCACCACATCGTCCTGTGCAACCGCTTCTTTTATCTCCTCACCCGTTAACCgttg                                                                                                                                                           <  2:237095/251‑1 (MQ=255)
                                                                                    ccgccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCGCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgttgt                                                                                                                                                          <  1:147250/251‑1 (MQ=255)
                                                                                         gggCGGCATTTGCGTGCTGGCTGCTGGTTATCTGGCCACAGGCAGCGGTAACGCCCCTGCGCCCAGGGCTGCGTCAGTACCATCACTTGACCGCCAGATCCGACCGCGCGGGCGGTGATAGGCGTCGTGTTGAGTGGCGCTCGGGCGGCATTTAGCTTCTGGCGAGGCTCCAGATTGGCGGGTCAGTCGAGCACCACAGCGTCCGTTGCAACCGCATCGTTTAACGCCTCACCCGTTAACCGTTGTTGTaa                                                                                                                                                     <  1:189564/249‑1 (MQ=255)
                                                                                             ggCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAACGTTGTTGTAATGCt                                                                                                                                                 >  1:215800/1‑251 (MQ=255)
                                                                                                    gCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCATGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGt                                                                                                                                           <  2:250884/250‑1 (MQ=255)
                                                                                                    gCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGCCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGt                                                                                                                                           <  2:99990/250‑1 (MQ=255)
                                                                                                     cGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCATTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTg                                                                                                                                         <  1:360943/251‑1 (MQ=255)
                                                                                                        tCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTTTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGaa                                                                                                                                       <  1:204992/250‑1 (MQ=255)
                                                                                                        tCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATAAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACCGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTTGAATGCTGTCAGTTGAAt                                                                                                                                      >  1:227395/1‑251 (MQ=255)
                                                                                                           ggCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGt                                                                                                                                    <  1:154615/250‑1 (MQ=255)
                                                                                                              tCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTGCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTTCCACGCAGGCGGCATTAACCTTCTGGAGCGGCGCCCTATTGTCGGTACAGTCGAGCACCACATCGGCCCGTCTCACCGCATCTTTTAACGCCTCTCCCGTTACCCGTTGTTGCAATGCTGGCAgcggcaaggagtg                                                                                                                                >  2:95172/1‑239 (MQ=255)
                                                                                                                     tATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGCCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGAGTGTTGCGGGCCACGCAGGCGGCATTAATCTTATGGCGAGACGCCATATTGTCGGTACAGTCGAGCACCCCATCGGCCCGTTCCACAGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAAAGCTGTCAGTTGAATGTCGGGATTCAAc                                                                                                                         >  1:383081/1‑251 (MQ=255)
                                                                                                                         tGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGGGAACCACATCGGCCCGTGCAACCGGCTCTTTTAACGCCTCCCCCGTTAAACGTTGTTGTATTGCTGTCAGTTGCATGGCGGGGTTCAACtgtg                                                                                                                     >  2:193748/1‑251 (MQ=255)
                                                                                                                          ggCCACAGGCAGCGGTAACACCCCTGCTCCCAGGTCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTgtgt                                                                                                                    <  1:216081/251‑1 (MQ=255)
                                                                                                                                    aGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATTTCGGGGTTCCACtggttcagccccttt                                                                                                          >  2:137260/1‑238 (MQ=255)
                                                                                                                                                      cccAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTc                                                                                        <  1:262741/251‑1 (MQ=255)
                                                                                                                                                       ccAGGGCGGCGTCAGTACCACCAAATGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGCGTTGAGTGCCACGCCGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCCCATCCGTCCCACCCACCGCCTCTTTTAACCCCTCCCCAGTTCACCGTTGGTGTAATGCTGACAGTTGAATGTCGGGCTTCAACTGCTTCAGTCGCTGGTGGCTGGCCttcggttttg                                                                                       >  2:26057/1‑242 (MQ=255)
                                                                                                                                                          gggCGGCGTCAGTACCAACAACTGACCGCCAAATCCGGACGCGCTGGCGGTGATAAGCGGCGGGTTGAGTGCCCCGCAGGGGGCATTAAACTTCTGGGGAGTCGCCATATTGTCGGGTCAGTCGAACACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                    <  1:137260/251‑1 (MQ=255)
                                                                                                                                                                     gTACCATCAACTGCCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCCTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATACTTTCAGTTGAATGTCGGGATTCAACTGTGTCATTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGCTATCt                                                                         >  1:73699/1‑251 (MQ=255)
                                                                                                                                                                         cATCAACTGACCGCCAAATCCGACCGCGCGGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAg                                                                     <  1:291126/251‑1 (MQ=255)
                                                                                                                                                                           acAACTGACCGCCAAACCCGATCGCGCTGTCGGGTATAAGCGGCGTTTTGATTGCCCCGCATGCTGCACTAATCTTCGGGAGTGGCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCGCGTGCAACCGCATCTTTTAACGCCTCTCCCGTTATCCGTTGTTGTAATTCTTTCAGTTTAATGTCGGGATTCAACGGTGTCAGTCGTGGTTGGCTGACCTGCGAGTTCGGGCGATCTATATCTTCattg                                                                   <  2:383081/250‑4 (MQ=255)
                                                                                                                                                                                tGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTaa                                                               >  1:394204/1‑250 (MQ=255)
                                                                                                                                                                                   ccGCTAAATCCGACCGCGCGGGCGTTGATAAGCGGCGTGTTGAGTTCCACGCAGGCGGCATTAATCTTCTGGCGGGGCGCCATATTGTCGGTACAGTCGAGCGCCACATCGGCCCGGGCAACCGCATCTTTTAACGCCTCGCCCGTTAACCGTTGTGGTAATGCTGGCAGTTTAATGTCGGGATGCAACTGTGTCAGTGGTTGTGGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAagag                                                           <  1:254359/251‑1 (MQ=255)
                                                                                                                                                                                            ccgtccGCGCTGGCGGTGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCCCCACATCGGCCCGTGCAACCGCATCTTTTATCGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCg                                                  <  2:73699/247‑1 (MQ=255)
                                                                                                                                                                                            ccGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCg                                                  <  1:273023/251‑1 (MQ=255)
                                                                                                                                                                                             cGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCAAGCAGGCGGCATTAATCTTCTGGCGAGTCGCCCTATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAAGCCTCACCCGTTAAACGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGAGTTCGGGCGATCGATATCTTCAGGGGTAAAGGGGATTTGttgt                                                 >  2:228080/1‑248 (MQ=255)
                                                                                                                                                                                                ccGCGCTGGCTTTGATAAGCGTCGTTTTGCGTGCCGCGCAGGCGGCATTAATCGTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTATCGCCTCACCCGTTAACCTTGGTTGTAATGCTTTCATTTGTATGTCGGGATTCATCTTTGTCAGTCGCTGTTGGCTGTCCTGCGAGTTCGGGCGATCGATATCTTCAGTGGTCAAGAGGATTTGTCGTTGc                                              <  2:179240/251‑1 (MQ=255)
                                                                                                                                                                                                  gcgcTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAg                                            <  1:109733/251‑1 (MQ=255)
                                                                                                                                                                                                    gcTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGAGTCGCCATATTGTCGGTACCGTCGAGCACCACATCGGCCCGTGCCACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTTGCTGACCTGCGATTTCGGGCGATCGCTATCTTCAGTGGTAAAGGGGGTTTTTCGTTGCCGAt                                          >  1:394262/1‑251 (MQ=255)
                                                                                                                                                                                                         cGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGGAAAGAGGATTTGTCGTTGCAGATTGCtt                                     >  1:249029/1‑251 (MQ=255)
                                                                                                                                                                                                          ggTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTTGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGcgta                                    <  2:322096/251‑4 (MQ=255)
                                                                                                                                                                                                                     gcgTGTTGAGTGCCACGCAGGCGGCATTACTCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGCTTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAACTGCACAtcg                         >  2:259464/1‑251 (MQ=255)
                                                                                                                                                                                                                       gtgtTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtc                       >  1:367342/1‑251 (MQ=255)
                                                                                                                                                                                                                            gAGTGCCACGCAGGCGGCATTCATCTTCTAGCGAGTCGCCATATTGTCGGAACAGTCGTGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTAAGCTGAATGTCGGGATTCAACTGTGTCCGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATCTCTTAAGTGGTAAAGAGGCTTTGTCGTTGCAGAGTGCTTAAATGCACATCGTCGTCAt                  >  2:181369/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                 >  1:21968/1‑251 (MQ=255)
                                                                                                                                                                                                                              gTGCCACGCAGGCTGCATTAATCTTCTGGCGAGTCGCCATATGGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCt                <  2:97044/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACCGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATGTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGcc             <  1:34844/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGcc             >  1:15182/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGACGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTAcc       >  1:273589/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gcATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACTCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGt  <  2:299680/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
CAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGT  >  NC_000913/4192931‑4193417

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