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,3770TC95.8% 207.3 / ‑3.7 73R161R (CGT→CGC) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (37/32);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

TGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCG  >  NC_000913/4625130‑4625618
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
tGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAACCTGCGCAGCCGcgtg                                                                                                                                                                                                                                                >  1:444659/1‑251 (MQ=255)
   cAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCg                                                                                                                                                                                                                                             >  1:473570/1‑251 (MQ=255)
       ttCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCACGCTGGGGAAAATCCCGCACCTGCGAACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGGTTGCATGTATTGATGTAATTGCCAAACA‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGc                                                                                                                                                                                                                                         >  1:37217/1‑251 (MQ=255)
          actcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTg                                                                                                                                                                                                                                       >  1:470659/1‑250 (MQ=255)
             cacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTTGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAGTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAg                                                                                                                                                                                                                                   <  1:273060/251‑1 (MQ=255)
                     gggCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGAc                                                                                                                                                                                                                           <  1:211941/251‑1 (MQ=255)
                                    cccATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa                                                                                                                                                                                                             >  1:441240/1‑250 (MQ=255)
                                    cccACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa                                                                                                                                                                                                             >  2:131327/1‑250 (MQ=255)
                                      catgAAGCGAAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACa                                                                                                                                                                                                          <  1:461755/248‑1 (MQ=255)
                                         gaatccaacaggtgtttCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGTTGATGGACATGGATTCCACTGCCATCAAGAGTGAATGTATTGATGAAATTGCCAAACT‑GGCAGGACCGGGTGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCCGGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGGTGACGCCAATATTCTGCAACAGGt                                                                                                                                                                                                       <  2:263774/235‑1 (MQ=255)
                                              gcaacTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTg                                                                                                                                                                                                  <  1:269002/247‑1 (MQ=255)
                                                  cTGGAAGTCGCCCGGGTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaaa                                                                                                                                                                                               >  1:680236/1‑250 (MQ=255)
                                                   tGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGACATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAAt                                                                                                                                                                                             >  1:115669/1‑251 (MQ=255)
                                                    ggATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCTCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCTAACTCGATTTTGCCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaa                                                                                                                                                                                                <  1:95409/247‑1 (MQ=255)
                                                     gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCt                                                                                                                                                                                           >  1:414571/1‑251 (MQ=255)
                                                      aTGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTg                                                                                                                                                                                          <  1:572740/251‑1 (MQ=255)
                                                          cGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCACAAGGAGCGTGAAAATCTGCCGc                                                                                                                                                                                      >  1:148680/1‑251 (MQ=255)
                                                              cctcgggggAAACCCCCGCCCCTGCGCACGCCGGGTTGGCGGGTGAGGGACATGGTTTCCACCCCCATCCAGATTGAATGGATTGAGGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGAt                                                                                                                                                                                  <  2:341839/246‑1 (MQ=255)
                                                                 gtGGGAAAAATCCCGCTCCTGCCCGCGCCGGGTTTGTGGGTGATGAACAGGGATTCCTCCGCCAACCAGATTGAATGTATGGATGAAATGGCCAACCT‑GGCCGGAACGGGGGAGATGGTGGCGTAAGTAACCGAACGGGCGAGGCGCGGCGACCTGGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGcc                                                                                                                                                                               <  2:91454/250‑1 (MQ=255)
                                                                   ggggAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTACCGCCAATATTCTGCAACAGGTGCGTGCAAATCTGCCGCTGATGcc                                                                                                                                                                               >  2:535579/1‑249 (MQ=255)
                                                                       aaaaCCCCGCCCCTGCGCCCGCCGGGTTTGCGGGTGATGGACAGGGATTCCACCGCCATCCAGATTGAAGTTATTGATGAAATTGCCAACCT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCtt                                                                                                                                                                         <  1:365558/251‑1 (MQ=255)
                                                                           cccaatcccggcgcCCGCCGGGTTTGCTGGTGAGGGACAGGTATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGACCGGGCGATGCACGGCGAACCCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACg                                                                                                                                                                     <  2:618006/241‑1 (MQ=255)
                                                                                  ccTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTgg                                                                                                                                                              <  1:223462/251‑1 (MQ=255)
                                                                                   cTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCACGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACg                                                                                                                                                                     >  1:618006/1‑243 (MQ=255)
                                                                                        ccctctgggTTTGCTGGTGATGGATATGGACGCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACG‑GGCAGGAATGGGCGAGATGGTGGCGGAAGTACCCGAACGGGCGATGCGCGGCAAACCCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGTAACAGGTGCGTGAAAATCTGCGGCTGATGCCAGGCTTAACGCAACTGGTGCTCa                                                                                                                                                        <  1:226429/245‑1 (MQ=255)
                                                                                        cACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGGATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGCGATGGTGGCGGAAGTAACCGACCGGACGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCTCGGGGCGACGATGGAAGACGCTGACGCCAATATTATGCAACAGGGGCGTGAAATTCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTc                                                                                                                                                         >  2:493669/1‑250 (MQ=255)
                                                                                        cACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAACGGCGCTGACGCCAAAATTCTGCAACAGGGGCGTGAACATCTGCCGCTGATGCCCGGCTTAACGCaaagggttctc                                                                                                                                                         >  2:44205/1‑241 (MQ=255)
                                                                                                       ggggATGGACATGGATTCCACCGCCATCCAGATTGAATGTTTTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTTCAACAGGTGCGTGAAAATCTGCCGCGGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgct                                                                                                                                           <  1:204293/246‑1 (MQ=255)
                                                                                                        gTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCAAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctggg                                                                                                                                        <  1:4910/251‑1 (MQ=255)
                                                                                                                     attCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCg                                                                                                                           <  1:116378/249‑1 (MQ=255)
                                                                                                                           ccGCCATCCAGAGTGAATGTATTGATGAAATTTCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGcc                                                                                                                     <  1:401571/251‑1 (MQ=255)
                                                                                                                                aTCCAGATTGAATGTATTGATGAAATTGCAAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcgg                                                                                                                >  1:202756/1‑251 (MQ=255)
                                                                                                                                  ccAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcggc                                                                                                               <  1:131327/250‑1 (MQ=255)
                                                                                                                                  ccAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGGGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCACCTGGTGCTCAAGCTGGAAACGCAGGGCTGGAAAGTGGCGATTGCCTCcgtcg                                                                                                              >  2:656269/1‑248 (MQ=255)
                                                                                                                                   taGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGGGGCGATTGCCTCcggcgg                                                                                                             >  2:254102/2‑251 (MQ=255)
                                                                                                                                   cAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcggcgg                                                                                                             >  2:441354/1‑251 (MQ=255)
                                                                                                                                       atgtatgttatttaggaagtggcAAACTTGGCCGGAACGGGTGAGATGTTGGCGGAAGTACCGGAACGGGGGAGGCGCGGCGAACTCTATTTTACCGCCAGCCTGCGCAGCCGCGTGTCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGAGGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCt                                                                                                           <  2:582855/228‑1 (MQ=255)
                                                                                                                                        taaattatttGATGAAATGGCCAAATT‑GGCCGGACCGGGGGAGAGGGTGGCGGAAGTACCCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCttt                                                                                                         <  2:114301/242‑1 (MQ=255)
                                                                                                                                               aTTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttctt                                                                                                 >  1:306397/1‑251 (MQ=255)
                                                                                                                                                ttGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGAGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttcttt                                                                                                >  1:35377/1‑251 (MQ=255)
                                                                                                                                                  gatgaAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGGGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGc                                                                                              <  1:535579/251‑1 (MQ=255)
                                                                                                                                                    tgaAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAAACGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGAGGCTGAAAGGCGCTGCCGCCCATATTCTGCAACCGGTGCGTGACAATCTGCCGCTGCTGCCAGGCTTAACGCAACTGGTGATCAGGCGGGAAACGATGGGCAGGAAAGTGGAGAATGCCTCCGGAGGCATTACCTTCATTgttt                                                                                            >  2:56611/1‑248 (MQ=255)
                                                                                                                                                          tGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGGTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGgcgc                                                                                                                                                                                                                               <  2:228600/114‑1 (MQ=255)
                                                                                                                                                          tGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGgcgc                                                                                                                                                                                                                               >  1:228600/1‑114 (MQ=255)
                                                                                                                                                           gCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATatc                                                                                      <  1:441354/250‑3 (MQ=255)
                                                                                                                                                               aaCT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGAAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGCAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGATTTCCTTTCTTTGCTGAATATCTgagc                                                                                 >  2:457544/1‑248 (MQ=255)
                                                                                                                                                                aCT‑GGCCGGAACGGGGGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCAAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcg                                                                                <  1:109859/251‑1 (MQ=255)
                                                                                                                                                                    ggCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAc                                                                              <  1:93958/250‑1 (MQ=255)
                                                                                                                                                                     gCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGAGTGACGACGCTGAAAGGCGCTGACGCCAATATTCAGCAACAGGCGCGTGAACATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcgtca                                                                             >  1:409534/1‑247 (MQ=255)
                                                                                                                                                                       cggcaCGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTTCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACa                                                                             >  2:5985/5‑248 (MQ=255)
                                                                                                                                                                       cGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAg                                                                           <  1:5985/250‑1 (MQ=255)
                                                                                                                                                                        ggAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGCCGCTGAAAGGCGCTGACGCCACTATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAATAGGCGATTGCCTCCGGCGGGTTTACATTCTTTGCTGattatatacgcgacatggt                                                                         >  2:337501/1‑233 (MQ=255)
                                                                                                                                                                        ggAACGGGCGAGATGGTGGCGGAAGTAACCGAAAGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCAGCGTGGCGAAGCTGAAAGGAGCTGACGCCAATATTATGCAAAAGGTGCGTGGAAATCTGCCGCTGCTGCCAGGCTTAACGAAACTGGTGCTCAAACTGGAAACGCTGGGGTGGAAAGCGGCGACTGCCTCAGGCGGGTTTACTTACTTTGGTGAATataggcgcgtcaaggt                                                                         >  2:450844/1‑236 (MQ=255)
                                                                                                                                                                        ggAACAGGCGAGATGGTAGCGGAAGTAACAGAAAGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCACAGACGCGTGGCGACGCTGAAAGGCGCTGGCGCCAATATTCTGAAACAGGTGCGTGAAAATATGCCGCTGATGACAGGCTTAACGCAACTGGTGCTCAAGCTGGAAATGCTGGGCTGGAAAGAGGCGCTGACCACCAGCGGATTTATTTTCTTTGCTGAATATCAGAGCGACAAgtt                                                                         >  2:433954/1‑249 (MQ=255)
                                                                                                                                                                            cGGGCGAGATGGTGGCGGAATTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGTCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTgcgc                                                                     <  1:65147/251‑1 (MQ=255)
                                                                                                                                                                              ggCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCt                                                                   >  1:552770/1‑251 (MQ=255)
                                                                                                                                                                                cGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAACGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGAGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATTTCTGCGGGCCAAGTTGCGcccgc                                                                 >  2:225306/1‑248 (MQ=255)
                                                                                                                                                                                  agaTGGAGGCGGAAGGAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTGCGCAGCCGCGTGGCGACGCTGACAGGCGCTGACGCCAATATTCAGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCGGGAAAGTGGCGCATGCCTCCGGCAGCTTTCCTTTCTTTGATGAATATCTGcgcggcgattttggccttacc                                                               >  2:544623/1‑234 (MQ=255)
                                                                                                                                                                                          gCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGACGCGAGGCGACGCTGAAAGGCGCTGACGCCACTATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGTATGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCGGCGCGCCAAGTTGCGCCTGACCGACGTgg                                                        >  2:526336/1‑250 (MQ=255)
                                                                                                                                                                                              aaGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGAGCGACAAGTTGCGCCTGACCGCCGTGGTAGcc                                                   >  1:70454/1‑251 (MQ=255)
                                                                                                                                                                                                   aCCGAACGGGCGATGCGCGGCGAACTCGAGTTTACCGCCAGCAAGCGCCACCGCGAGGCGCCGCTGACAGGCGCGGACGACAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGACACGCTGGGCGGGAAAGTGGCGATTGCCTACGGCGGCTTTACTTTCTTTGCTGAATAGCTGCGCGACAAGTTGCGCCTGACAGCCGTGGTGGCCAATgg                                              >  2:68265/1‑250 (MQ=255)
                                                                                                                                                                                                       aacggtcgaggcgggtcgAACTCGGTTTCACCGCCAGCCTGCGCAGCAGCGTGGCGACGCTGAAAGGCTATGACGCCAATATTCTGCAACAGGGTCGTTAAAATCTGCCGCTGATGCCAGGCTTAACCCAACTGGTGCTCAAGCTGGAAACGCTGGTCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCGGACCGCCGTGGTAGCCAATGAACTg                                          <  2:519704/235‑1 (MQ=255)
                                                                                                                                                                                                          gggCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGaa                                                                                                                                                >  1:328597/1‑146 (MQ=255)
                                                                                                                                                                                                          gggCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                         >  1:588152/1‑249 (MQ=255)
                                                                                                                                                                                                                     ggcgAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACgg                            >  1:357561/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaaCTCGTTTTTCCCGCCAGCCTGCGCAGCCGCGTGGCGCCGCTGAAAGGCGCTGACGCAAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTa                          <  2:357561/249‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCaa             <  1:544003/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cccacccttcgcAGCCGCGGGGCGACGCTGAAAGGCGTGGACGCCAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAtgt          <  2:202756/242‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGAt       >  1:317025/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATc      <  1:656269/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGCGCAGCCTCGTGGCGCAGCTGAAAGGCGCTGACGCCCATATTCTGCAAAAGGTGCGTGAAATTCTGCTGCTGATGACAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGCGATCATGGACGGTCACTTTACAGGCAATGTGAtaggg   >  2:463642/1‑247 (MQ=255)
                                                                                                                                                                                                                                               gcggaGCCGCGTGGCGACGCTGAAAGGCGCTGACTCCAATATGCGTCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTTCGCCTGACCGCTGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGc   <  1:382543/246‑1 (MQ=255)
                                                                                                                                                                                                                                               gcgcAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCgg    <  1:567451/249‑1 (MQ=255)
                                                                                                                                                                                                                                                cgcAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCg  >  1:335640/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
TGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACT‑GGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCG  >  NC_000913/4625130‑4625618

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