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

ACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGT  >  NC_000913/2816950‑2817440
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
acagccACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGcaacg                                                                                                                                                                                                                                                   <  1:400695/247‑1 (MQ=255)
     cACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCACCGACGGCGGCGTCAAAATAATGGCGTAATCCCAGATGCGCCCGCAATGCCTCAGCGATGGAGCTTTCAC‑TCCCCGTTCCTCCCGCCATTGGGTAACGACCAt                                                                                                                                                                                                                                              >  2:117716/1‑251 (MQ=255)
        ggCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGcc                                                                                                                                                                                                                                           >  1:401875/1‑251 (MQ=255)
          cGTCCAGGCCTGTCGCACGTGCCGCCTTAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTGGCACGCCCATACGCTGCGCGCACAACAAAAATGTGGCTGGCGCGGGGTTATGGTGTTTTACGTGATCGGCAGCGTCGACGGCGTCTAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCGCAGCGAAGGCGATTTCAC‑GCCCCTTTCCGACAGCCGTTGGGCAACGACCATGCCaa                                                                                                                                                                                                                                         <  2:346020/251‑1 (MQ=255)
              cATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACTCACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGCCGCCGGCGCCAAAATATTGGCTTAATCCCAGGGGCCCCAGCCATGCCTCCGCGTGGGCGCTTTCAC‑TCCCCGTTCCTACACCAATTGGGCAACCCCCTTGCCAACttt                                                                                                                                                                                                                                     >  2:359539/1‑251 (MQ=255)
                 gCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTca                                                                                                                                                                                                                                  <  2:77927/251‑1 (MQ=255)
                  ccTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTcac                                                                                                                                                                                                                                 <  2:238412/251‑1 (MQ=255)
                         cACGGGCCGCCTGTCTACCGAAATCGGCATCTTCAAAGACCACACACAGCGTCGGTTGCACTCCCATCCGCTGCGCGCACAACAAAAATGTGTCCGGCGCGGGTTTATGGTGTTTGACGTCTTCGGCAGCGTCGACGGCCTCAAACTAATGTGGGCATCCCAGGGGCGCCAGAACTGCCTCAGTGCTGGATCTT‑CACATCCACGATCCTACAGCCATTGGCCAACGCCAATGCCATCTTTTCACCTCATc                                                                                                                                                                                                                           >  2:195913/1‑250 (MQ=255)
                             ggCCGCCTGAATTCCGAAATCGGCATCTTCAAAGACCACCCACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCTTAATCCCAGTTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATc‑aacaa                                                                                                                                                                                                                     <  1:241162/251‑1 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTCTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCAATTGGGCGACGACCATGCCAACTTTTCACCACATC‑AACCAGAGGAAGc                                                                                                                                                                                                             >  2:216410/1‑251 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCAGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGaa                                                                                                                                                                                                               <  1:7649/249‑1 (MQ=255)
                                     gAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAGGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAg                                                                                                                                                                                                              <  2:401875/250‑1 (MQ=255)
                                      aaTACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAA‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCg                                                                                                                                                                                                            >  1:237151/1‑251 (MQ=255)
                                         aCCGAAATCGGCATCTTCAAAGACCACACCCTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCATCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGtt                                                                                                                                                                                                         <  1:20339/251‑1 (MQ=255)
                                                 cGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCt                                                                                                                                                                                                  <  1:299911/250‑1 (MQ=255)
                                                       cctCAAAGACCACACACCTCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAAGTGTCTGGCGCGGGTTTATGGTTTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCcag                                                                                                                                                                                            <  1:30394/248‑1 (MQ=255)
                                                         tCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAAGGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACAGC‑AACAAGAGGAAGCGGTTCGACGCTATCcag                                                                                                                                                                                            <  2:62128/248‑1 (MQ=255)
                                                         tCAAAGACCACACACTGAGTCGGCTGCACGCCCATACGCTGCGCGCACACCAACAATGTGTCTCGCGCGGGATTATGGTGTTTGACGTGATCGCCAGCGACGACGGCGTCCAAATAATGGCGTAATCCCAGGCGCGCCAGCCATGCCTCAGCGATGGCGCTTTTAC‑TCCCCGCTCCACCAGCCCTTGGGACACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGAAGCTATCcagc                                                                                                                                                                                           >  2:245497/1‑249 (MQ=255)
                                                          cAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                    >  1:193137/1‑239 (MQ=255)
                                                          cAAAGACCACACACGGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACg                                                                                                                                                                                                    <  2:193137/239‑1 (MQ=255)
                                                             aGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATa                                                                                                                                                                                     <  1:476012/251‑1 (MQ=255)
                                                             aGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGAAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATa                                                                                                                                                                                     <  1:159099/251‑1 (MQ=255)
                                                              gACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGAGCCCAACAAAACTATGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGCAGGCGTCAAAATAATGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTCCC‑TCCCCAGTCCTACCGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAATAGGAAGCGGTTAGCAGCTATCCGCCAGCATAc                                                                                                                                                                                    >  2:473039/1‑251 (MQ=255)
                                                              gACCACACACTGCGTCGGTTGCAAGCCCATACCCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGCGTTTGACGTGAACGGCAGCGACGACGGCGTCAAAATAATGGCGGAATCCCCGGAGCGCCAGCAATGCCTAAGCGATGGCGCTTTCAC‑TCCACGTGCATACAGCCAATGGGCAACGACCATGCCAACTTTTCACCACATCCAACA‑GAGCAAGCGGTTCGGCGCAATCCAGCAGCATa                                                                                                                                                                                     >  2:28087/1‑250 (MQ=255)
                                                                     caccgcggcGGTTGCACGCCCATACGCTGCGCGCACAAAAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTAc                                                                                                                                                                             <  2:444120/243‑1 (MQ=255)
                                                                         ggttCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGc                                                                                                                                                                          <  2:438154/247‑1 (MQ=255)
                                                                          ggTCGGTTGCACGCCCATACGCTGCGAGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCtt                                                                                                                                                                        <  2:237300/250‑1 (MQ=255)
                                                                           gTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGTCGCTATCCAGCAGCATACTTCTTACTGCTTc                                                                                                                                                                       <  1:143029/251‑1 (MQ=255)
                                                                                 tGCACGCCCATACGCTGCGCGCACCACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATAGGGCAACGACCCTGCCACCTTTTACCCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttt                                                                                                                                                                 >  2:293729/1‑251 (MQ=255)
                                                                                                   gcgcACAACAAAAATGTGGCTGGCGCGGGTTTATGGGGTTTGACGGGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATg                                                                                                                                               <  2:362973/251‑1 (MQ=255)
                                                                                                    cgcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                                                                                                                                              >  1:155893/1‑251 (MQ=255)
                                                                                                     gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCg                                                                                                                                             <  1:211234/251‑1 (MQ=255)
                                                                                                        caacaaAGATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCTGCAGCGACTACGGCGTCAAAATAATGGCTTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCCC‑GCCCCGTTTCTACATCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTTTCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGt                                                                                                                                          <  1:206463/251‑1 (MQ=255)
                                                                                                           caaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCg                                                                                                                                        >  1:411728/1‑250 (MQ=255)
                                                                                                            aacaaTGGGTCTGGCGCGGGTTTATGGTGTTTGAAGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTATTCCCAGGTGCGCCAGCAATGCCTCAGCGAAGGAGCTTTCAC‑TCCCCGATCCTACAGCCATTGGGCAACGACCATTCCAACTTTTCACCACATC‑CACAAGAGGAAGGGGATCGAGGCTATCTAGCAGCATACTTCTTACTGCTTTGGTTTTTTCACGCGCTaccacctgcagagcg                                                                                                                                        >  2:392644/4‑235 (MQ=255)
                                                                                                                aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCg                                                                                                                                        >  1:180821/1‑245 (MQ=255)
                                                                                                                aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCg                                                                                                                                  <  1:285235/251‑1 (MQ=255)
                                                                                                                             gggTTTATTGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCt                                                                                                                     >  1:154195/1‑251 (MQ=255)
                                                                                                                             gggTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCt                                                                                                                     >  2:156213/1‑251 (MQ=255)
                                                                                                                                 ttATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGCTTCAGCTCaat                                                                                                                 >  2:500058/1‑251 (MQ=255)
                                                                                                                                         tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCt                                                                                                         >  1:422886/1‑251 (MQ=255)
                                                                                                                                           tGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTg                                                                                                        >  1:324727/1‑250 (MQ=255)
                                                                                                                                               gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑CCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCaa                                                                                                   <  2:175549/251‑1 (MQ=255)
                                                                                                                                                     ggCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGcc                                                                                             >  2:14876/1‑251 (MQ=255)
                                                                                                                                                        aGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAgg                                                                                          <  2:180821/251‑1 (MQ=255)
                                                                                                                                                        aGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAgg                                                                                          <  1:160370/251‑1 (MQ=255)
                                                                                                                                                               acgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCg                                                                                    >  1:188029/1‑250 (MQ=255)
                                                                                                                                                                   cGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATcc                                                                                >  1:110507/1‑250 (MQ=255)
                                                                                                                                                                    gTCAAAATAATGGCGTAATCCCAGGTGCGCCCGCAATGCCTCAGCGATGGCGCTTTCAC‑GCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACGTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCAt                                                                              <  1:448061/251‑1 (MQ=255)
                                                                                                                                                                        aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTa                                                                            >  1:352874/1‑249 (MQ=255)
                                                                                                                                                                         aataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCCGGTGGGCGATCCATTAAgc                                                                         >  2:381744/1‑251 (MQ=255)
                                                                                                                                                                                gCGTAAGCCCAGGTGCGCCAGCAATGCCTCAGCGTTGGCGCTTTCAC‑CCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAatca                                                                  <  2:110507/251‑1 (MQ=255)
                                                                                                                                                                                       cccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGCTCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTg                                                           >  2:158914/1‑251 (MQ=255)
                                                                                                                                                                                        ccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTg                                                           >  1:476342/1‑250 (MQ=255)
                                                                                                                                                                                                 ccAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGcagatcggaagagc                                                 <  1:26467/251‑14 (MQ=255)
                                                                                                                                                                                                  cAGCAATGCCCCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACt                                                <  2:352874/251‑1 (MQ=255)
                                                                                                                                                                                                         gCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATAGGCCTGATTCAGCTCAATAATTGCATGAGCAATACGCCAGGGGGGCGCTCCATTAAGCGCAATCATCGCCGGAATATCGTACGGAAGAcc                                         >  2:140006/1‑251 (MQ=255)
                                                                                                                                                                                                              aGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGt                                    <  1:291354/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGcc                                 >  2:122919/1‑251 (MQ=255)
                                                                                                                                                                                                                       gcTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATa                           <  1:87137/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttCAC‑CCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATAc                          <  2:496755/249‑1 (MQ=255)
                                                                                                                                                                                                                           ccAC‑TCCCCGTTCCTACAGCCAGTGGGCAACGACCATGCCAACTTTTCACCACACC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATATTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTc                       <  2:400170/250‑1 (MQ=255)
                                                                                                                                                                                                                            cAC‑TCCCCGTTCCTACAGCCATTGGGCACCGACCATGCCAACTTTTCACCACAGC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTcg                      <  2:411728/251‑1 (MQ=255)
                                                                                                                                                                                                                            cAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGCCCGTAGTGCCCTAATACTTcg                      >  2:260301/1‑251 (MQ=255)
                                                                                                                                                                                                                            cAC‑TCCCCGTTCCTACAGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTcg                      >  1:34993/1‑251 (MQ=255)
                                                                                                                                                                                                                                         caagaGCCTTTGGGCAACTAACATTCCAACTTTTCCCCACAGC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATTCTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGTCGATCCATGAAGCGCAATCAGCGCCTGAATATCTTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCTCGCTTTa          <  2:312617/247‑1 (MQ=255)
                                                                                                                                                                                                                                             aGCCATTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGt      <  2:168566/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTGGGCAACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGt  >  1:258297/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ACATCCACGGCGTCCATGCCTGCTGCACGGGCCGCCTGAATACCGAAATCGGCATCTTCAAAGACCACACACTGCGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCAC‑TCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATC‑AACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCCTGATTCAGCTCAATAATTGCCTGAGCAATACGCCAGGTGGGCGATCCATTAAGCGCAATCATCGCCTGAATATCGTACTGAAGACCGTAGTGCCCTAATACTTCGCGCCACGCTTTACGGTGCGT  >  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: ‑