Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 908,290 T→C 15.0% intergenic (‑8/+3) ybjT ← / ← ltaE putative NAD‑dependent oxidoreductase/L‑allo‑threonine aldolase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,2900TC15.0% 124.8 / 2.6 55intergenic (‑8/+3)ybjT/ltaEputative NAD‑dependent oxidoreductase/L‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (30/14);  new base C (3/6);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 6.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGC  >  NC_000913/908066‑908509
                                                                                                                                                                                                                                |                                                                                                                                                                                                                            
gAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGGGTGCGCACCATATGCTGACCACTGTAGCCACTGGCCCCGGGAACTAAAAAGCGTTGCCGGACGTCTCTCCCCt                                                                                                                                                                                                                           >  1:88852/1‑228 (MQ=255)
   attagcacgttacGATTTCCGTCAACAGGCCCGGTAGGTTATCCGGCCAGCTGAGACTGTGTTTATGGCAATGGACGGTTGCCAGTTGACATTTGGCTAGCAGGGTGCCATGACGGGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCGGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTaa                                                                                                                                                                                                                        <  1:501780/216‑1 (MQ=255)
        cccGTATCGATATCCTTCAAAAGGGCGGGTTGGTTACCCGGCCAGTTGAGATCGCTTTTATGGCAACTGACGTTTGGCAGTTGAAGCTTTGCAAGCCTGTCGCATTGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACgcgc                                                                                                                                                                                                                   <  1:381211/227‑1 (MQ=255)
        aCCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGCTGCCCTTGCTGGCTGAGTGTGCGCACCAGCTGCTGACCAATGTAGCCACTGGCACCGAGAACTACAATCCATTGCGGCCCGTTTCCCTCCCTAAcccc                                                                                                                                                                                                                    >  1:531814/1‑224 (MQ=255)
        aCCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATCCTGACCAATGTAGCCACTGGCCCCGAGAACTAAAATCCGTTGCGGCCCGTCTCTCTCCTTAAcgcg                                                                                                                                                                                                                    >  2:652725/1‑227 (MQ=255)
             aTCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGAGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTTTCCTTAACGCCCCAGGa                                                                                                                                                                                                              >  2:368255/1‑228 (MQ=255)
                  tttccggaagcAGGGCCGGTGGGTTGTGCGGCCATGGCGGATCGACTTTAGGCCAACTGCCGTTGGCCAGTTGCATCTTTGCAGGCGTGTCGAACTGACGTCCCGCTGCCAGTATCTGATGCCCGTGCTGGCTGAGTGGGCGCCCCAGATGCTGACCAATGTAGCCACTGGCCCCGAGATCTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCa                                                                                                                                                                                                         <  1:190830/218‑1 (MQ=255)
                    tCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTACTGCAACTGACGTTTGCCAGTTGCAGCCTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGCCGCCCTTGCTGGCTGAGTGTGCGCCCCAGATCCTGCCCCATTTACCCCCCGGCCCCGCGACCTAAAATCAGTTGCGGCCCCCCTTTCTCCCTacaccccccgacatgccc                                                                                                                                                                                                        >  1:171098/1‑210 (MQ=255)
                    tCCTGAAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCAAGTTGCGGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCa                                                                                                                                                                                                         <  2:52781/226‑1 (MQ=255)
                       tGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATCCTGACCAATGTAGCCACTGGCACCGGGAACTCAAATCCGTTGCGGCCCGTCTCTCTCCTTAACCCCCCAGGCAAGCCCGca                                                                                                                                                                                                     >  1:655531/1‑226 (MQ=255)
                       tGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGGGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCCCCTTAACGCGCCAGGAATGCACGcc                                                                                                                                                                                                     >  2:473543/1‑227 (MQ=255)
                       tGCAACAGGGCCGGTAGGTTAACCGGCCAGCTGAGATCGTCTTTATGGCACCTGACGTTTGCCAGTTTCAGTTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCa                                                                                                                                                                                                    <  2:627253/228‑1 (MQ=255)
                                  cgttcggTTACCCGGCCGGCTGAGTTCGCCTTTATGGCACCTGACGTTTGCCAGTTGCAGCTTTGCAAGCTTGTCGACATGACGTGCCGCCGCCAGGAGCTGAGGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGa                                                                                                                                                                                         <  2:419326/223‑1 (MQ=255)
                                      aGGTTACCCGGCCAGCTGAGATCGACTTTATGGCAACGGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTggctacattggtc                                                                                                                                                                                     <  2:88852/228‑12 (MQ=255)
                                           accgggCCAGCTGGATTCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGcc                                                                                                                                                                                 <  2:531814/223‑1 (MQ=255)
                                               ggCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGGGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAACGCACGCCAGTGGCCGGCCACTTCCGccgctt                                                                                                                                                                             >  2:223582/1‑223 (MQ=255)
                                                 ccAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAAGGCGTTGCGCCACGTCCCTCTCCTTAACGCGCCAGGAATGCACCCCAGTGGGCGGCCACTTCCCCCCgttgt                                                                                                                                                                           >  2:448745/1‑227 (MQ=255)
                                                                         aCTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCACCAGGAATGCACGCCAGTGGGCGGCGTCTTCCGCCAGTTGTTCGCGCGAGACGTAACGATGCGTCa                                                                                                                                                  >  2:345879/1‑228 (MQ=255)
                                                                         aCTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACCGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCAAGTGGGCGGCGACTTCCGCCAGTCGTTCGCCCGAGAGGTCAAGATCCGCCa                                                                                                                                                  >  2:654061/1‑228 (MQ=255)
                                                                          cTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGCAGCCACCGGCACCGAGACCTAAAATGCGTTGCGGCCCGCCTCTCTCCTTAACGCCCCAGGCTTGCACGCCAGGGGGCGGGCAATCCCCCCGCTTGCTCGCGCGAGACGTccaactgcgcca                                                                                                                                                  >  1:601446/1‑216 (MQ=255)
                                                                                              cTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCGGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCCCTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAgg                                                                                                                             <  2:520404/228‑1 (MQ=255)
                                                                                                                aTGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGTTGCGTCACCAGGCGGCCAATCGCCCAGCCGTTAATCAGCACGTTtc                                                                                                           >  2:608825/1‑228 (MQ=255)
                                                                                                                        ccgccgccAGGACTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGAATCGTACCCAGGCGGACAATCGCCCAGGCGTTAATAGGCACGTTTCTCGCtttac                                                                                                  >  2:594534/1‑226 (MQ=255)
                                                                                                                               ctggctCTGGTTCCCTTGCGGGCTGAGTTGGCGCCCCAGATGTTGACCAATTTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTAt                                                                                             <  1:431787/222‑1 (MQ=255)
                                                                                                                               cAGGATCAGATTCCCTGGCGGGCTGGGTGTACGGAGCAGATGCTGACCAATGTAGCCACTGGCACCGAAAACTAAAATGCGTTGCGGCACGTCTCTCGCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTAt                                                                                             <  2:171098/227‑1 (MQ=255)
                                                                                                                                        aTGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCCAGGCGTTAATCAGCACGTTTTTCGCTTTCATGGATTCGCCTAAca                                                                                   >  2:487986/1‑227 (MQ=255)
                                                                                                                                           cccTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCCGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGGCACCAGGCGGACAATCGGCGTGGCTCTAATCAGATCGTTTTCCCCTTTCAGGGATTCGCCTAAagcca                                                                                >  1:644989/1‑226 (MQ=255)
                                                                                                                                             cTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCCCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGTGGGGTTAATCCGCACGTTTCTCGCGTTCATGTAATCGCCTAACGCGgc                                                                               >  2:193511/1‑227 (MQ=255)
                                                                                                                                               tGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAACATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCACGTTTTTCTCGCGCGACGGCCAGATGCGTCCCCAGGCGGACAATCGGCGAGCAGTTACGCACACCGATTCTCGTCTTCAGGTATTAACCTaaaggcgaag                                                                             >  1:386790/1‑219 (MQ=255)
                                                                                                                                               tGCTGGCTGAGTGTGCGCACCAGATGCCGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGgcagc                                                                           <  2:655531/228‑1 (MQ=255)
                                                                                                                                                              cgcACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGAGTTCCGCCAGTTGTTCGCGCGAGACGTCAAAAAGCGTCGCCAGGCGGACAATCGCCGAAGCGTTAATCAGCACGTTTCTCGCTTGCAGGTATTCGCCTAACCCGGCAGGCATTTTTTCCCCCAc                                                            >  1:361318/1‑228 (MQ=255)
                                                                                                                                                              cgcACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCCTTTTCTTCCCCGa                                                             >  2:32921/1‑227 (MQ=255)
                                                                                                                                                               gcccaaGATGCTGACCAATGTAGCCACTGGCAACGAGAACTAAAATGCGTTGCGGCACGTCTTTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGAcg                                                           <  2:146046/223‑1 (MQ=255)
                                                                                                                                                                        gCTGACCAATTTAGCCCCTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCGCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGaacaaa                                                   <  2:517811/227‑1 (MQ=255)
                                                                                                                                                                           gACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAg                                                >  2:520821/1‑227 (MQ=255)
                                                                                                                                                                           gACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTCCTCGCTTTCATGTATTCGCCTAACGCGGCACCATTTCCTTCCCCCACCCCAACAAACAg                                                >  1:432503/1‑227 (MQ=255)
                                                                                                                                                                            cccatttttgcccgtGGCACCGAGGACAAAAATGCGTTGCGGCACGTCCCTCCCTTTACCGCGCCAGGAATGCCCGCCAGTGGGCGGCGACTCCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCa                                              <  1:520821/214‑1 (MQ=255)
                                                                                                                                                                                      ccACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTGAATCAGCACGTTTCTCGATTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAACCAGCATATTgggtt                                     >  2:340384/1‑224 (MQ=255)
                                                                                                                                                                                        aCTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCAAGTTTCTCGCTTTCATGAATCTCCCTACGCCGGCAGCATTTCCTTCCCCCCCCCGCACAAACCGCATATTGGTGTcc                                   >  1:275374/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGTCAGCATTTTCTTCCCCGCCGCGCACACAACGCATATTTGTGTCCTGa                                >  2:414249/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGa                                >  2:104732/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGAAGCATTTTCTTCCCCGCCGCGAACAAACAGCATATTGGGGTCCCGa                                >  1:513626/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGAATTCGCCTAACGCGGCCGCATTTTCTTCCCCGCCGCGCACAAACCCCAAATTggggcccgtac                               >  1:540941/1‑219 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCCCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGCCGCGAACAACCAGCATATTGGTGTCCTgc                                >  2:70287/1‑226 (MQ=255)
                                                                                                                                                                                               ggcacaaCTAAAATTGTTTGGGGCCCGTCTCTCCCCCTAACGCGCCAGGAATGCACGCCAGTGGGGGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCAt                           <  1:440299/223‑1 (MQ=255)
                                                                                                                                                                                               acGAGAACTAAGATGCGTTGCGGCACGTCTCCCTCCTTAACGGGCCAGGAATGGCCGCCGGTCGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCCGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCAt                           <  2:532257/227‑1 (MQ=255)
                                                                                                                                                                                                         aaaTGGCTTGCGGCACGGCGCTCCCCTTAACGCGCCAGGAATGCCCGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCCCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCcgc                  <  2:275374/227‑1 (MQ=255)
                                                                                                                                                                                                         aaaTGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCCTGTATTCCCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATGGGTGTCCTGCCGCCTCACATcccc                  >  1:370380/1‑225 (MQ=255)
                                                                                                                                                                                                          attggtttGCGGGCCGTCTCCCCCTGTAACGGGCCCGGAACGCACCCCGTGGGGCGGCGACTTCCGCCAGGTGGTCGCGCGAGGCGGCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTGCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCgcgc                <  1:193511/222‑1 (MQ=255)
                                                                                                                                                                                                          acggccttGCGGCCGGTCTCGCCCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCgcgc                <  2:540941/222‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCCCCTgtcc          >  2:440299/1‑224 (MQ=255)
                                                                                                                                                                                                                 tGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGCCGTTAATCAGCACGTTTCTCGCttt                                                                                                    >  2:222964/1‑138 (MQ=255)
                                                                                                                                                                                                                     gcccGTCTCCTTCCTTAACGCGCCAGGAATGCACGACAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGTCGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGc     <  2:644989/225‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTCTCTCCCCTAACCGCGCCAGGAAGGCCCGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAg   <  2:361318/228‑1 (MQ=255)
                                                                                                                                                                                                                        cGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGGCAATCGGCTAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCTGACAAACAGCAGATTGGCGCCCTGACGCCTCACCCCCCGCCCTGCTTAACGccta  >  1:654570/1‑225 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                            
GAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCAGCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGC  >  NC_000913/908066‑908509

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