Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,848,239 G→A F146F (TTC→TTT ydjA ← putative oxidoreductase YdjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,848,2390GA98.7% 254.4 / NA 74F146F (TTC→TTTydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/34);  minor base C (1/0);  total (40/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 = 1.00e+00

CGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCC  >  NC_000913/1848002‑1848483
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
cGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGcc                                                                                                                                                                                                                                         >  1:547689/1‑250 (MQ=255)
 gCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCAc                                                                                                                                                                                                                                        >  1:14741/1‑251 (MQ=255)
  cTCATGGACATCTCCTGTCAATGCAATCCGTGTATGGTAACCCTCATTCCGTGACGTGATAAGTGAGAGTGTCTGAATTCCTGCGACTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGATACAAACGGCGTCGGGTCCGGGACGTTAATCGGCGAAGATGATTTCAGCTGCGGCGTACCGCGGTAGGGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCAcc                                                                                                                                                                                                                                       >  2:533956/1‑250 (MQ=255)
           aTCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACtt                                                                                                                                                                                                                              >  2:342149/1‑251 (MQ=255)
                 tGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCCGtt                                                                                                                                                                                                                        >  1:240637/1‑251 (MQ=255)
                  gTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGGCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGt                                                                                                                                                                                                                         <  1:48893/249‑1 (MQ=255)
                      aTGCAATCCGGTTATGGTAAACCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTAAGCTGCGGCGTACCGAGGTAGAGAAAATCGACAATTTTATCCTGCTCACGGCAACCAACTGCTTCACGCACTACCGGACTTTCAGTTAATgc                                                                                                                                                                                                                   >  2:492243/1‑251 (MQ=255)
                      aTGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATgc                                                                                                                                                                                                                   <  1:247220/251‑1 (MQ=255)
                             ccGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAAGCCAGACAGTTGCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAAGCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACGACCGGACTTTCAGTTAATGCGCCACt                                                                                                                                                                                                             <  2:482958/250‑1 (MQ=255)
                              cGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTgc                                                                                                                                                                                                           <  1:599491/251‑1 (MQ=255)
                              cGGGTATGGTAAACCTCATTCCGTGAAGTGATAAGTGAGCGTGTCTGCATTCCTGCGCCTTTGCTCACAATACAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTCCCGGCCTTTCAGTTAATGCGCCACTgc                                                                                                                                                                                                           >  2:68344/1‑251 (MQ=255)
                                 gTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGcc                                                                                                                                                                                                        <  1:101333/251‑1 (MQ=255)
                                   aTGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAg                                                                                                                                                                                                      <  2:148731/251‑1 (MQ=255)
                                                   cGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAAGCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACGTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCt                                                                                                                                                                                      <  2:240637/251‑1 (MQ=255)
                                                    gTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTg                                                                                                                                                                                     >  1:163581/1‑251 (MQ=255)
                                                       aaGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGc                                                                                                                                                                                  >  1:691544/1‑251 (MQ=255)
                                                       aaGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGc                                                                                                                                                                                  >  1:574105/1‑251 (MQ=255)
                                                           gATAAGTGAGAGTGTCTGAATTCCTGCGCCTTGGCTCCCAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCATTGCGCCAGATGCCGCCAATCCCCTGGGCAAct                                                                                                                                                                              <  2:525417/251‑1 (MQ=255)
                                                              aaGTGAGAGTGTCTGAATGCCTGCGCCTGTGCTTACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAATGGCGTCGGGTCCGGGACGTTAATCGGCGTAGATGCTTTCAGCTGCGGCGGACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAActgct                                                                                                                                                                           <  2:629753/251‑1 (MQ=255)
                                                                gTGAGAGTGTCTTAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAATCGGCGTCGGGTCCGGGACGTTAATCGACGTATATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGGCTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGACAACggctgc                                                                                                                                                                         <  2:677318/251‑1 (MQ=255)
                                                                   agagTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACACTTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCAt                                                                                                                                                                      >  1:480097/1‑251 (MQ=255)
                                                                            gAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGGTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGTCATTTGCATc                                                                                                                                                               <  2:673417/249‑1 (MQ=255)
                                                                             aaTTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTTCTGCCATTTGCATCGcc                                                                                                                                                            <  1:192896/251‑1 (MQ=255)
                                                                                    gcgcCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGAc                                                                                                                                                       <  2:691544/249‑1 (MQ=255)
                                                                                    gcgcCGTTGCTCCCAATCCAGACAGTGTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGAcc                                                                                                                                                      <  2:572023/250‑1 (MQ=255)
                                                                                     cgcCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACcgc                                                                                                                                                    >  1:47723/1‑251 (MQ=255)
                                                                                         tttGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCAt                                                                                                                                                >  2:58471/1‑251 (MQ=255)
                                                                                          ttGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCAt                                                                                                                                                <  1:263651/250‑1 (MQ=255)
                                                                                           tGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATc                                                                                                                                               <  2:171008/250‑1 (MQ=255)
                                                                                                  aaTCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGAc                                                                                                                                       >  1:599869/1‑251 (MQ=255)
                                                                                                  aaTCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGAc                                                                                                                                       >  2:491288/1‑251 (MQ=255)
                                                                                                     ccAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACAtt                                                                                                                                    >  2:126752/1‑251 (MQ=255)
                                                                                                        gACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACAtttac                                                                                                                                 <  1:203377/251‑3 (MQ=255)
                                                                                                         aCAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCt                                                                                                                                >  1:102087/1‑251 (MQ=255)
                                                                                                         aCAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCt                                                                                                                                <  2:574105/251‑1 (MQ=255)
                                                                                                                cgcgACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTccc                                                                                                                          <  2:269274/250‑1 (MQ=255)
                                                                                                                 gcgACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTccc                                                                                                                          >  1:659827/1‑249 (MQ=255)
                                                                                                                    aCAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAgcgc                                                                                                                     >  1:317509/1‑251 (MQ=255)
                                                                                                                       aTTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCggg                                                                                                                  <  2:689490/251‑1 (MQ=255)
                                                                                                                        ttATCAGAAATAAGTTACAACCGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGATTTCAGCTGCAGCGTACCGAGGTAAAGAACACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGCTGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGGCCGCGCATCCGGCCACCATTTCCTGTGCCCAGcgcgtg                                                                                                                  >  2:601338/1‑248 (MQ=255)
                                                                                                                            cAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                             >  1:134998/1‑251 (MQ=255)
                                                                                                                             aGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTa                                                                                                            <  1:491288/251‑1 (MQ=255)
                                                                                                                                   aaGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATtc                                                                                                      <  1:161651/251‑1 (MQ=255)
                                                                                                                                       tacaaaacgagacgggcccgggccGTTACTCGACGTAGATAATTTCAGCTGCGACGCACCGAGGTAGAGAAAACCGAAAATTTTACCATACTCACAGCAACCCACTGCATCACGCACTACCGGACTTTCAGATAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGTCAACTGCTGCCATTTGAATCGCCATCACAGCGCTTCGGGCAGACATTTCCTGTTCCCAGCGCGTGACTTTATGATTCCCtt                                                                                                  >  2:441781/13‑251 (MQ=255)
                                                                                                                                         gaaaCGGCGTCGGGTCCGAGACGTTAATCGAAGTAGATGCTTTAAGCTGCGGCGTACCGAGGGAGAGCATACCGCCACTTTTATCCTGCTCCCGGCAACCAAATGCTTCACGCACTCCCGGAATTTTAGTTAATGCGCCACTGCGCCAGCTGCCGCCCAACAACTGGGCAACTGCTGCCATTTGCCTCGACCTGACCGAGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTCTGCTTCTCTTCg                                                                                                >  2:653994/2‑251 (MQ=255)
                                                                                                                                         cAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTACTGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCg                                                                                                >  2:491307/1‑251 (MQ=255)
                                                                                                                                               gCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTc                                                                                          >  1:370641/1‑251 (MQ=255)
                                                                                                                                                             aCGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGCTTCTCTTCGCATTTCGCCACCACCGtgat                                                                            >  1:760/1‑251 (MQ=255)
                                                                                                                                                                 tAATCGACGTAGATGCTTTCAGCCGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCACCCAAATGATTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTACCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGtgatgatg                                                                        >  2:648766/1‑251 (MQ=255)
                                                                                                                                                                       aCGTAGATGCTTTCAGCTGCGGCGTTCCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACTACCGTGATGATGAGCGGt                                                                  <  1:285472/251‑1 (MQ=255)
                                                                                                                                                                            gATGCTTTAAGCGGCGGCGTACAGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTGTATGATTCTCTTCGCATTTCGCCACCACCGTGAGGAAGAGCGGTGcgcg                                                             >  2:319674/1‑251 (MQ=255)
                                                                                                                                                                                  ttCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAAc                                                         <  1:16166/249‑1 (MQ=255)
                                                                                                                                                                                   tCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACgg                                                       <  2:760/250‑1 (MQ=255)
                                                                                                                                                                                           ggcgTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACg                                              >  1:195302/1‑251 (MQ=255)
                                                                                                                                                                                                 ccGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACTGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCttt                                        <  1:432497/251‑1 (MQ=255)
                                                                                                                                                                                                 ccGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCttt                                        <  1:252515/251‑1 (MQ=255)
                                                                                                                                                                                                  cGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTg                                       >  2:472083/1‑251 (MQ=255)
                                                                                                                                                                                                    aGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTc                                     <  1:698948/251‑1 (MQ=255)
                                                                                                                                                                                                      gTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATg                                                                                                                                                         >  1:36967/1‑133 (MQ=255)
                                                                                                                                                                                                      gTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATg                                                                                                                                                         <  2:36967/133‑1 (MQ=255)
                                                                                                                                                                                                             aaaCCGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATTACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCttt                            >  1:223009/1‑251 (MQ=255)
                                                                                                                                                                                                               cccGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTTCTTCCATTTGCATCGCCATGACCGCTCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATGCTCTGCGCATTTCGCCACCACCGTGATGATGAGCGTTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGt                          <  1:424986/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtc                         >  1:13744/1‑250 (MQ=255)
                                                                                                                                                                                                                 cGACAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtc                         >  2:195087/1‑250 (MQ=255)
                                                                                                                                                                                                                   aCAATTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCAATGCCCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCCGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtcatc                      >  1:471177/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAATTTTATCCTGCTCAAGGCCACCAAATGCTTCCCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtcatc                      >  1:570390/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaTTTTATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCAc                    <  1:195087/251‑1 (MQ=255)
                                                                                                                                                                                                                       ttttATCCTGCTCACGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACt                   >  2:217607/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGGGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCg     <  1:68344/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATAGCCTTGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGGGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGc    >  2:244275/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcc   <  1:58471/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGGCATCACTACCGGCAGCAATCGcc   >  1:544878/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccc  <  2:195302/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gCAACCAAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTCGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGccc  >  2:118590/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
CGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCC  >  NC_000913/1848002‑1848483

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