Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 615,736 A→C E527A (GAA→GCA)  entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913615,7360AC98.5% 211.5 / ‑6.1 65E527A (GAA→GCA) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (25/39);  minor base G (1/0);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGAGATTCCAGAAACCACGCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTT  >  NC_000913/615504‑615962
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
gTTGAGATTCCAGAAACCACGCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGGCCTACCGCGCTCGGTCTTTTTGAACCCGGCACTCCATGCGATAGTTGCAGc                                                                                                                                                                                                                                 >  1:193443/1‑236 (MQ=255)
     aatcccAGACACCCAGCTTAGCGCGCGTGTGGCAGACAGAGCGGCAAAATCACCGGATGCTCCGGCGCTGGCAGATGCGCGTTGCCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGa                                                                                                                                                                                                                             <  2:497385/231‑1 (MQ=255)
          cAGAACCCGCGCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGc                                                                                                                                                                                                                        <  1:179868/235‑1 (MQ=255)
                 cgctcTTGGCGCGCTGGTGCCACAACGAGCTGCAAAACCACCGGATGCTCCGGCGCTGGCAGATGGTCGTTCCCAGTTCGTGCTGCGGAAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCGGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGCCCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGc                                                                                                                                                                                                                   <  2:3680/229‑1 (MQ=255)
                   cGCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCg                                                                                                                                                                                                              <  1:538424/236‑1 (MQ=255)
                    gCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACCCCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGc                                                                                                                                                                                                             <  2:220415/236‑1 (MQ=255)
                                ggtggCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGGCTTTTTGACCCTGGCACTCCATGCGCTAGTTTGACATGGAGCGGCCTGGCTTCCCCTGGATACCGGc                                                                                                                                                                                                  >  1:674666/1‑235 (MQ=255)
                                ggtggCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGCTAAACCAGGGGACAGCGTGGCGGTGGCACTACCCCGCTCGGTCTTTTTGACCCTGGCACTCCATTCGATAGTTGCAGCTGGAGCGGCCTTGCTTCCGCTGGGTACCgcc                                                                                                                                                                                                  >  1:97238/1‑233 (MQ=255)
                                           aaGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGGTGCAGCTGGAGCGGCCTGGCTACCGCTGGATAACGGCCATCCGGACGa                                                                                                                                                                                       >  1:707580/1‑235 (MQ=255)
                                           aaGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGGCTTTTTGAACCTGGCACTCCATGCGATAGATGCACCTGGAGCGGCCTGGGTACCGCCGGGGACCGGCTAGTCGGACGa                                                                                                                                                                                       >  2:197500/1‑235 (MQ=255)
                                                                           gCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGa                                                                                                                                                                                                                             <  2:514096/165‑1 (MQ=255)
                                                                           gCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGa                                                                                                                                                                                                                             >  1:514096/1‑165 (MQ=255)
                                                                             agtggcagtttACCTGGTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGGGAGGGCGGCGTTAAACCAGGGGGAAACGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgt                                                                                                                                                    <  2:344421/227‑1 (MQ=255)
                                                                                   gtgTTGCCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGGGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGc                                                                                                                                                 <  2:547904/231‑1 (MQ=255)
                                                                                   gcgTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCg                                                                                                                                                  <  1:515309/232‑1 (MQ=255)
                                                                                            gTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTacca                                                                                                                                     >  1:423401/1‑236 (MQ=255)
                                                                                             ttCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCCCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGGCCGTCCCTGTTAATTacca                                                                                                                                     >  1:614324/1‑235 (MQ=255)
                                                                                                  gcttacGGGAAATGCGCCGTCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGGGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACcga                                                                                                                                 <  2:579613/230‑1 (MQ=255)
                                                                                                        gggAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTg                                                                                                                                                                                                                     >  2:221883/1‑145 (MQ=255)
                                                                                                        gggAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTg                                                                                                                                                                                                                     <  1:221883/145‑1 (MQ=255)
                                                                                                               gcgcGAGTAAGTGGTGGCGACGGCGAATCTGATGCGTGAGCGCTGCGTTAAACCGGGGGACCGCGAGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCcg                                                                                                                    <  2:360646/235‑1 (MQ=255)
                                                                                                               gcgcGAGCAGGTGGTGGCGCTGGCGAATCTGGTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCcg                                                                                                                    <  1:712265/235‑1 (MQ=255)
                                                                                                               gcccaaGCAGGTGGTAGCCATGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCcg                                                                                                                    <  2:515538/230‑1 (MQ=255)
                                                                                                                cgcgAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCCGCCATCCGGACGATCGCCTTGAAATGATTCTGGAAGCAGCGCGGCCCGCGCTGCTACTTTCCACCGACGATCAAATGcccag                                                                                                                  >  2:71184/1‑233 (MQ=255)
                                                                                                                            gtggCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCcg                                                                                                                    <  1:549580/222‑1 (MQ=255)
                                                                                                                             tggCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGCACAGCGTGGCGGTGGCACTGCCGCGATCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCg                                                                                                          <  2:284644/231‑1 (MQ=255)
                                                                                                                             aggCGCTGGGGAATCTGCTGCTGGAGCCGGGCGTTAACCCAGGGGACAGCGTGGCGGTGGCCCTACCGCGCTCGGTCTTTTTGACCCTGGCCCTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCg                                                                                                          <  2:193443/230‑1 (MQ=255)
                                                                                                                              ggCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCCATACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGtt                                                                                                     <  1:579150/235‑1 (MQ=255)
                                                                                                                              ggCACTGGCGAATCTGCTGCGTTAGCGGGGCGTTAAACCAGGGGACAGCGTGGCGGGGGCACTACCGCGTTCGGTCTTTTGGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGtt                                                                                                     <  2:614324/235‑1 (MQ=255)
                                                                                                                                   tgggaaaTCTGCTGCGAGCGCACGGCGTTAAACCAGGGGACTGAGTGGCGTTGGCTCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCaa                                                                                                <  2:199645/230‑1 (MQ=255)
                                                                                                                                   tGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgt                                                                                                                                                    >  1:320715/1‑183 (MQ=255)
                                                                                                                                   tGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGcgtccgt                                                                                                                                                    <  2:320715/183‑1 (MQ=255)
                                                                                                                                       taagctgctgCGTGAGCGTGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTaa                                                                                           <  1:753536/232‑1 (MQ=255)
                                                                                                                                                           gcgTTAAACCAGGGGACAGCGTGGCGGTGGCACCACCGCGCTCGGGCTTTTTTACCCCGGCACTCCCTGCGATAGTTGCAGCTGGGGCGGCCTGGCTCCCGCCGGGTACCGGCTATCCGGACGAACGCCCGAAAATGATGCTGGAAGATGCGCGTCCGCCGCCGTTAAATTCCACCGACTATCAACAGCCGCGCTTTTGCGCGGTTCCCAATTTATAAAGCCTCTGCtctacatc                                                                        >  2:474770/1‑231 (MQ=255)
                                                                                                                                                                  gccAGGGGACAGCGTGGCGGTGGCACTTCCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTa                                                                <  1:250939/235‑1 (MQ=255)
                                                                                                                                                                  acgaGGGGAAAGCATAGCGGTGGCCCTACCGCGCTCGTTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTa                                                                <  2:488204/233‑1 (MQ=255)
                                                                                                                                                                   ccAGGGGACAGGGTGGGGGTGGCACTACCGCGGTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTa                                                                <  1:199529/235‑1 (MQ=255)
                                                                                                                                                                            aGCGTGGCGGTGGCACTACCCCGCTCGGTCTTTTTTACCCTGGCACTCCCTGCGGTAGTTGCCGCTGGAGCGGCCTGGCTACCGCTGGCTACCGGCCATCCGGACGAACGCCTGAAAATGATGCTTGAAGATTCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCCTTTGCGATGTTTCCCATTTAACAAGCCCTTGGTATAACACCCCGCGTTCACCGcctgg                                                      >  2:69659/1‑232 (MQ=255)
                                                                                                                                                                             gCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaa                                                                                                                                                                              >  2:638409/1‑115 (MQ=255)
                                                                                                                                                                             gCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGaa                                                                                                                                                                              <  1:638409/115‑1 (MQ=255)
                                                                                                                                                                              cGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGa                                                                                                                                                                  <  2:738249/126‑1 (MQ=255)
                                                                                                                                                                              cGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGa                                                                                                                                                                  >  1:738249/1‑126 (MQ=255)
                                                                                                                                                                                        gCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCCTGCGATAGTTGCCGCTGGAGCGGCCTGGCTACCCCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGGAGATGCGCGTCCGTCGCTGTTCATTACCACCGACGATCAACTGCCGCGCTTTAGCGATTTGCCCCATTTTACAAGCCCTTGCTCTCACGGCCCGCTTACACCTCCGGGCAGTGcgccgc                                           >  2:494812/1‑235 (MQ=255)
                                                                                                                                                                                        gCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGc                                                                        <  1:648147/206‑1 (MQ=255)
                                                                                                                                                                                        gCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGc                                                                        >  2:648147/1‑206 (MQ=255)
                                                                                                                                                                                            tACCGCGTTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCa                                       <  1:71184/235‑1 (MQ=255)
                                                                                                                                                                                            tACCGCGTGCGGTCTGTGTGACCCTGGCACACAGTGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCa                                       <  2:36785/235‑1 (MQ=255)
                                                                                                                                                                                              ccGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAAc                                                                          <  2:21659/198‑1 (MQ=255)
                                                                                                                                                                                              ccGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAAc                                                                          >  1:21659/1‑198 (MQ=255)
                                                                                                                                                                                                   cTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTca                                >  1:279252/1‑235 (MQ=255)
                                                                                                                                                                                                      ggTCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAAc                            >  1:112365/1‑236 (MQ=255)
                                                                                                                                                                                                        tCTTTTTGACCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTAACAAcc                           >  1:627544/1‑235 (MQ=255)
                                                                                                                                                                                                               gACACTGGCACTCCAGGCGATGGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATcaca                   <  2:391633/236‑1 (MQ=255)
                                                                                                                                                                                                                aCCCTGGCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCacac                  <  1:197500/236‑1 (MQ=255)
                                                                                                                                                                                                                      gCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTTCGCCGCTGCAACTTTCACAACCGCATCACACGGCtt             >  1:257257/1‑235 (MQ=255)
                                                                                                                                                                                                                      gCACTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCtt             >  2:157182/1‑235 (MQ=255)
                                                                                                                                                                                                                        aCTCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTtata          <  1:496669/236‑1 (MQ=255)
                                                                                                                                                                                                                          tCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat         <  1:489708/235‑1 (MQ=255)
                                                                                                                                                                                                                          tCCATGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat         <  1:69659/235‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATcac                    <  2:550039/220‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATcac                    >  1:550039/1‑220 (MQ=255)
                                                                                                                                                                                                                               gCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCa                                                         >  2:719108/1‑182 (MQ=255)
                                                                                                                                                                                                                               gCGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCa                                                         <  1:719108/182‑1 (MQ=255)
                                                                                                                                                                                                                                cGATAGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCttt  <  1:474770/236‑1 (MQ=255)
                                                                                                                                                                                                                                    aGTTGCAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCCCACGGCTTATATCATCttt  >  1:263978/1‑232 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
GTTGAGATTCCAGAAACCACGCTTAGCGCGCTGGTGGCAGAACAAGCGGCAAAAACACCGGATGCTCCGGCGCTGGCAGATGCGCGTTACCTGTTCAGCTATCGGGAAATGCGCGAGCAGGTGGTGGCGCTGGCGAATCTGCTGCGTGAGCGCGGCGTTAAACCAGGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTT  >  NC_000913/615504‑615962

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