Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,778 C→G 6.4% S62S (TCC→TCG) ‡ ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7780CG 6.4% 433.5 / 14.1 134S62S (TCC→TCG) ‡ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base C (59/65);  new base G (6/3);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 137 total aligned reads displayed.

GAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGG  >  NC_000913/2863629‑2863927
                                                                                                                                                     |                                                                                                                                                     
gaggagATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTcc                                                                                                                                                       <  2:1512036/150‑1 (MQ=255)
  ggagATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaa                                                                                                                                                     >  2:692636/1‑150 (MQ=255)
    agaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAg                                                                                                                                                   <  1:997703/150‑1 (MQ=255)
    agaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAg                                                                                                                                                   <  2:339674/150‑1 (MQ=255)
    agaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTAGGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAg                                                                                                                                                   <  1:962067/150‑1 (MQ=255)
      aTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTc                                                                                                                                                 <  2:1741814/150‑1 (MQ=255)
        gACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                               >  1:929399/1‑150 (MQ=255)
         aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgcc                                                                                                                                              >  1:450991/1‑150 (MQ=255)
           gcgGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTGCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgccgc                                                                                                                                            <  2:1097876/150‑1 (MQ=255)
            cgGATTGCCAGTTCATTCTTTCACCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCCACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCtgcg                                                                                                                                           <  2:1154605/150‑1 (MQ=255)
             ggATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCgg                                                                                                                                          >  1:1377745/1‑150 (MQ=255)
             ggATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCgg                                                                                                                                          >  1:1140518/1‑150 (MQ=255)
                 tGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATgg                                                                                                                                      <  1:1006138/150‑1 (MQ=255)
                  gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                     <  1:480619/150‑1 (MQ=255)
                  gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                     >  2:1620033/1‑150 (MQ=255)
                  gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                     >  1:526844/1‑150 (MQ=255)
                  gCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGGCGCGGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                     <  1:1125378/150‑1 (MQ=255)
                   ccAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCg                                                                                                                                    <  2:873028/150‑1 (MQ=255)
                   ccAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCGGACGGGAATTTACTGGCGACCCCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCg                                                                                                                                    >  2:510106/1‑150 (MQ=255)
                          aTTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGtt                                                                                                                             >  1:860045/1‑150 (MQ=255)
                             cTTTCAGCGCGGCTATGCAACCGGTTCGGCAGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCTCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                                                                          <  1:1316443/150‑1 (MQ=255)
                              tttCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgt                                                                                                                         <  1:1370415/150‑1 (MQ=255)
                                 cAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGCCAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtggt                                                                                                                      >  2:270198/1‑150 (MQ=255)
                                 cAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtggt                                                                                                                      <  2:1162414/150‑1 (MQ=255)
                                  aGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtggtg                                                                                                                     >  2:1314497/1‑150 (MQ=255)
                                  aGCGCGGCTATGCAACCGGTTAGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACTGGTTCATGCCTGGTCAATCTCGATCCGCAGCGTCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtggtg                                                                                                                     <  2:546629/150‑1 (MQ=255)
                                    cgcgGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGAc                                                                                                                   <  2:1702866/150‑1 (MQ=255)
                                    cgcgGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGAc                                                                                                                   >  2:950909/1‑150 (MQ=255)
                                        gCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAc                                                                                                               <  1:29152/150‑1 (MQ=255)
                                        gCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAc                                                                                                               <  2:921398/150‑1 (MQ=255)
                                         cTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAcc                                                                                                              >  1:1905897/1‑150 (MQ=255)
                                            tGCAACCGGTGCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTc                                                                                                           <  2:655514/150‑1 (MQ=255)
                                                   ggTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGACTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                    >  2:474259/1‑150 (MQ=255)
                                                    gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTTCTGGCGACACCGACAGGTTCATGCATGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAgg                                                                                                   <  2:926887/150‑1 (MQ=255)
                                                        ggCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGct                                                                                               <  1:1510344/150‑1 (MQ=255)
                                                        ggCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGct                                                                                               <  2:1912410/150‑1 (MQ=255)
                                                           tGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCTGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCtt                                                                                            >  1:293007/1‑150 (MQ=255)
                                                            ggCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCttt                                                                                           >  1:798019/1‑150 (MQ=255)
                                                            ggCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCttt                                                                                           >  1:1075004/1‑150 (MQ=255)
                                                             gCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTc                                                                                          >  1:690766/1‑150 (MQ=255)
                                                              cAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                                                                         <  1:1597686/150‑1 (MQ=255)
                                                              cAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                                                                         <  2:779241/150‑1 (MQ=255)
                                                                aTCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCCTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGCACCTCTTTcgag                                                                                       <  1:1105465/150‑4 (MQ=255)
                                                                 tCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGCACCTCTTTcgagg                                                                                      <  1:476929/150‑5 (MQ=255)
                                                                   tGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                                                                    >  1:886704/1‑150 (MQ=255)
                                                                      ggctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgcgc                                                                                 <  2:1843953/149‑1 (MQ=255)
                                                                        ctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTg                                                                               <  2:1377745/150‑1 (MQ=255)
                                                                        ctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGAAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTg                                                                               <  1:86130/150‑1 (MQ=255)
                                                                           ctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTAt                                                                            >  1:1080412/1‑150 (MQ=255)
                                                                            ttttACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGCAAGAGGTGCTCTTTCATCTGGCGCTGTATc                                                                           >  2:235065/1‑150 (MQ=255)
                                                                             tttACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCg                                                                          <  1:983375/150‑1 (MQ=255)
                                                                                 ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                                                                      <  2:1909115/150‑1 (MQ=255)
                                                                                 ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                                                                      >  1:795475/1‑150 (MQ=255)
                                                                                    gACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCGAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCCTCTGGCGCTGTATCGCCACAAt                                                                   >  1:123142/1‑150 (MQ=255)
                                                                                    gACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAAt                                                                   <  1:824259/150‑1 (MQ=255)
                                                                                    gACGAGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAAt                                                                   >  2:646309/1‑150 (MQ=255)
                                                                                      cGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATcc                                                                 >  1:37419/1‑150 (MQ=255)
                                                                                         gAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgcg                                                              <  1:1470799/150‑1 (MQ=255)
                                                                                           aTTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCTCAGCGGCTTTCCAAATTCGCCGCGGATGGCGAATGGTTATGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCt                                                            <  1:714630/150‑1 (MQ=255)
                                                                                           aTTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAATGTCGCCGCGGATGGCGTATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTTCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCt                                                            <  1:1324460/150‑1 (MQ=255)
                                                                                            tttACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTg                                                           <  1:1900264/150‑1 (MQ=255)
                                                                                               aCTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTaa                                                        >  1:541167/1‑150 (MQ=255)
                                                                                                  ggCGACACCGAGAGGTTCATGCCTGGGCAATCTCGATCCTCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc                                                     <  2:293007/150‑1 (MQ=255)
                                                                                                  ggCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc                                                     >  2:469440/1‑150 (MQ=255)
                                                                                                    cGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCgg                                                   >  2:1445944/1‑150 (MQ=255)
                                                                                                     gACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggt                                                  <  1:751275/150‑1 (MQ=255)
                                                                                                      acacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGGCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtg                                                 <  1:1808600/150‑1 (MQ=255)
                                                                                                       cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtgg                                                <  2:593893/150‑1 (MQ=255)
                                                                                                          cGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGc                                             >  2:1902853/1‑150 (MQ=255)
                                                                                                             cAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTATAGCGGTGGTGCAtt                                          >  1:1687952/1‑150 (MQ=255)
                                                                                                              aGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAttt                                         <  1:1256710/150‑1 (MQ=255)
                                                                                                               ggTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCAGCGCTGTAAAGCGGTGGTGCATTTg                                        >  1:951702/1‑150 (MQ=255)
                                                                                                                 ttCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                                      <  1:1470141/150‑1 (MQ=255)
                                                                                                                 ttCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                                      <  1:692636/150‑1 (MQ=255)
                                                                                                                   cATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcaca                                    >  1:1180843/1‑150 (MQ=255)
                                                                                                                       ccTGGGCAATCTCCATCAGCAGCGGCTGTCCAAAGTCGCCGCGGATGGCGAATGTTTAAGTGGGGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCTCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagcac                                <  1:1006693/150‑1 (MQ=255)
                                                                                                                       catGGGCAATCTCGATCCGCAGAGGCTTTCCAAACGCGCCGCTGATGGCGAAGGGTTCAGGGGGGACAAACCCTCGAAAGAGGTGCGCTTTCATCTGGCGCGGNATCGCAACTCTCCGCGCTGTAAAGCGGTGGTGCATTTgcacagcac                                <  1:1000241/148‑1 (MQ=255)
                                                                                                                         tGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTGgcccggcactg                              >  1:264377/1‑148 (MQ=255)
                                                                                                                          gggCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATg                             >  2:1233346/1‑150 (MQ=255)
                                                                                                                          gggCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATg                             >  2:777188/1‑150 (MQ=255)
                                                                                                                            gCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGt                           <  1:976860/150‑1 (MQ=255)
                                                                                                                               aTCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGa                        >  2:624634/1‑150 (MQ=255)
                                                                                                                                tctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGAc                       <  2:691660/150‑1 (MQ=255)
                                                                                                                                  gccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTGTCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACgg                     >  2:591039/4‑150 (MQ=255)
                                                                                                                                   ccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                    >  2:933425/3‑150 (MQ=255)
                                                                                                                                   ccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                    >  2:1371920/3‑150 (MQ=255)
                                                                                                                                   cGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                    >  2:1877708/1‑150 (MQ=255)
                                                                                                                                     aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCGACAATCCGCGCTGTAATGCGGTGGTGCATTTGCACAGCACATGGTCGACGgcgc                  >  1:277711/1‑150 (MQ=255)
                                                                                                                                     aTCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgcgc                  <  1:1314497/150‑1 (MQ=255)
                                                                                                                                         gCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTc              >  2:664228/1‑150 (MQ=255)
                                                                                                                                         gCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTc              <  1:1584131/150‑1 (MQ=255)
                                                                                                                                          cAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGCTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcc             <  1:844096/150‑1 (MQ=255)
                                                                                                                                           aGCGTCTTTCCAAAGTCGCCGCGGATGGCGAATTGTTAAGTGGTTACAAACCATCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGAGCTGTAAAGCGGTGGTGCATTTGCGCAGCAAATGGTCGACGTCCCTATcct            >  2:1901430/1‑150 (MQ=255)
                                                                                                                                           aGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcct            >  1:1716501/1‑150 (MQ=255)
                                                                                                                                             cGGCTTTCCAAAGTCGCCGCGGAAGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCCCACTATCGCAACAATACGCGCTTTAAAACGGTGGTGCAAATACACAGCACATGGTCGACCGCGCtatacttc          >  2:1724100/1‑143 (MQ=255)
                                                                                                                                                cTTTCCAAAGTCACCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctgcctg       >  1:1367907/1‑150 (MQ=255)
                                                                                                                                                 tttCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGGGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCctgcctgc      <  2:263114/150‑1 (MQ=255)
                                                                                                                                                  ttCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCa     <  2:798019/150‑1 (MQ=255)
                                                                                                                                                  ttCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCa     <  2:1075004/150‑1 (MQ=255)
                                                                                                                                                     cAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAgg  >  1:1726690/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGG  >  NC_000913/2863629‑2863927

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