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

GAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGC  >  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)
                                                                                                                        ccgccgccAGGA‑CTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGAATCGTACCCAGGCGGACAATCGCCCAGGCGTTAATAGGCACGTTTCTCGCtttac                                                                                                  >  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)
                                                                                                                                               tGCTGGCTGAGTGTGCGCACCAGATGCCGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGgca‑gc                                                                           <  2:655531/228‑1 (MQ=255)
                                                                                                                                                              cgcACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGAGTTCCGCCAGTTGTTCGCGCGAGACGTCAAAAAGCGTCGCCAGGCGGACAATCGCCGAAGCGTTAATCAGCACGTTTCTCGCTTGCAGGTATTCGCCTAACCCGGCAGGCATTTT‑TTCCCCCAc                                                            >  1:361318/1‑228 (MQ=255)
                                                                                                                                                              cgcACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCCTTTTCTTCCCCGa                                                             >  2:32921/1‑227 (MQ=255)
                                                                                                                                                               gcccaaGATGCTGACCAATGTAGCCACTGGCAACGAGAACTAAAATGCGTTGCGGCACGTCTTTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGAcg                                                           <  2:146046/223‑1 (MQ=255)
                                                                                                                                                                        gCTGACCAATTTAGCCCCTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCGCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGaacaaa                                                   <  2:517811/227‑1 (MQ=255)
                                                                                                                                                                           gACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAg                                                >  2:520821/1‑227 (MQ=255)
                                                                                                                                                                           gACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTCCTCGCTTTCATGTATTCGCCTAACGCGGCA‑CCATTTCCTTCCCCCACCCCAACAAACAg                                                >  1:432503/1‑227 (MQ=255)
                                                                                                                                                                            cccatttttgcccgtGGCACCGAGGACAAAAATGCGTTGCGGCACGTCCCTCCCTTTACCGCGCCAGGAATGCCCGCCAGTGGGCGGCGACTCCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCa                                              <  1:520821/214‑1 (MQ=255)
                                                                                                                                                                                      ccACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTGAATCAGCACGTTTCTCGATTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAACCAGCATATTgggtt                                     >  2:340384/1‑224 (MQ=255)
                                                                                                                                                                                        aCTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCAAGTTTCTCGCTTTCATGAATCTCCCTACGCCGGCA‑GCATTTCCTTCCCCCCCCCGCACAAACCGCATATTGGTGTcc                                   >  1:275374/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGTCA‑GCATTTTCTTCCCCGCCGCGCACACAACGCATATTTGTGTCCTGa                                >  2:414249/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGa                                >  2:104732/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGAA‑GCATTTTCTTCCCCGCCGCGAACAAACAGCATATTGGGGTCCCGa                                >  1:513626/1‑227 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGAATTCGCCTAACGCGGCC‑GCATTTTCTTCCCCGCCGCGCACAAACCCCAAATTggggcccgtac                               >  1:540941/1‑219 (MQ=255)
                                                                                                                                                                                           ggCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCCCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGCCGCGAACAACCAGCATATTGGTGTCCTgc                                >  2:70287/1‑226 (MQ=255)
                                                                                                                                                                                               ggcacaaCTAAAATTGTTTGGGGCCCGTCTCTCCCCCTAACGCGCCAGGAATGCACGCCAGTGGGGGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCAt                           <  1:440299/223‑1 (MQ=255)
                                                                                                                                                                                               acGAGAACTAAGATGCGTTGCGGCACGTCTCCCTCCTTAACGGGCCAGGAATGGCCGCCGGTCGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCCGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCAt                           <  2:532257/227‑1 (MQ=255)
                                                                                                                                                                                                         aaaTGGCTTGCGGCACGGCGCTCCCCTTAACGCGCCAGGAATGCCCGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCCCCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCcgc                  <  2:275374/227‑1 (MQ=255)
                                                                                                                                                                                                         aaaTGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCCTGTATTCCCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATGGGTGTCCTGCCGCCTCACATcccc                  >  1:370380/1‑225 (MQ=255)
                                                                                                                                                                                                          attggtttGCGGGCCGTCTCCCCCTGTAACGGGCCCGGAACGCACCCCGTGGGGCGGCGACTTCCGCCAGGTGGTCGCGCGAGGCGGCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTGCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCgcgc                <  1:193511/222‑1 (MQ=255)
                                                                                                                                                                                                          acggccttGCGGCCGGTCTCGCCCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCgcgc                <  2:540941/222‑1 (MQ=255)
                                                                                                                                                                                                                 tGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCCCCTgtcc          >  2:440299/1‑224 (MQ=255)
                                                                                                                                                                                                                 tGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGCCGTTAATCAGCACGTTTCTCGCttt                                                                                                    >  2:222964/1‑138 (MQ=255)
                                                                                                                                                                                                                     gcccGTCTCCTTCCTTAACGCGCCAGGAATGCACGACAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGTCGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGc     <  2:644989/225‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTCTCTCCCCTAACCGCGCCAGGAAGGCCCGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAg   <  2:361318/228‑1 (MQ=255)
                                                                                                                                                                                                                        cGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGGCAATCGGCTAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCTGACAAACAGCAGATTGGCGCCCTGACGCCTCACCCCCCGCCCTGCTTAACGccta  >  1:654570/1‑225 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                            
GAAAATAGACCGTATCGATATCCTGCAACAGGGCCGGTAGGTTATCCGGCCAGCTGAGATCGACTTTATGGCAACTGACGTTTGCCAGTTGCAGCTTTGCAAGCCTGTCGACATGACGTGCCGCCGCCAGGATCTGATGCCCTTGCTGGCTGAGTGTGCGCACCAGATGCTGACCAATGTAGCCACTGGCACCGAGAACTAAAATGCGTTGCGGCACGTCTCTCTCCTTAACGCGCCAGGAATGCACGCCAGTGGGCGGCGACTTCCGCCAGTTGTTCGCGCGAGACGTCAAGATGCGTCACCAGGCGGACAATCGGCGAGGCGTTAATCAGCACGTTTCTCGCTTTCATGTATTCGCCTAACGCGGCA‑GCATTTTCTTCCCCGACGCGAACAAACAGCATATTGGTGTCCTGACGCATCACATCCGCGCCTGCTTCACGCAGC  >  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: ‑