Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,903 G→A E523K (GAA→AAA)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,9030GA98.5% 225.9 / NA 67E523K (GAA→AAA) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/44);  minor base C (1/0);  total (23/44)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTA  >  NC_000913/1660‑2140
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
ggggTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACt                                                                                                                                                                                                                                        <  2:376010/249‑1 (MQ=255)
gTGGTAAATAACGATGATGCGACCACTGTCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACt                                                                                                                                                                                                                                        <  2:374658/251‑1 (MQ=255)
   gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGgc                                                                                                                                                                                                                                     <  1:437575/251‑1 (MQ=255)
         aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGc                                                                                                                                                                                                                               <  1:22915/251‑1 (MQ=255)
          aCGATGATGCGACCACTGGCGTGCGCGTTACGCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAAGGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGcc                                                                                                                                                                                                                              <  1:272787/251‑1 (MQ=255)
            gatgatGCTACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATGCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCGGGAAAACTGGCAGAAAGAACTGGCGCAAGCCaa                                                                                                                                                                                                                            <  2:92595/251‑1 (MQ=255)
                atgCGGCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAagag                                                                                                                                                                                                                        <  2:413692/251‑1 (MQ=255)
                atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAagag                                                                                                                                                                                                                        <  1:249757/251‑1 (MQ=255)
                         cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCCAGCCaa                                                                                                                                                                                                                            >  2:249473/1‑238 (MQ=255)
                         cTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCaa                                                                                                                                                                                                                            <  1:249473/238‑1 (MQ=255)
                            gagtgggcgTTACTCATCAGATGCTGTTCAATACCGATCAGGTGATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTGTAAtctc                                                                                                                                                                                                            <  2:7281/245‑1 (MQ=255)
                              gTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGCAGGGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAACCTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                                                                                                                          >  2:23447/1‑251 (MQ=255)
                                  gcgTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGgcgc                                                                                                                                                                                                      >  1:508581/1‑251 (MQ=255)
                                    gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTGCATGGCCTTAATCTGGAAATCTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCtt                                                                                                                                                                                                    <  1:378614/251‑1 (MQ=255)
                                       acacatcaGATGCTGTTCAATACTGATCAGGTTATCGAAGGGTGTGTGATTGGCGTCGGTGGCGTTGGCTGTGTGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCTCCAATGTACATGGCCTTAATCTGGAAAACTGTCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCTCTTAAt                                                                                                                                                                                                 <  2:385090/248‑1 (MQ=255)
                                           atcagctgcagctcactACCGATCCGGTTGTAGAAGTGTTTGTGATTGGCGCCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCCGCAAAGCTGGCTGCAGACTAAACCTATCGCCTTACGTGTCTGCGGCGTTGCCAACTCGCAGGCTCTGCTCACCACTGTAAATGGCCTTAATCTGGCAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGAGCTTAATTCGc                                                                                                                                                                                             >  2:576660/17‑251 (MQ=255)
                                            tcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTCATTCGcc                                                                                                                                                                                            >  2:206615/1‑251 (MQ=255)
                                             caGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACGGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCt                                                                                                                                                                                           <  2:290127/251‑1 (MQ=255)
                                               gATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCg                                                                                                                                                                                         <  1:352289/251‑1 (MQ=255)
                                               gATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCg                                                                                                                                                                                         >  1:519538/1‑251 (MQ=255)
                                               gATGCTGTTCAATACCGAGCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCg                                                                                                                                                                                         <  2:243938/251‑1 (MQ=255)
                                                        cattACCGAGCAGGTTATCGAAGGGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGTAAAGCTGGCTGAAGAATAATCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGACAACTCGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAat                                                                                                                                                                                <  1:481563/248‑1 (MQ=255)
                                                            aCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCGCGGGCGCTTAATTCGCCTCGTGAAAGAATatca                                                                                                                                                                            <  1:461494/251‑1 (MQ=255)
                                                            aCCGATCAGGTGATCGAAGTGTTTGTGATGGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATatc                                                                                                                                                                             <  1:533874/250‑1 (MQ=255)
                                                                aTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAAGCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTGTAATCTCGGGCTCTTAATTCGCCTCGTGAAAGAATATCATctg                                                                                                                                                                        <  1:531466/251‑1 (MQ=255)
                                                                 tccggTTATCGAAGTGATTGTAATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTGGAGCAACTGACGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGGGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCACGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgc                                                                                                                                                                       >  2:89976/4‑251 (MQ=255)
                                                                            aaGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTGTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTc                                                                                                                                                            <  2:199109/251‑1 (MQ=255)
                                                                                      gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt                                                                                                                                                  >  1:144720/1‑251 (MQ=255)
                                                                                       aTTGGCGTCGGCGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAGGCTGGCTGAAGAATAAAAATATCGACTTACGTGTCTGCGGCATTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGAGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTg                                                                                                                                                 >  1:66021/1‑251 (MQ=255)
                                                                                          ggcggcGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAAATAACATAGCGACTTACGTGACTGCGGTGTTGCCAACTCCAAGGCGCTGCTCACCAATGTACATGGCCGTAATCTGGAAAGCTGGCAGATAGAGCTGGTTCAAGCGAAAGAGTCGTGTAATCTCGGGCGCGTAATTCGCCGCGTTTAAGAATAACAGCTGCTGATCCCTGTCATGGTTGACTGCAc                                                                                                                                              <  2:311696/251‑1 (MQ=255)
                                                                                           gCGTCGGTGGCGTTGGCGGTGTGCTGCTGGAGCAACTGAAGCTTCAGCAAAGCTGGCTGAAGAATAAACAGATCGGCTTGCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACAGGGCCTTAATCTGGAAAACTGGCATAAAGAACTGGCGCAAGCGAAAGGGCCGTTTAATCTCGGGCGCTTAAGTCTCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt                                                                                                                                             <  1:42361/251‑1 (MQ=255)
                                                                                                   ggCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATGGTTGACTGCACTTccagcca                                                                                                                                     <  1:224828/251‑1 (MQ=255)
                                                                                                               gcgcTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCgg                                                                                                                          <  1:396708/250‑1 (MQ=255)
                                                                                                                        gAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCTGATCAATATGc                                                                                                                <  2:392362/251‑1 (MQ=255)
                                                                                                                                    ggTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcg                                                                                                    <  2:244714/250‑1 (MQ=255)
                                                                                                                                     gTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcgc                                                                                                   >  1:264366/1‑251 (MQ=255)
                                                                                                                                                   ctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACAGGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGa                                                                                                 <  1:323710/239‑1 (MQ=255)
                                                                                                                                                    tgAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                                                                                <  1:342938/239‑1 (MQ=255)
                                                                                                                                                    tgAAGAATAAACATATCGACTTACGCGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACAAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGACTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                                                                                >  2:342938/1‑239 (MQ=255)
                                                                                                                                                     gaagaaTAAACATTTCGACTTACGTGGCTGCGGTGTTGCCAACTCGCAGGCTCTGCACACCAATGTACATGGCCGTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATAGCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGTCTTCCTGCGCGAAGGTTTCCTCGTTg                                                                                   <  2:318028/251‑1 (MQ=255)
                                                                                                                                                          aTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACg                                                                              <  1:301210/251‑1 (MQ=255)
                                                                                                                                                                atatCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAAc                                                                        >  2:152365/1‑251 (MQ=255)
                                                                                                                                                                       cTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAACACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCCTGCGCGGAGGTTTCCACGTTGTCACGCCGACCAAAAAgg                                                                 >  2:540205/1‑251 (MQ=255)
                                                                                                                                                                           cGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGAAGTGGCGGCTCAAAATGCCGACTTCCTGAGCGAAGGTTTCCACGGTGTCCCGCCGACCAACAAGGCCa                                                              >  2:194543/1‑250 (MQ=255)
                                                                                                                                                                            gtggcTGCGGTGTTGCCAACTCGAAGGCTCTGCGCACCCATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCGCGCCGAACAAAAAGGCCAac                                                            <  2:161993/247‑1 (MQ=255)
                                                                                                                                                                                  gCGGTGTTGACAACTCGAAGGCTCTGCTCACCAATGTACATGGCCGTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGTGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTGCCAGCCAGTCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCTAACAAAAAGGCCAACACCtc                                                       <  2:353841/250‑1 (MQ=255)
                                                                                                                                                                                       gtTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCGCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACGCGGTCATTGTTGACTGCACTTCCAGCCAGGGAGTGGCGGATCAATATGCCGACGTCCTGGGCGAAGGTTTCCACGTTGTCACGCCGTACAAAAATGCCAGCAGCTCGTCGAt                                                 <  2:284983/251‑1 (MQ=255)
                                                                                                                                                                                        ttGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATg                                                >  1:480386/1‑251 (MQ=255)
                                                                                                                                                                                                 cGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGGAGGTTTCCACGTTGTAACGCCGAACAAAAAGGCCCACACCTCGTCGATGGATtactac                                       >  2:452471/1‑251 (MQ=255)
                                                                                                                                                                                                     ggctcttatCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGAGCAAGCCAACGAGTCGTGTAATCTCGGGCGCTTAGTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTACCTGCGCGAATGTTTCCACGTTGTCACGCCGATCAAAAAGGCCAACACCGCGTCGATGGATTCCTATCAta                                   <  2:121484/243‑2 (MQ=255)
                                                                                                                                                                                                       ctctGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATc                                   >  1:342232/1‑249 (MQ=255)
                                                                                                                                                                                                           gCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCCGGCAGTGGCGGACAAATATGCCGACTTCCTGAGCGACGGCTTCAACGTTATCACTCCGGACAAAAAGGCCAACACCACGCCGATGGATTACTACCATCAGTTGc                             >  1:550127/1‑251 (MQ=255)
                                                                                                                                                                                                            cTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGCTTCCACGTTGTCACGCCGAACAAAAAGGCCCACCCCTCGTCGATGGATTACTACCATCAGTTGCg                            >  2:142326/1‑251 (MQ=255)
                                                                                                                                                                                                              cACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACCGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATGGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACGTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGTAGTACTACCATCAGTTGca                            <  2:387715/249‑2 (MQ=255)
                                                                                                                                                                                                                 cAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCTTTATg                       >  2:66817/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaTGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATgc                      <  1:142326/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaTGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATgc                      <  2:206744/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAACAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCTTCAGTTGCGTTATgcg                     <  2:395089/251‑1 (MQ=255)
                                                                                                                                                                                                                     gTACATGGCCTTAATCGGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATTGTTGACGGCACTTCCTGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCTAAGGTTTCCGCGTTGTCACGCCGATCAAAAAGGCCAACACCGCGTCGATGGATTACGACTAGCAGGTGCGTTATgcggc                   <  2:534504/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGcggcgg                 >  1:81725/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggCCTTACTCTGGAAAACTGGCAGCAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGaaa              >  2:235050/1‑250 (MQ=255)
                                                                                                                                                                                                                            gccctactcttgAAAACTTGCAGAAAGAACAGGCGCAAGCCAAAGAGCCGTTTAAACTCGGGCGCTTATATCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTGTCCGCTTTTTCGCTCCGAACAAAAAGGCCAACGCATCGGCGATGGATTACTACCAGCAGTTGCGTTATGCGGCGCaaaac            <  2:550127/240‑2 (MQ=255)
                                                                                                                                                                                                                             ccTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAAt           <  1:540205/251‑1 (MQ=255)
                                                                                                                                                                                                                             ccTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATAGCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAAt           <  2:414425/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTCCCCTCAGTTGCGTTATGCGGCGGAAAAATCgcggcg    >  1:146380/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggAAAACTGGCAGAAAGAACTTGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCTCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTa  <  2:178327/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTGCCTTCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCAGCAGTTGCGTTATGCGGCGGAAAAATCgcgg      <  1:577993/247‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             
GTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTA  >  NC_000913/1660‑2140

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