Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 619,358 G→A D369N (GAT→AAT)  fepE → polysaccharide co‑polymerase family protein FepE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,3580GA98.5% 211.0 / NA 67D369N (GAT→AAT) fepEpolysaccharide co‑polymerase family protein FepE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/31);  minor base C (1/0);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACTGACGTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACC  >  NC_000913/619114‑619603
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ggTCACTGACGTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGcc                                                                                                                                                                                                                                                  >  1:468723/1‑250 (MQ=255)
ggTCACTGACGTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGAGGGGGGCTTGTGGTGGCGTGCTGTTGGGCTATGCGATGGTATCCAGAAAAAAGCATGccc                                                                                                                                                                                                                                                 >  1:154468/1‑250 (MQ=255)
  tCACTGACGTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGGGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGGTGGCATCCAGAAAAAAGAATGCCatg                                                                                                                                                                                                                                               >  1:254844/1‑251 (MQ=255)
          gTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGc                                                                                                                                                                                                                                          >  1:40129/1‑248 (MQ=255)
               ggAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCACTTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATACAGAAAACAGAATGCCATGATGGACGACCACt                                                                                                                                                                                                                                  >  1:229370/1‑251 (MQ=255)
                     gAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATTCCATGATGGCAGACGGCTTAGtt                                                                                                                                                                                                                             <  2:327018/250‑1 (MQ=255)
                        cGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAAc                                                                                                                                                                                                                         <  1:434464/251‑1 (MQ=255)
                                tACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGGGGCGTGGTGTTGCGCTATGCGCTGGCATCCAGAAAACAGAATGACATGATGGCAGACCACTTAGTTTAACTCTTTGCg                                                                                                                                                                                                                 >  2:220348/1‑251 (MQ=255)
                                        cGGCAGCATCTTGTCGCGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAACGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTTTGGTGGCGTGCTGGTGCGCTCTGCGGAGGCATTCAGCAAACACAATGCCCTGATGGCTGACCCCTTCGTTTAAATATTTTCGGtttcagg                                                                                                                                                                                                         >  2:542858/1‑246 (MQ=255)
                                           cAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGTGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGc                                                                                                                                                                                                      <  2:69114/251‑1 (MQ=255)
                                                   tGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTc                                                                                                                                                                                                 <  1:453600/248‑1 (MQ=255)
                                                   tGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCg                                                                                                                                                                                                <  2:74618/249‑1 (MQ=255)
                                                          cAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGGGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCgg                                                                                                                                                                                       >  1:159204/1‑251 (MQ=255)
                                                                  aaaaGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACcaca                                                                                                                                                                                >  1:174767/1‑250 (MQ=255)
                                                                  aaaaGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACaa                                                                                                                                                                               <  1:30956/251‑1 (MQ=255)
                                                                     aGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcg                                                                                                                                                                            <  2:55287/251‑1 (MQ=255)
                                                                     aGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGGTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcg                                                                                                                                                                            >  1:444904/1‑251 (MQ=255)
                                                                     aGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGGGCTGTTGCGCTATGCGCTGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTTGTTTAACTATTTGCGGGTGAAGGTGCCGTTCGTTCGAGCGGccccccaatcc                                                                                                                                                                            >  2:264318/1‑241 (MQ=255)
                                                                       cacaTGTCCAAGATGTGACTTTTACGCCGTTTACGTCTCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGAGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGCGGCGCGCTGTTGCGCTATCCGATGGCATCCAGAAAAAAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCCCCACAAgcggg                                                                                                                                                                          >  2:267624/1‑250 (MQ=255)
                                                                       cacaTGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAACAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGGGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTTAAGGTTCCGTTCGTCCGAGCGGCACCACAAgcggc                                                                                                                                                                          >  2:495509/1‑251 (MQ=255)
                                                                            gTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTcc                                                                                                                                                                     <  2:294540/251‑1 (MQ=255)
                                                                            gTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTcc                                                                                                                                                                     <  1:563743/251‑1 (MQ=255)
                                                                            gTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGGGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTcc                                                                                                                                                                     >  1:401031/1‑251 (MQ=255)
                                                                                       gAATTTTACGCCGTTTACATATCAGTTAAGCCCGTCATTGCCAGTGAAAAACGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGGGCCGTTCGTCCGCGCGGCACCACAAACGGGGTTCCGGCAACTGGAt                                                                                                                                                          >  2:37396/1‑251 (MQ=255)
                                                                                        aaTTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATc                                                                                                                                                         <  1:219412/251‑1 (MQ=255)
                                                                                          ttttACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGGTGGTGGCTTGGGGTGGAGGGCTGGTGTGCTATGCGATGGAATTCAGGAAACAAAATGCCCAGATGGCAGACCACTTAGGTTAACTCTTTGCGGTTGAAGGTGCCGTTCGTCCGGTCGGCACCACAAGCGGCGTTCCGGCAACTgggatgg                                                                                                                                                       >  2:513558/1‑246 (MQ=255)
                                                                                          ttttACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGAAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCCCAAGCGGCGTTCCGGCAACTGGATCGt                                                                                                                                                       >  1:498093/1‑251 (MQ=255)
                                                                                                gCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATg                                                                                                                                                  >  1:300908/1‑250 (MQ=255)
                                                                                                        aaTATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAg                                                                                                                                         <  1:146597/251‑1 (MQ=255)
                                                                                                          tatCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGGGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcg                                                                                                                                       >  1:424530/1‑251 (MQ=255)
                                                                                                           atCAGTTAAGCCCGTCATTGCCCGTGAAACAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGGTGGTGGGTTGGGGTGGGGGGCTGTTGCGCTATGCGGTGGCATCCAGGAAACAGAATGCCATGGTGGCAGAACACTTAGTTTAACTATTTGCGGGTGAAGGGGCCGTTCGTCCGAGCGGCAACACACAAGGCGTTCCGGCAACCGGagcggcaaggaacaagccacgc                                                                                                                                      >  2:562938/1‑230 (MQ=255)
                                                                                                           atCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcgc                                                                                                                                      <  2:373959/251‑1 (MQ=255)
                                                                                                                   aGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCAt                                                                                                                               >  1:225772/1‑250 (MQ=255)
                                                                                                                    gCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                                                             <  1:93483/251‑1 (MQ=255)
                                                                                                                    gCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCCATGCGATGGGATCCAGAAAACAGAATGCCATGTTGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGGTACGTTCGGCCGGGCGGCACCCCAAGCGGCGTTCCGGCAACTGGATCGTCAATGGTCCTGCAgcccaaacccta                                                                                                                              >  2:159802/1‑240 (MQ=255)
                                                                                                                     cccGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGGGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGATTGCATCCAGAAAACAGAATGGCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAAGTGCCGTTCGGCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAAAGATAATGCAGCGCAGACCATaa                                                                                                                             >  2:2263/1‑250 (MQ=255)
                                                                                                                       cGTCAATGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTGCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGCCCATAAATg                                                                                                                          <  2:481705/251‑1 (MQ=255)
                                                                                                                         tCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGAGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACCGAATGCCATGATGGCAGACCACTTAGTTTTACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCCGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCCATGATCATGCAGCGCAGACCATAACTgcg                                                                                                                        >  1:368894/1‑251 (MQ=255)
                                                                                                                           aTTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGGTGGAGGCTTGGGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCt                                                                                                                      >  2:133901/1‑251 (MQ=255)
                                                                                                                                    gAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGtt                                                                                                             <  2:255941/251‑1 (MQ=255)
                                                                                                                                     aaaaaaGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGGGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                                                            >  1:271707/1‑251 (MQ=255)
                                                                                                                                       aaaaGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcag                                                                                                          <  2:185605/251‑1 (MQ=255)
                                                                                                                                           gACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCCTGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCCCAAGCGGCGTTCCGGCAACTGGATTGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGt                                                                                                      >  1:313834/1‑251 (MQ=255)
                                                                                                                                               gTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGTCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACa                                                                                                  <  1:196535/251‑1 (MQ=255)
                                                                                                                                                     gTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGGGGGGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCCCTTAGTTTAAATATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGGTCCGGCCACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTcc                                                                                            >  1:252898/1‑251 (MQ=255)
                                                                                                                                                       aaGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCtt                                                                                          <  1:220348/251‑1 (MQ=255)
                                                                                                                                                           cGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgc                                                                                      <  2:401031/251‑1 (MQ=255)
                                                                                                                                                            gATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGGGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCCCTTAGTTTAACTATTTGCGGTTGAAAGGGCCGTTCGTCCGAGCGGCCCCCCCAGCGGCGTTTCTGCCACTGGATCGTCAATGATCATTCAGCGCAGACCCTAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCgccg                                                                                     >  2:199651/1‑248 (MQ=255)
                                                                                                                                                                    tgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGTGGGCTGTTGCTCTATTCGATTGTATCCCGGAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCCGTCGGCCGAGCGGCACCACAAACGGCGTTCCGGCAACTGGATCGTCAATGATCATGCCGCGCAGACCATACATGCGCTCCATCCGTTACGCAGTGACCATCTCCTTCGgcgcctcctg                                                                              >  2:415809/1‑244 (MQ=255)
                                                                                                                                                                           ttCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaaca                                                                       <  2:448896/250‑1 (MQ=255)
                                                                                                                                                                               gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGGGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtttt                                                                  >  1:572319/1‑251 (MQ=255)
                                                                                                                                                                                 gTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                                                                <  1:132667/251‑1 (MQ=255)
                                                                                                                                                                                    gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttccctt                                                             >  1:437420/1‑251 (MQ=255)
                                                                                                                                                                                    gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGGTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCGCCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCAGGCAGCGCGGACCATAAATGCGCTCAATCAGTGCAGCAGTGACAGGCTCCTTCGGCGTGCCTTGAGCAAAAATTgtccct                                                              <  2:78405/250‑1 (MQ=255)
                                                                                                                                                                                     aTCGGCGGGATGGTGGCTTGAGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttc                                                            <  1:83867/251‑1 (MQ=255)
                                                                                                                                                                                       cggcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttccc                                                          <  2:250227/251‑1 (MQ=255)
                                                                                                                                                                                        ggcggGATGGTGGCTTGTGGAGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCACGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACCATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCg                                                         >  2:1210/1‑251 (MQ=255)
                                                                                                                                                                                              aTGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaattt                                                   <  2:183621/251‑3 (MQ=255)
                                                                                                                                                                                                              gCGTGCTGTTGCGCTATGCGATGGCATGCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTGAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCGCAAGCGGCGTTCCGGCAACTGGATCGGCAATGATCATGCAGCGCAGACCATAAATGCGCGCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTGCCTGAGCAGCAATTTGCCCTTCCCGCAATGCAAGCAAAGGGCTGGCg                                   <  1:276630/251‑1 (MQ=255)
                                                                                                                                                                                                              gCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCg                                   <  1:102005/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTGCGATGGCATCCAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCt                      <  1:175336/249‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGATGGCATCCAGACAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGAt                   >  1:464725/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ccAGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTgc        >  1:24355/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aGAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTgcag      <  2:174767/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gAAAACAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTgcagc     <  1:61763/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaCAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCAcc  <  1:59157/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaCAGAATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCAcc  >  1:515189/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
GGTCACTGACGTTGCGGAACTGAACGGTGAATTACGTAATCGGCAGTATCTTGTCGAGCAATTAACAAAAGCACATGTCAACGATGTGAATTTTACGCCGTTTAAATATCAGTTAAGCCCGTCATTGCCAGTGAAAAAAGACGGTCCGGGTAAGGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACC  >  NC_000913/619114‑619603

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