Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,038 C→T L568L (CTG→TTG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,0380CT96.5% 165.0 / ‑4.5 57L568L (CTG→TTG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (20/35);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGC  >  NC_000913/1801‑2271
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGttt                                                                                                                                                                                                                               <  1:373944/251‑1 (MQ=255)
          agaaTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATTAGGCGAATTCTTGCGC‑GAAGGTTTCCACGTTGTc                                                                                                                                                                                                                     >  1:476267/1‑251 (MQ=255)
               aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGcc                                                                                                                                                                                                                <  1:463057/251‑1 (MQ=255)
               aaaCATAACGACTTACGTGTCGGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGcc                                                                                                                                                                                                                <  2:163702/251‑1 (MQ=255)
                       cgattgaCGTGTCGTCGGTGGTGCCAACTAGAAGGAACTTCGAACCAATGTGCCTGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGTAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTTAGTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCTGAGCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGGCTCGCCGAACaaa                                                                                                                                                                                                         <  2:289873/244‑1 (MQ=255)
                        gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGGTGTCACGCCGAACcaaa                                                                                                                                                                                                        >  1:378867/1‑250 (MQ=255)
                         aCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACaaaa                                                                                                                                                                                                        >  1:329781/1‑249 (MQ=255)
                           ttACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGc                                                                                                                                                                                                    >  1:113085/1‑251 (MQ=255)
                             aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                                                                                                                                                                                                   <  1:110532/250‑1 (MQ=255)
                                    tGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCAACCCGCACAAAAAAGCCAACAccac                                                                                                                                                                                           >  2:599376/1‑249 (MQ=255)
                                          gtTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCGGGAAAACTGGCAGGAAGAACTGGCGAAAGCCAAAGTGCCGCTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCAGCTGCCGAACCCGGTCATTGTTGACTGCACTTCCAGCCATGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCCAACACCTCGTCGAt                                                                                                                                                                                     <  2:41477/251‑1 (MQ=255)
                                          gtTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGAt                                                                                                                                                                                     <  1:176529/251‑1 (MQ=255)
                                                    cGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATtactac                                                                                                                                                                           <  1:132325/251‑1 (MQ=255)
                                                              gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                  <  1:217660/250‑1 (MQ=255)
                                                              gCTCAACAATGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                  <  1:17525/250‑1 (MQ=255)
                                                                  aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGtc                                                                                                                                                              >  1:274721/1‑249 (MQ=255)
                                                                  aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTAGTCGATGGATTACTACCATCTGTTGCGtc                                                                                                                                                              >  1:125735/1‑249 (MQ=255)
                                                                    cAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATg                                                                                                                                                           <  1:599376/251‑1 (MQ=255)
                                                                                           aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                    <  1:497827/251‑1 (MQ=255)
                                                                                           aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGATGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaa                                                                                                                                     >  1:157664/1‑250 (MQ=255)
                                                                                                 ggCAGTAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCct                                                                                                                               <  1:489880/250‑1 (MQ=255)
                                                                                                  gCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTCACATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCctct                                                                                                                             >  1:222564/1‑251 (MQ=255)
                                                                                                    aGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAt                                                                                                                           <  1:477947/251‑1 (MQ=255)
                                                                                                           aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGACAGCATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGAGGCTCAATATGCCGGCTTCTTGCGC‑GAAGGCTTCCACGTTGTAACGCCGACCAAAAAGGCCCACACCTCGTCGCTGGATTACTACCATATGTTGCGTCTTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACAccc                                                                                                                   >  1:616665/1‑250 (MQ=255)
                                                                                                           aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCa                                                                                                                   <  1:667552/251‑1 (MQ=255)
                                                                                                            ctagcgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACCAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGGAAAATCGCGGCGTAAATTTCTCTAT‑GACACCa                                                                                                                   >  1:363332/4‑250 (MQ=255)
                                                                                                            cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTCAATTCCTCTAT‑GACACCa                                                                                                                   >  1:363351/1‑250 (MQ=255)
                                                                                                                cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGtt                                                                                                              <  1:386768/251‑1 (MQ=255)
                                                                                                                  cAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTgg                                                                                                            <  2:274721/251‑1 (MQ=255)
                                                                                                                    aGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCTGTCAATGTTTACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTgggg                                                                                                          <  2:329781/251‑1 (MQ=255)
                                                                                                                         aatagctgTTTAATCTCGGGTGCTTAATTCACCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAGTGTTGACTGCACTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGC‑GAACGTTTCCACGTTGTCACGCCGAACAAAACGGCCAACACCTCGTCGATGTATTACTACCATCTGTTGCGTCATGCGGCGGATAAATCGCGGCGTAATTTCCTCTAT‑GACTACAACGTTGGGGCTGGa                                                                                                     <  2:634879/244‑1 (MQ=255)
                                                                                                                          agtgCCGTTTAATCGCGGGCGCTTAACTCGCCTCGTGCAAGAATATCATCTGCTGAACCCTGTCTTTGTTGACTGCACGTCCAGCCAGGCAGTGGTGGATCAATATGCCGACTTCTTGAGC‑GAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACTTCGGCGATCGATGACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACAAACGTTGGGGCTGGAt                                                                                                    <  2:378867/248‑1 (MQ=255)
                                                                                                                           gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCCACGTTGGGGCTGGAtt                                                                                                   >  1:566851/1‑251 (MQ=255)
                                                                                                                              ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTAcc                                                                                                <  1:384125/251‑1 (MQ=255)
                                                                                                                                gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCgg                                                                                              <  1:280962/251‑1 (MQ=255)
                                                                                                                                                     tCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAAt                                                                          <  1:576554/250‑1 (MQ=255)
                                                                                                                                                        ccTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGc                                                                      <  1:378470/251‑1 (MQ=255)
                                                                                                                                                              gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                    <  1:460972/184‑1 (MQ=255)
                                                                                                                                                              gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAAAATGCCGACTTCTTGCGC‑GAAGGTTTCCCCGTTGTCCCGCCGAAAAAAAAGGCCAACACCTCGTCGCTGGATTACTACCAACTGGTGCGTCATGCTGCGGAAAAATCGCGGCGTaaa                                                                                                                                    >  2:460972/1‑184 (MQ=255)
                                                                                                                                                                 aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGctg                                                             <  1:17787/251‑3 (MQ=255)
                                                                                                                                                                      atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAACTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                        <  1:391179/251‑1 (MQ=255)
                                                                                                                                                                       tcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCTGGCAGTGGCGGATCAATTTGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTTACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATTGCGGCGTTAATTATTCCAT‑GCCACCAACGTTGGGGATTGCTTACCGGTTATTATGAACCTGCAAAAACTTCTCAATGCTGGtgt                                                         >  2:37647/1‑248 (MQ=255)
                                                                                                                                                                             gctgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACAAACAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGAGAATTCCTCTATTCCCACCA‑CGTTGGGGCTGGATTACCGGTTATTTAGAACCTGCAAAATCTGCGCaaagatgggggagacttgg                                                 >  2:196533/1‑234 (MQ=255)
                                                                                                                                                                                     cGGTCACTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATAAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCATCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGtgc                                         <  2:58181/251‑3 (MQ=255)
                                                                                                                                                                                                     cACTTCCAGCCAGGCAGTGGCTGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCCCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                          <  1:540709/250‑2 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGGCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTGGCGTCATGCGGCGTAAAAATCGCGGCGTAAATTCCTCTAT‑GTCACCAACGTTGGGGCTGGAGTACCGGTTATTGAGaat                                                                                   <  2:628942/192‑2 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGCACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTCT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                                                                   >  1:628942/1‑192 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGAGGCGTTAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTTAGGACCTGCAAAATCTGCCCAATGCTGGTGATGAATTGATGATGTTCTCCCGCATTCTTTTTgg                        >  2:264255/1‑251 (MQ=255)
                                                                                                                                                                                                         tCCAGCCAGTCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCTAGGTGGTCACGCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCCCGGCGTAAATTCCTCTAT‑GACTCCAACGTTGGGGCTGGATTTCCTGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATTAATTTCGCCGGcgtggtttcaggtt                      <  2:112668/250‑14 (MQ=255)
                                                                                                                                                                                                              ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                <  1:673986/251‑1 (MQ=255)
                                                                                                                                                                                                              ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                <  1:563405/251‑1 (MQ=255)
                                                                                                                                                                                                                  gtaGTGGCGGATCAATCTGCCGTCCTCCTGCGC‑GAATGTTTCCATGTTGTCTCGTCGAACTAAAAGGCCAACACCTCCTCGATGGATTACTACCAACTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATCCCTCTGT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTta            <  2:240711/249‑1 (MQ=255)
                                                                                                                                                                                                                        gCGGATCAATAAGCCGACTTCTTG‑GCAGAAGGTTTCCACGTTGTCACTCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCtt      <  2:344015/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCg    <  1:163745/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc  >  1:459564/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGAGGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc  <  1:291906/251‑1 (MQ=255)
                                                                                                                                                                                                                              cAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATACGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAAt                                                                 >  1:96955/1‑186 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGC  >  NC_000913/1801‑2271

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