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→CGCserBphosphoserine 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

TGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCG  >  NC_000913/4625130‑4625618
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
tGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAACCTGCGCAGCCGcgtg                                                                                                                                                                                                                                                >  1:444659/1‑251 (MQ=255)
   cAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCg                                                                                                                                                                                                                                             >  1:473570/1‑251 (MQ=255)
       ttCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCACGCTGGGGAAAATCCCGCACCTGCGAACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGGTTGCATGTATTGATGTAATTGCCAAACAGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGc                                                                                                                                                                                                                                         >  1:37217/1‑251 (MQ=255)
          actcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTg                                                                                                                                                                                                                                       >  1:470659/1‑250 (MQ=255)
             cacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTTGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAGTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAg                                                                                                                                                                                                                                   <  1:273060/251‑1 (MQ=255)
                     gggCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGAc                                                                                                                                                                                                                           <  1:211941/251‑1 (MQ=255)
                                    cccATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa                                                                                                                                                                                                             >  1:441240/1‑250 (MQ=255)
                                    cccACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa                                                                                                                                                                                                             >  2:131327/1‑250 (MQ=255)
                                      catgAAGCGAAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACa                                                                                                                                                                                                          <  1:461755/248‑1 (MQ=255)
                                         gaatccaacaggtgtttCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGTTGATGGACATGGATTCCACTGCCATCAAGAGTGAATGTATTGATGAAATTGCCAAACTGGCAGGACCGGGTGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCCGGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGGTGACGCCAATATTCTGCAACAGGt                                                                                                                                                                                                       <  2:263774/235‑1 (MQ=255)
                                              gcaacTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTg                                                                                                                                                                                                  <  1:269002/247‑1 (MQ=255)
                                                  cTGGAAGTCGCCCGGGTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaaa                                                                                                                                                                                               >  1:680236/1‑250 (MQ=255)
                                                   tGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGACATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAAt                                                                                                                                                                                             >  1:115669/1‑251 (MQ=255)
                                                    ggATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCTCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCTAACTCGATTTTGCCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaa                                                                                                                                                                                                <  1:95409/247‑1 (MQ=255)
                                                     gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCt                                                                                                                                                                                           >  1:414571/1‑251 (MQ=255)
                                                      aTGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTg                                                                                                                                                                                          <  1:572740/251‑1 (MQ=255)
                                                          cGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCACAAGGAGCGTGAAAATCTGCCGc                                                                                                                                                                                      >  1:148680/1‑251 (MQ=255)
                                                              cctcgggggAAACCCCCGCCCCTGCGCACGCCGGGTTGGCGGGTGAGGGACATGGTTTCCACCCCCATCCAGATTGAATGGATTGAGGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGAt                                                                                                                                                                                  <  2:341839/246‑1 (MQ=255)
                                                                 gtGGGAAAAATCCCGCTCCTGCCCGCGCCGGGTTTGTGGGTGATGAACAGGGATTCCTCCGCCAACCAGATTGAATGTATGGATGAAATGGCCAACCTGGCCGGAACGGGGGAGATGGTGGCGTAAGTAACCGAACGGGCGAGGCGCGGCGACCTGGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGcc                                                                                                                                                                               <  2:91454/250‑1 (MQ=255)
                                                                   ggggAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTACCGCCAATATTCTGCAACAGGTGCGTGCAAATCTGCCGCTGATGcc                                                                                                                                                                               >  2:535579/1‑249 (MQ=255)
                                                                       aaaaCCCCGCCCCTGCGCCCGCCGGGTTTGCGGGTGATGGACAGGGATTCCACCGCCATCCAGATTGAAGTTATTGATGAAATTGCCAACCTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCtt                                                                                                                                                                         <  1:365558/251‑1 (MQ=255)
                                                                           cccaatcccggcgcCCGCCGGGTTTGCTGGTGAGGGACAGGTATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGACCGGGCGATGCACGGCGAACCCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACg                                                                                                                                                                     <  2:618006/241‑1 (MQ=255)
                                                                                  ccTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTgg                                                                                                                                                              <  1:223462/251‑1 (MQ=255)
                                                                                   cTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCACGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACg                                                                                                                                                                     >  1:618006/1‑243 (MQ=255)
                                                                                        ccctctgggTTTGCTGGTGATGGATATGGACGCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACGGGCAGGAATGGGCGAGATGGTGGCGGAAGTACCCGAACGGGCGATGCGCGGCAAACCCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGTAACAGGTGCGTGAAAATCTGCGGCTGATGCCAGGCTTAACGCAACTGGTGCTCa                                                                                                                                                        <  1:226429/245‑1 (MQ=255)
                                                                                        cACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGGATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGCGATGGTGGCGGAAGTAACCGACCGGACGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCTCGGGGCGACGATGGAAGACGCTGACGCCAATATTATGCAACAGGGGCGTGAAATTCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTc                                                                                                                                                         >  2:493669/1‑250 (MQ=255)
                                                                                        cACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAACGGCGCTGACGCCAAAATTCTGCAACAGGGGCGTGAACATCTGCCGCTGATGCCCGGCTTAACGCaaagggttctc                                                                                                                                                         >  2:44205/1‑241 (MQ=255)
                                                                                                       ggggATGGACATGGATTCCACCGCCATCCAGATTGAATGTTTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTTCAACAGGTGCGTGAAAATCTGCCGCGGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgct                                                                                                                                           <  1:204293/246‑1 (MQ=255)
                                                                                                        gTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCAAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctggg                                                                                                                                        <  1:4910/251‑1 (MQ=255)
                                                                                                                     attCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCg                                                                                                                           <  1:116378/249‑1 (MQ=255)
                                                                                                                           ccGCCATCCAGAGTGAATGTATTGATGAAATTTCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGcc                                                                                                                     <  1:401571/251‑1 (MQ=255)
                                                                                                                                aTCCAGATTGAATGTATTGATGAAATTGCAAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcgg                                                                                                                >  1:202756/1‑251 (MQ=255)
                                                                                                                                  ccAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcggc                                                                                                               <  1:131327/250‑1 (MQ=255)
                                                                                                                                  ccAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGGGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCACCTGGTGCTCAAGCTGGAAACGCAGGGCTGGAAAGTGGCGATTGCCTCcgtcg                                                                                                              >  2:656269/1‑248 (MQ=255)
                                                                                                                                   taGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGGGCTCAAGCTGGAAACGCTGGGCTGGAAAGGGGCGATTGCCTCcggcgg                                                                                                             >  2:254102/2‑251 (MQ=255)
                                                                                                                                   cAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCcggcgg                                                                                                             >  2:441354/1‑251 (MQ=255)
                                                                                                                                       atgtatgttatttaggaagtggcAAACTTGGCCGGAACGGGTGAGATGTTGGCGGAAGTACCGGAACGGGGGAGGCGCGGCGAACTCTATTTTACCGCCAGCCTGCGCAGCCGCGTGTCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGAGGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCt                                                                                                           <  2:582855/228‑1 (MQ=255)
                                                                                                                                        taaattatttGATGAAATGGCCAAATTGGCCGGACCGGGGGAGAGGGTGGCGGAAGTACCCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCttt                                                                                                         <  2:114301/242‑1 (MQ=255)
                                                                                                                                               aTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttctt                                                                                                 >  1:306397/1‑251 (MQ=255)
                                                                                                                                                ttGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGAGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTActttcttt                                                                                                >  1:35377/1‑251 (MQ=255)
                                                                                                                                                  gatgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGGGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGc                                                                                              <  1:535579/251‑1 (MQ=255)
                                                                                                                                                    tgaAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAAACGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGAGGCTGAAAGGCGCTGCCGCCCATATTCTGCAACCGGTGCGTGACAATCTGCCGCTGCTGCCAGGCTTAACGCAACTGGTGATCAGGCGGGAAACGATGGGCAGGAAAGTGGAGAATGCCTCCGGAGGCATTACCTTCATTgttt                                                                                            >  2:56611/1‑248 (MQ=255)
                                                                                                                                                          tGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGGTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGgcgc                                                                                                                                                                                                                               <  2:228600/114‑1 (MQ=255)
                                                                                                                                                          tGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGgcgc                                                                                                                                                                                                                               >  1:228600/1‑114 (MQ=255)
                                                                                                                                                           gCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATatc                                                                                      <  1:441354/250‑3 (MQ=255)
                                                                                                                                                               aaCTGGCCGGAACGGGCGAGATGGTGGCGGAAGAAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGCAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGATTTCCTTTCTTTGCTGAATATCTgagc                                                                                 >  2:457544/1‑248 (MQ=255)
                                                                                                                                                                aCTGGCCGGAACGGGGGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCAAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcg                                                                                <  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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
TGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCG  >  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: