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

AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGC  >  NC_000913/1801‑2271
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGttt                                                                                                                                                                                                                               <  1:373944/251‑1 (MQ=255)
          agaaTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATTAGGCGAATTCTTGCGCGAAGGTTTCCACGTTGTc                                                                                                                                                                                                                     >  1:476267/1‑251 (MQ=255)
               aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                                                                                                                                                                                                <  1:463057/251‑1 (MQ=255)
               aaaCATAACGACTTACGTGTCGGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                                                                                                                                                                                                <  2:163702/251‑1 (MQ=255)
                       cgattgaCGTGTCGTCGGTGGTGCCAACTAGAAGGAACTTCGAACCAATGTGCCTGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGTAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTTAGTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCTGAGCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGGCTCGCCGAACaaa                                                                                                                                                                                                         <  2:289873/244‑1 (MQ=255)
                        gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACcaaa                                                                                                                                                                                                        >  1:378867/1‑250 (MQ=255)
                         aCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACaaaa                                                                                                                                                                                                        >  1:329781/1‑249 (MQ=255)
                           ttACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGc                                                                                                                                                                                                    >  1:113085/1‑251 (MQ=255)
                             aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                                                                                                                                                                                                   <  1:110532/250‑1 (MQ=255)
                                    tGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCAACCCGCACAAAAAAGCCAACAccac                                                                                                                                                                                           >  2:599376/1‑249 (MQ=255)
                                          gtTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCGGGAAAACTGGCAGGAAGAACTGGCGAAAGCCAAAGTGCCGCTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCAGCTGCCGAACCCGGTCATTGTTGACTGCACTTCCAGCCATGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCCAACACCTCGTCGAt                                                                                                                                                                                     <  2:41477/251‑1 (MQ=255)
                                          gtTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGAt                                                                                                                                                                                     <  1:176529/251‑1 (MQ=255)
                                                    cGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATtactac                                                                                                                                                                           <  1:132325/251‑1 (MQ=255)
                                                              gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                  <  1:217660/250‑1 (MQ=255)
                                                              gCTCAACAATGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                  <  1:17525/250‑1 (MQ=255)
                                                                  aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGtc                                                                                                                                                              >  1:274721/1‑249 (MQ=255)
                                                                  aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTAGTCGATGGATTACTACCATCTGTTGCGtc                                                                                                                                                              >  1:125735/1‑249 (MQ=255)
                                                                    cAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATg                                                                                                                                                           <  1:599376/251‑1 (MQ=255)
                                                                                           aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                    <  1:497827/251‑1 (MQ=255)
                                                                                           aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGATGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaa                                                                                                                                     >  1:157664/1‑250 (MQ=255)
                                                                                                 ggCAGTAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCct                                                                                                                               <  1:489880/250‑1 (MQ=255)
                                                                                                  gCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTCACATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCctct                                                                                                                             >  1:222564/1‑251 (MQ=255)
                                                                                                    aGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAt                                                                                                                           <  1:477947/251‑1 (MQ=255)
                                                                                                           aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGACAGCATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGAGGCTCAATATGCCGGCTTCTTGCGCGAAGGCTTCCACGTTGTAACGCCGACCAAAAAGGCCCACACCTCGTCGCTGGATTACTACCATATGTTGCGTCTTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAccc                                                                                                                   >  1:616665/1‑250 (MQ=255)
                                                                                                           aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCa                                                                                                                   <  1:667552/251‑1 (MQ=255)
                                                                                                            ctagcgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACCAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGGAAAATCGCGGCGTAAATTTCTCTATGACACCa                                                                                                                   >  1:363332/4‑250 (MQ=255)
                                                                                                            cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTCAATTCCTCTATGACACCa                                                                                                                   >  1:363351/1‑250 (MQ=255)
                                                                                                                cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGtt                                                                                                              <  1:386768/251‑1 (MQ=255)
                                                                                                                  cAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTgg                                                                                                            <  2:274721/251‑1 (MQ=255)
                                                                                                                    aGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCTGTCAATGTTTACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTgggg                                                                                                          <  2:329781/251‑1 (MQ=255)
                                                                                                                         aatagctgTTTAATCTCGGGTGCTTAATTCACCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAGTGTTGACTGCACTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGCGAACGTTTCCACGTTGTCACGCCGAACAAAACGGCCAACACCTCGTCGATGTATTACTACCATCTGTTGCGTCATGCGGCGGATAAATCGCGGCGTAATTTCCTCTATGACTACAACGTTGGGGCTGGa                                                                                                     <  2:634879/244‑1 (MQ=255)
                                                                                                                          agtgCCGTTTAATCGCGGGCGCTTAACTCGCCTCGTGCAAGAATATCATCTGCTGAACCCTGTCTTTGTTGACTGCACGTCCAGCCAGGCAGTGGTGGATCAATATGCCGACTTCTTGAGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACTTCGGCGATCGATGACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACAAACGTTGGGGCTGGAt                                                                                                    <  2:378867/248‑1 (MQ=255)
                                                                                                                           gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCCACGTTGGGGCTGGAtt                                                                                                   >  1:566851/1‑251 (MQ=255)
                                                                                                                              ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTAcc                                                                                                <  1:384125/251‑1 (MQ=255)
                                                                                                                                gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCgg                                                                                              <  1:280962/251‑1 (MQ=255)
                                                                                                                                                     tCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAAt                                                                          <  1:576554/250‑1 (MQ=255)
                                                                                                                                                        ccTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGc                                                                      <  1:378470/251‑1 (MQ=255)
                                                                                                                                                              gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                    <  1:460972/184‑1 (MQ=255)
                                                                                                                                                              gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAAAATGCCGACTTCTTGCGCGAAGGTTTCCCCGTTGTCCCGCCGAAAAAAAAGGCCAACACCTCGTCGCTGGATTACTACCAACTGGTGCGTCATGCTGCGGAAAAATCGCGGCGTaaa                                                                                                                                    >  2:460972/1‑184 (MQ=255)
                                                                                                                                                                 aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGctg                                                             <  1:17787/251‑3 (MQ=255)
                                                                                                                                                                      atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAACTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                        <  1:391179/251‑1 (MQ=255)
                                                                                                                                                                       tcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCTGGCAGTGGCGGATCAATTTGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTTACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATTGCGGCGTTAATTATTCCATGCCACCAACGTTGGGGATTGCTTACCGGTTATTATGAACCTGCAAAAACTTCTCAATGCTGGtgt                                                         >  2:37647/1‑248 (MQ=255)
                                                                                                                                                                             gctgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAACAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGAGAATTCCTCTATTCCCACCACGTTGGGGCTGGATTACCGGTTATTTAGAACCTGCAAAATCTGCGCaaagatgggggagacttgg                                                 >  2:196533/1‑234 (MQ=255)
                                                                                                                                                                                     cGGTCACTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATAAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCATCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGtgc                                         <  2:58181/251‑3 (MQ=255)
                                                                                                                                                                                                     cACTTCCAGCCAGGCAGTGGCTGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCCCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                          <  1:540709/250‑2 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGGCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTGGCGTCATGCGGCGTAAAAATCGCGGCGTAAATTCCTCTATGTCACCAACGTTGGGGCTGGAGTACCGGTTATTGAGaat                                                                                   <  2:628942/192‑2 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGCACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTCTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                                                                   >  1:628942/1‑192 (MQ=255)
                                                                                                                                                                                                      aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGAGGCGTTAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTTAGGACCTGCAAAATCTGCCCAATGCTGGTGATGAATTGATGATGTTCTCCCGCATTCTTTTTgg                        >  2:264255/1‑251 (MQ=255)
                                                                                                                                                                                                         tCCAGCCAGTCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCTAGGTGGTCACGCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCCCGGCGTAAATTCCTCTATGACTCCAACGTTGGGGCTGGATTTCCTGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATTAATTTCGCCGGcgtggtttcaggtt                      <  2:112668/250‑14 (MQ=255)
                                                                                                                                                                                                              ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                <  1:673986/251‑1 (MQ=255)
                                                                                                                                                                                                              ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                <  1:563405/251‑1 (MQ=255)
                                                                                                                                                                                                                  gtaGTGGCGGATCAATCTGCCGTCCTCCTGCGCGAATGTTTCCATGTTGTCTCGTCGAACTAAAAGGCCAACACCTCCTCGATGGATTACTACCAACTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATCCCTCTGTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTta            <  2:240711/249‑1 (MQ=255)
                                                                                                                                                                                                                        gCGGATCAATAAGCCGACTTCTTGGCAGAAGGTTTCCACGTTGTCACTCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCtt      <  2:344015/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCg    <  1:163745/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc  >  1:459564/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGAGGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc  <  1:291906/251‑1 (MQ=255)
                                                                                                                                                                                                                              cAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATACGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAAt                                                                 >  1:96955/1‑186 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
AGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGC  >  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: