Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,377 T→C R161R (CGT→CGC)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,3770TC96.1% 145.5 / ‑6.3 51R161R (CGT→CGC) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (21/28);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGAC  >  NC_000913/4625133‑4625620
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ccggTTCACTCACCGCACGGGCTACGCGCCTGGCCCGAGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAGCCCGCCCCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCg                                                                                                                                                                                                                                                <  1:208754/249‑1 (MQ=255)
        ctcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGACTGGGCAGCCG‑CGAGGCGACGcggac                                                                                                                                                                                                                                        >  2:272965/1‑247 (MQ=255)
                              gggCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAACCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCCCCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTg                                                                                                                                                                                                                 <  1:161219/250‑1 (MQ=255)
                               ggCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CCATATTCTGc                                                                                                                                                                                                                >  2:481666/1‑251 (MQ=255)
                                gcccatgaggcgcacctggagtgcGCCCCTCGGGGGAAAACCCCGCACCTGCGCACGCGGGGTTGGCTGGTGATGGACATGGTTTCCACCGCCATCCAGATTGAAGGTATTGATGAAATTGCCAACCTGGCCGGAACGGGCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCa                                                                                                                                                                                                               <  2:148632/228‑1 (MQ=255)
                                   catgaagcgcaatgggaggTCGCCCCGCTGGGGAAAATCCCGCACCGGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCCCCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACGCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACa                                                                                                                                                                                                            <  2:215092/233‑1 (MQ=255)
                                           gcaacTGGAGGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTg                                                                                                                                                                                                    <  1:457862/247‑1 (MQ=255)
                                                          cccGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTg                                                                                                                                                                                      >  1:284756/1‑250 (MQ=255)
                                                                  ggAAAATACCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTCTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGc                                                                                                                                                                             >  1:329291/1‑251 (MQ=255)
                                                                    aaaaTCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCtt                                                                                                                                                                           >  2:132375/1‑251 (MQ=255)
                                                                        ccccGCACCAGCGCACGGGGGGTTTGCTGTTGATGGACTTGGATTCCACCGCCCTCCAGATTGAAGGTATTGATGAAATGACCAAACTGGCCGGAACGGGCGAGATGGTGGGGGAAGTAACCGAACGGGCGATGCGCGGCGAACCCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGTTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTCACg                                                                                                                                                                       <  2:93015/250‑1 (MQ=255)
                                                                            gCACCTGCGCACGCCGGGTTTGCTGGGGATGGACATGGATTCCACCGCCATCCAGATTGCATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGAGGCGGAAGTAACCGAACGGGCGATGCGCGCCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGACGGGCGCTGACGC‑CCATATTCTGCAACAGGGGCGTGAAAATCTTCCGCGGATGCCCGGCTTAACGCAAc                                                                                                                                                                   >  2:69211/1‑251 (MQ=255)
                                                                             ccccTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACt                                                                                                                                                                  <  1:490294/249‑1 (MQ=255)
                                                                              aCCTGCGAACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGGACGGGCGATGCGAGGCGAACTCGATTTTACCGCGGGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGACGCTGACGC‑CAATATTCTGCCAAAGGTACCTGAAAATCTGCCGCAGATGcccgcagtaaggcaacgg                                                                                                                                                                 >  2:15678/1‑235 (MQ=255)
                                                                                   ggcACGCCGGGTTGGCTGGTGATGGAGATGGATTCCACCGCCATCCAGATTTAATGTATTGATGAAATGGCAAACCGGGCGGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCt                                                                                                                                                            <  2:200911/250‑1 (MQ=255)
                                                                                     cACGACGGGTGTGCTGGTGATGGACCTGGATTCCACCGCCATCCAGATTGAATGAATTGATGAAATTGCCAAACTGGCCGGCACGGGCGAGGTTGGGGCGGAAGTAACCGAACGAGCGATGCGCAGCGAACTCGATTTTACCGCCCGACTGCGCAGCCG‑CGTGGCGACGCTGAAAGGAGCTGACGC‑CAATATTCTGCAACAGGTGCGTGACAATCTGCAGCTGATGCCAGGcatacctcgactggtacac                                                                                                                                                           >  2:75997/1‑232 (MQ=255)
                                                                                                 gCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgc                                                                                                                                              >  2:412162/1‑251 (MQ=255)
                                                                                                        aTGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggctg                                                                                                                                       >  1:380973/1‑251 (MQ=255)
                                                                                                                 gattCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGACCGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGc                                                                                                                              <  1:112432/248‑1 (MQ=255)
                                                                                                                 gattCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGc                                                                                                                              <  1:69579/248‑1 (MQ=255)
                                                                                                                     ccaccgcatcacAGATTGAATGTTTTGAAGAAATCGCCAATCGGGCGGGAAGGGGCGAGAGGGGGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCTATTTTTCCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAAGCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGAtt                                                                                                                          <  2:126360/240‑1 (MQ=255)
                                                                                                                     ccACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggctgagatcggaagagc                                                                                                                          <  1:481666/251‑14 (MQ=255)
                                                                                                                        ccGCCCTCCAGATTGAATTGAGTGATGAAATGGCCAAGTGGGCCGGAACGGGCGAGATGGTGGCGGAAGTACCCGTACGGGCGATGCGCGGCGAACTCGATTTTTCCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGcc                                                                                                                       <  2:491286/251‑1 (MQ=255)
                                                                                                                           ccATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTcc                                                                                                                    <  1:102668/251‑1 (MQ=255)
                                                                                                                             aTCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcgg                                                                                                                  >  1:132516/1‑251 (MQ=255)
                                                                                                                                   aTTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTACCGCTGATGCCAGGCTTAACACAACTGGTGCTCAAGCTGGAAACGCGGGGCTGACAAGTGGCGAGTGACTCCGGCGGCtt                                                                                                            >  1:15323/1‑251 (MQ=255)
                                                                                                                                            aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttctt                                                                                                   >  1:378895/1‑251 (MQ=255)
                                                                                                                                              tgatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttcttt                                                                                                  <  1:154791/250‑1 (MQ=255)
                                                                                                                                                atgaAATTGCCAACCTGGCCGGCACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTCACGC‑CAATATTCTGCAACAGGTGCGTGAACATCTGCCGCTGCTGCCAGGCTTAACGCACCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGGAGCGGTTGCATCCGACGGCTTTACTTTCTTTGCt                                                                                               >  2:210669/1‑251 (MQ=255)
                                                                                                                                                  gaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTg                                                                                              <  1:383030/250‑1 (MQ=255)
                                                                                                                                                   aaaTTGCCAAATGGGCCGGAGCGGGCGGGAGGGTGGCGTAAGGAACCGAAGGGGCGAGGCGCGGCGAACTCGATTGTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTGAACGCAACTGGTGCTCAATCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGaa                                                                                            <  1:43485/251‑1 (MQ=255)
                                                                                                                                                        gCCAAACTGGCCGGAACGGGCGATATGGTGGCGGAAGTACCCAACCGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATatct                                                                                       <  1:452312/251‑4 (MQ=255)
                                                                                                                                                         ccAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTg                                                                                      <  1:185881/251‑1 (MQ=255)
                                                                                                                                                          aaaaCTGGCCGGAACGGGGGAGAGGGGGGGGGAAGTACCCGAACGGGGGAGGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTgc                                                                                     <  2:60887/250‑1 (MQ=255)
                                                                                                                                                                        cGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTgcgc                                                                       >  1:245572/1‑251 (MQ=255)
                                                                                                                                                                         gggCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGcc                                                                      <  1:449862/251‑1 (MQ=255)
                                                                                                                                                                             gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGAc                                                                  >  2:189852/1‑251 (MQ=255)
                                                                                                                                                                                         gaagaaccctaaCGGGCGATGCGCGGCGACTTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATCCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCGGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGc                                                      <  1:485913/241‑1 (MQ=255)
                                                                                                                                                                                                      gggggattcgcgggcaacccGATTTTACCGCAAGCCTGCGCAGCGG‑TGTGGCGACGCTGAAAGGCGCTGCCGC‑CAATATTCTGCAACAGGTGCGAGAAAATCGGCCGCTGATGCCATGCTTAACGCATCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGgag                                         <  2:312086/232‑1 (MQ=255)
                                                                                                                                                                                                      gggCGAGGCGCGGAAACCTCGTTTTTCCCGCCAGCCGGCGCAGCCG‑CGTGGCGACGCGGAAATGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTACCGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCTGTGGCCTCCGGCGGCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGgag                                         <  2:33711/251‑1 (MQ=255)
                                                                                                                                                                                                         cGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                           >  1:92325/1‑246 (MQ=255)
                                                                                                                                                                                                              cgcgGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTATGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATgg                                  >  1:22900/1‑250 (MQ=255)
                                                                                                                                                                                                                 ggcgAACTCGTTTTTCCCGCCAGCCTGCGCAGCGGCTG‑GGCGACGCTGAAAGGCGTTGA‑GCACAATATTCTGAAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGTTTAACGCAACTGGTGCCTAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGTTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACgg                              <  2:436564/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggcgAACTCGATTTTCCAGCCAGCCGGCGCAGCCG‑CGTGGCGACGCGGAAAGGCGTTGACGC‑CAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCGGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCCGAAGATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCgg                              <  2:153471/251‑1 (MQ=255)
                                                                                                                                                                                                                    gAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACggaa                            >  1:372213/1‑248 (MQ=255)
                                                                                                                                                                                                                         cGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGGGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAtt                        <  1:321889/249‑1 (MQ=255)
                                                                                                                                                                                                                          gtttttACCGCCAGCCTGCGAAGCCG‑CGTGGCGACGCTGAAAGTGGCGGACGC‑CAATTTTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAc                     <  2:157111/249‑1 (MQ=255)
                                                                                                                                                                                                                           aTTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                    >  2:346183/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaCCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                 <  1:410689/249‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGCCAGTCTGCGCAGCCT‑CGGGGCGACGCGGAAAGGCGCTGACGC‑CAATATCCTGCAACAGGTGCGTGAAACTCTGCCGATTATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCTCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAt              <  1:312712/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gcAGCCG‑CGTGGAGACGATGAAAGGAGCTGACGC‑CCATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGATGGAAAGTGGCGATTGCCTCCGGCGGATTTACTTTCTTTGCTGAATATCTGCGCGAACGGTTGAGCCTGAACGCAGTGGTAGACAATGAACTGGGGATCATGGACGGTAAATTTCCCGGGAATGTGAGCGGGGAc  >  2:276662/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGAC  >  NC_000913/4625133‑4625620

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