Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,817,197 G→A R103C (CGC→TGC)  yqaB ← fructose‑1‑phosphate phosphatase YqaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,817,1970GA98.5% 200.1 / ‑6.0 67R103C (CGC→TGC) yqaBfructose‑1‑phosphate phosphatase YqaB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGT  >  NC_000913/2816950‑2817440
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
acagccACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGcaacg                                                                                                                                                                                                                                                   <  1:400695/247‑1 (MQ=255)
     cACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCACCGACGGCGGCGTCAAAATAATGGCGTAATCCCAGATGCGCCCGCAATGCCTCAGCGATGGAGCTTTCACTCCCCGTTCCTCCCGCCATTGGGTAACGACCAt                                                                                                                                                                                                                                              >  2:117716/1‑251 (MQ=255)
        ggCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGcc                                                                                                                                                                                                                                           >  1:401875/1‑251 (MQ=255)
          cGTCCAGGCCTGTCGCACGTGCCGCCTTAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTGGCACGCCCATACGCTGCGCGCACAACAAAAATGTGGCTGGCGCGGGGTTATGGTGTTTTACGTGATCGGCAGCGTCGACGGCGTCTAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCGCAGCGAAGGCGATTTCACGCCCCTTTCCGACAGCCGTTGGGCAACGACCATGCCaa                                                                                                                                                                                                                                         <  2:346020/251‑1 (MQ=255)
              cATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACTCACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGCCGCCGGCGCCAAAATATTGGCTTAATCCCAGGGGCCCCAGCCATGCCTCCGCGTGGGCGCTTTCACTCCCCGTTCCTACACCAATTGGGCAACCCCCTTGCCAACttt                                                                                                                                                                                                                                     >  2:359539/1‑251 (MQ=255)
                 gCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTca                                                                                                                                                                                                                                  <  2:77927/251‑1 (MQ=255)
                  ccTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTcac                                                                                                                                                                                                                                 <  2:238412/251‑1 (MQ=255)
                         cACGGGCCGCCTGTCTACCGAAATCGGCATCTTCAAAGACCACACACAGCGTCGGTTGCACTCCCATCCGCTGCGCGCACAACAAAAATGTGTCCGGCGCGGGTTTATGGTGTTTGACGTCTTCGGCAGCGTCGACGGCCTCAAACTAATGTGGGCATCCCAGGGGCGCCAGAACTGCCTCAGTGCTGGATCTTCACATCCACGATCCTACAGCCATTGGCCAACGCCAATGCCATCTTTTCACCTCATc                                                                                                                                                                                                                           >  2:195913/1‑250 (MQ=255)
                             ggCCGCCTGAATTCCGAAATCGGCATCTTCAAAGACCACCCACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCTTAATCCCAGTTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATcaacaa                                                                                                                                                                                                                     <  1:241162/251‑1 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTCTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCAATTGGGCGACGACCATGCCAACTTTTCACCACATCAACCAGAGGAAGc                                                                                                                                                                                                             >  2:216410/1‑251 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCAGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGaa                                                                                                                                                                                                               <  1:7649/249‑1 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAGGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAg                                                                                                                                                                                                              <  2:401875/250‑1 (MQ=255)
                                      aaTACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAATCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCg                                                                                                                                                                                                            >  1:237151/1‑251 (MQ=255)
                                         aCCGAAATCGGCATCTTCAAAGACCACACCCTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCATCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                                                                                                                                                         <  1:20339/251‑1 (MQ=255)
                                                 cGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCt                                                                                                                                                                                                  <  1:299911/250‑1 (MQ=255)
                                                       cctCAAAGACCACACACCTCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAAGTGTCTGGCGCGGGTTTATGGTTTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcag                                                                                                                                                                                            <  1:30394/248‑1 (MQ=255)
                                                         tCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAGGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACAGCAACAAGAGGAAGCGGTTCGACGCTATCcag                                                                                                                                                                                            <  2:62128/248‑1 (MQ=255)
                                                         tCAAAGACCACACACTGAGTCGGCTGCACGCCCATACGCTGCGCGCACACCAACAATGTGTCTCGCGCGGGATTATGGTGTTTGACGTGATCGCCAGCGACGACGGCGTCCAAATAATGGCGTAATCCCAGGCGCGCCAGCCATGCCTCAGCGATGGCGCTTTTACTCCCCGCTCCACCAGCCCTTGGGACACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAAGCTATCcagc                                                                                                                                                                                           >  2:245497/1‑249 (MQ=255)
                                                          cAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                    >  1:193137/1‑239 (MQ=255)
                                                          cAAAGACCACACACGGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                    <  2:193137/239‑1 (MQ=255)
                                                             aGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATa                                                                                                                                                                                     <  1:476012/251‑1 (MQ=255)
                                                             aGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGAAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATa                                                                                                                                                                                     <  1:159099/251‑1 (MQ=255)
                                                              gACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGAGCCCAACAAAACTATGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGCAGGCGTCAAAATAATGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTCCCTCCCCAGTCCTACCGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAATAGGAAGCGGTTAGCAGCTATCCGCCAGCATAc                                                                                                                                                                                    >  2:473039/1‑251 (MQ=255)
                                                              gACCACACACTGCGTCGGTTGCAAGCCCATACCCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGCGTTTGACGTGAACGGCAGCGACGACGGCGTCAAAATAATGGCGGAATCCCCGGAGCGCCAGCAATGCCTAAGCGATGGCGCTTTCACTCCACGTGCATACAGCCAATGGGCAACGACCATGCCAACTTTTCACCACATCCAACAGAGCAAGCGGTTCGGCGCAATCCAGCAGCATa                                                                                                                                                                                     >  2:28087/1‑250 (MQ=255)
                                                                     caccgcggcGGTTGCACGCCCATACGCTGCGCGCACAAAAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTAc                                                                                                                                                                             <  2:444120/243‑1 (MQ=255)
                                                                         ggttCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGc                                                                                                                                                                          <  2:438154/247‑1 (MQ=255)
                                                                          ggTCGGTTGCACGCCCATACGCTGCGAGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCtt                                                                                                                                                                        <  2:237300/250‑1 (MQ=255)
                                                                           gTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGTCGCTATCCAGCAGCATACTTCTTACTGCTTc                                                                                                                                                                       <  1:143029/251‑1 (MQ=255)
                                                                                 tGCACGCCCATACGCTGCGCGCACCACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATAGGGCAACGACCCTGCCACCTTTTACCCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttt                                                                                                                                                                 >  2:293729/1‑251 (MQ=255)
                                                                                                   gcgcACAACAAAAATGTGGCTGGCGCGGGTTTATGGGGTTTGACGGGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATg                                                                                                                                               <  2:362973/251‑1 (MQ=255)
                                                                                                    cgcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                                                                                                                                              >  1:155893/1‑251 (MQ=255)
                                                                                                     gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCg                                                                                                                                             <  1:211234/251‑1 (MQ=255)
                                                                                                        caacaaAGATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCTGCAGCGACTACGGCGTCAAAATAATGGCTTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCCCGCCCCGTTTCTACATCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTTTCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGt                                                                                                                                          <  1:206463/251‑1 (MQ=255)
                                                                                                           caaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCg                                                                                                                                        >  1:411728/1‑250 (MQ=255)
                                                                                                            aacaaTGGGTCTGGCGCGGGTTTATGGTGTTTGAAGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTATTCCCAGGTGCGCCAGCAATGCCTCAGCGAAGGAGCTTTCACTCCCCGATCCTACAGCCATTGGGCAACGACCATTCCAACTTTTCACCACATCCACAAGAGGAAGGGGATCGAGGCTATCTAGCAGCATACTTCTTACTGCTTTGGTTTTTTCACGCGCTaccacctgcagagcg                                                                                                                                        >  2:392644/4‑235 (MQ=255)
                                                                                                                aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCg                                                                                                                                        >  1:180821/1‑245 (MQ=255)
                                                                                                                aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCg                                                                                                                                  <  1:285235/251‑1 (MQ=255)
                                                                                                                             gggTTTATTGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCt                                                                                                                     >  1:154195/1‑251 (MQ=255)
                                                                                                                             gggTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCt                                                                                                                     >  2:156213/1‑251 (MQ=255)
                                                                                                                                 ttATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGCTTCAGCTCaat                                                                                                                 >  2:500058/1‑251 (MQ=255)
                                                                                                                                         tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCt                                                                                                         >  1:422886/1‑251 (MQ=255)
                                                                                                                                           tGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                        >  1:324727/1‑250 (MQ=255)
                                                                                                                                               gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACCCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCaa                                                                                                   <  2:175549/251‑1 (MQ=255)
                                                                                                                                                     ggCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGcc                                                                                             >  2:14876/1‑251 (MQ=255)
                                                                                                                                                        aGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAgg                                                                                          <  2:180821/251‑1 (MQ=255)
                                                                                                                                                        aGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAgg                                                                                          <  1:160370/251‑1 (MQ=255)
                                                                                                                                                               acgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCg                                                                                    >  1:188029/1‑250 (MQ=255)
                                                                                                                                                                   cGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATcc                                                                                >  1:110507/1‑250 (MQ=255)
                                                                                                                                                                    gTCAAAATAATGGCGTAATCCCAGGTGCGCCCGCAATGCCTCAGCGATGGCGCTTTCACGCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCAt                                                                              <  1:448061/251‑1 (MQ=255)
                                                                                                                                                                        aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTa                                                                            >  1:352874/1‑249 (MQ=255)
                                                                                                                                                                         aataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCCGGTGGGCGATCCATTAAgc                                                                         >  2:381744/1‑251 (MQ=255)
                                                                                                                                                                                gCGTAAGCCCAGGTGCGCCAGCAATGCCTCAGCGTTGGCGCTTTCACCCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAatca                                                                  <  2:110507/251‑1 (MQ=255)
                                                                                                                                                                                       cccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGCTCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTg                                                           >  2:158914/1‑251 (MQ=255)
                                                                                                                                                                                        ccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTg                                                           >  1:476342/1‑250 (MQ=255)
                                                                                                                                                                                                 ccAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGcagatcggaagagc                                                 <  1:26467/251‑14 (MQ=255)
                                                                                                                                                                                                  cAGCAATGCCCCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACt                                                <  2:352874/251‑1 (MQ=255)
                                                                                                                                                                                                         gCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATAGGCCTGATTCAGCTCAATAATTGCATGAGCAATACGCCAGGGGGGCGCTCCATTAAGCGCAATCATCGCCGGAATATCGTACGGAAGAcc                                         >  2:140006/1‑251 (MQ=255)
                                                                                                                                                                                                              aGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGt                                    <  1:291354/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGcc                                 >  2:122919/1‑251 (MQ=255)
                                                                                                                                                                                                                       gcTTTCACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATa                           <  1:87137/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttCACCCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATAc                          <  2:496755/249‑1 (MQ=255)
                                                                                                                                                                                                                           ccACTCCCCGTTCCTACAGCCAGTGGGCAACGACCATGCCAACTTTTCACCACACCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATATTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTc                       <  2:400170/250‑1 (MQ=255)
                                                                                                                                                                                                                            cACTCCCCGTTCCTACAGCCATTGGGCACCGACCATGCCAACTTTTCACCACAGCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTcg                      <  2:411728/251‑1 (MQ=255)
                                                                                                                                                                                                                            cACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGCCCGTAGTGCCCTAATACTTcg                      >  2:260301/1‑251 (MQ=255)
                                                                                                                                                                                                                            cACTCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTcg                      >  1:34993/1‑251 (MQ=255)
                                                                                                                                                                                                                                         caagaGCCTTTGGGCAACTAACATTCCAACTTTTCCCCACAGCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATTCTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGTCGATCCATGAAGCGCAATCAGCGCCTGAATATCTTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCTCGCTTTa          <  2:312617/247‑1 (MQ=255)
                                                                                                                                                                                                                                             aGCCATTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGt      <  2:168566/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTGGGCAACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGt  >  1:258297/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGT  >  NC_000913/2816950‑2817440

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