Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,533 C→T Y213Y (TAC→TAT serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5330CT96.4% 162.1 / ‑2.3 57Y213Y (TAC→TATserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/25);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

CCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625286‑4625776
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
ccAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTg                                                                                                                                                                                                                                                  <  1:185881/251‑1 (MQ=255)
 aaaaCTGGCCGGAACGGGGGAGAGGGGGGGGGAAGTACCCGAACGGGGGAGGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTgc                                                                                                                                                                                                                                                 <  2:60887/250‑1 (MQ=255)
               cGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTgcgc                                                                                                                                                                                                                                   >  1:245572/1‑251 (MQ=255)
                gggCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGcc                                                                                                                                                                                                                                  <  1:449862/251‑1 (MQ=255)
                    gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGAc                                                                                                                                                                                                                              >  2:189852/1‑251 (MQ=255)
                                gaagaaccctaaCGGGCGATGCGCGGCGACTTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATCCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCGGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGc                                                                                                                                                                                                                  <  1:485913/241‑1 (MQ=255)
                                             gggggattcgcgggcaacccGATTTTACCGCAAGCCTGCGCAGCGGTGTGGCGACGCTGAAAGGCGCTGCCGCCAATATTCTGCAACAGGTGCGAGAAAATCGGCCGCTGATGCCATGCTTAACGCATCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGgag                                                                                                                                                                                                     <  2:312086/232‑1 (MQ=255)
                                             gggCGAGGCGCGGAAACCTCGTTTTTCCCGCCAGCCGGCGCAGCCGCGTGGCGACGCGGAAATGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTACCGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCTGTGGCCTCCGGCGGCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGgag                                                                                                                                                                                                     <  2:33711/251‑1 (MQ=255)
                                                cGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                                                                                                                                                                                       >  1:92325/1‑246 (MQ=255)
                                                     cgcgGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTATGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATgg                                                                                                                                                                                              >  1:22900/1‑250 (MQ=255)
                                                        ggcgAACTCGTTTTTCCCGCCAGCCTGCGCAGCGGCTGGGCGACGCTGAAAGGCGTTGAGCACAATATTCTGAAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGTTTAACGCAACTGGTGCCTAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGTTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACgg                                                                                                                                                                                          <  2:436564/251‑1 (MQ=255)
                                                        ggcgAACTCGATTTTCCAGCCAGCCGGCGCAGCCGCGTGGCGACGCGGAAAGGCGTTGACGCCAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCGGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCCGAAGATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCgg                                                                                                                                                                                          <  2:153471/251‑1 (MQ=255)
                                                           gAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACggaa                                                                                                                                                                                        >  1:372213/1‑248 (MQ=255)
                                                                cGATTTTACCGCCAGCCTGCGCAGCCGCGTGGGGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAtt                                                                                                                                                                                    <  1:321889/249‑1 (MQ=255)
                                                                 gtttttACCGCCAGCCTGCGAAGCCGCGTGGCGACGCTGAAAGTGGCGGACGCCAATTTTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAc                                                                                                                                                                                 <  2:157111/249‑1 (MQ=255)
                                                                  aTTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                                >  2:346183/1‑251 (MQ=255)
                                                                      aaCCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                             <  1:410689/249‑1 (MQ=255)
                                                                        ccGCCAGTCTGCGCAGCCTCGGGGCGACGCGGAAAGGCGCTGACGCCAATATCCTGCAACAGGTGCGTGAAACTCTGCCGATTATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCTCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAt                                                                                                                                                                          <  1:312712/251‑1 (MQ=255)
                                                                                    gcAGCCGCGTGGAGACGATGAAAGGAGCTGACGCCCATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGATGGAAAGTGGCGATTGCCTCCGGCGGATTTACTTTCTTTGCTGAATATCTGCGCGAACGGTTGAGCCTGAACGCAGTGGTAGACAATGAACTGGGGATCATGGACGGTAAATTTCCCGGGAATGTGAGCGGGGAc                                                                                                                                                              >  2:276662/1‑251 (MQ=255)
                                                                                                 cGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAcgcg                                                                                                                                                  >  2:257626/1‑251 (MQ=255)
                                                                                                   aCGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACgcgc                                                                                                                                                 >  1:407433/1‑250 (MQ=255)
                                                                                                     gCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATAAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGt                                                                                                                                              >  2:53541/1‑251 (MQ=255)
                                                                                                               gcTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGa                                                                                                                                     >  1:305748/1‑250 (MQ=255)
                                                                                                                  gACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAGCGTAGACGCGCAGTACAAAGCGaaa                                                                                                                                   <  1:412162/249‑1 (MQ=255)
                                                                                                                   ttgaccatttttCTGCAACGGGTGCGTTAAATACGGCCGCGGATGCCAGGTTTAACGCAACTGTGCCTCAAGCTGGAAACGCTGTGTTGGAAAGTGGCTATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCGGCGCGACAAGTTGCCCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAGGCGaaaa                                                                                                                                  <  2:347024/240‑1 (MQ=255)
                                                                                                                      ccAATATTCTGCAGCAGGTGCGTGAAAATCGGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGTCTGGAAAGTGGCGCTTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGATGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGGACAAAGCGAAAACTCTg                                                                                                                             >  2:263859/1‑251 (MQ=255)
                                                                                                                         aaTATTCTGCACCAGGCGCGTGAAAATCTGCCGCTAATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTAATTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGAGCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAACACTCTGAc                                                                                                                           >  2:248169/1‑251 (MQ=255)
                                                                                                                         aaTATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAAGCGCAGTACAAAGCGAAAACTCTGAc                                                                                                                           >  1:239874/1‑251 (MQ=255)
                                                                                                                         aaTATTCCTCAACAGGTTCGGGAAAATCTGCCGCTGAGGCCAGGCTTAACGCAGCTGGGGCTCAAGCGTGAAGCGCTGGGCTGGAAAGTGGCGATTGCCGCCGGCGGTTTTACTTTGTTGGCTGAATGTCTGCGCGGTAGGTTGCGCCTGGCCGCCGTGGTAGCCAATGAACTGGATATCATGGACGGTAAATTTGCCTGCAATGTGATCGGCGACATCTTAGACGCGCAGTACAAAGCGAAAACTCTGAc                                                                                                                           <  1:210669/251‑1 (MQ=255)
                                                                                                                            atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACaaa                                                                                                                                         >  1:104976/1‑234 (MQ=255)
                                                                                                                            atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActcg                                                                                                                        <  1:272965/251‑1 (MQ=255)
                                                                                                                            atTCGGCAACAGGTGCGGGAAAATCGGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACaaa                                                                                                                                         <  2:104976/234‑1 (MQ=255)
                                                                                                                                      aggtcggtgcaaaTCTGCCGCTGATGCCAGGCTTAACGCACCTTGTGCTCAAGCTGGAAACGCTTGGCTTGAAAGTGGCTATTGCCTCCGGCGGCTTTTCTTTCTTTGCTGAATATCTGCGCTACAAGTTGCGCCTGTCCGCCGTGGTAGCCAATGAACTGTAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAg                                                                                                              <  1:266960/241‑1 (MQ=255)
                                                                                                                                      aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                                >  2:478142/1‑249 (MQ=255)
                                                                                                                                            gTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                        >  1:72649/1‑251 (MQ=255)
                                                                                                                                                               gCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGa                                                                                     <  1:132375/251‑1 (MQ=255)
                                                                                                                                                                 aaGGCTTAACGCAACTGGGGTTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGAcc                                                                                   <  1:12268/250‑1 (MQ=255)
                                                                                                                                                                   ggCTTAACGCAACTGGTGCTCAAGCTGGAACCGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGCACTGGAGATCATGGACGGTAAATTTACCGGCACTGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCATCGCGCAGGGGTATGAAATCCCGCTGgcgacgaccgg                                                                                 >  2:491361/1‑243 (MQ=255)
                                                                                                                                                                            cAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGGAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattgg                                                                        >  2:9054/1‑251 (MQ=255)
                                                                                                                                                                                 ggTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGGAAATTTACCGGCAATGAGATCGGCGCCATCGTAGACGCGCAGTACAAATCGAAAAATCTGACACGCCTCGCGCCGGAGTATGCAACCCCGCCGgcgccgcgcgaggcgattagcgagg                                                                   >  2:408774/1‑230 (MQ=255)
                                                                                                                                                                                  ggggcgcaggCTGGAACCGCGGGGCTGGAAGGGGGCGATGGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATg                                                                   <  2:108666/242‑1 (MQ=255)
                                                                                                                                                                                     cTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGc                                                               >  1:261114/1‑251 (MQ=255)
                                                                                                                                                                                         aaGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCa                                                             >  1:384471/1‑250 (MQ=255)
                                                                                                                                                                                         aaGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                       >  1:136438/1‑240 (MQ=255)
                                                                                                                                                                                         aaGCGGGAAACGCTGGGCGGGAAAGGGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                       <  2:136438/240‑1 (MQ=255)
                                                                                                                                                                                          agctggaaacgttgggcgagaggtgggCGATGGCCTCCGGCGGCTTTCCTTGCTTTGCTAAATCTCTGCGCGGCAAGTTGCGCCTGACCGCCGGGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAATCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAAGCCCGCTGGCGCAGACCGTGGCGATTGGCTATGGAGCCAAt                                                           <  2:284756/226‑1 (MQ=255)
                                                                                                                                                                                           gCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATg                                                          <  1:362486/251‑1 (MQ=255)
                                                                                                                                                                                              ggAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGa                                                         <  1:478142/249‑1 (MQ=255)
                                                                                                                                                                                               gAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGAGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGAACGTGGCGATTGGCGATGGAGCCAATAAcc                                                       >  1:220057/1‑250 (MQ=255)
                                                                                                                                                                                                      tgggctggAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAACTGCCgatg                                               >  1:498824/1‑251 (MQ=255)
                                                                                                                                                                                                                              ccGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAGAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGACATCCCGCAGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCGGGTCTGGGGAtt                       >  2:8860/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGACAATAACCTGCAGATGATCAAAGCGGAAGGGCGgggggatgccaaccc               >  2:447705/1‑240 (MQ=255)
                                                                                                                                                                                                                                        tACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATg             >  1:176259/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ctttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCCTgcc           >  2:231920/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttctttGCTGAATAGCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgc            <  1:151433/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa    >  1:103094/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     gAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCaaaa  <  1:346183/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
CCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625286‑4625776

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