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.4% 189.5 / NA 63E523K (GAA→AAA) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (29/31);  minor base C (1/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTC  >  NC_000913/1663‑2145
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gtaaaacaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCCGGTTATAGAAGTGTTTGTGATTGGAGTCGGTGGCGTTGGCGGTGCGCTGATGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATAGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGTTATGCTCACCACTGTAATAGGCCTTAATTTGGAACACTAGCagaaagaaacgg                                                                                                                                                                                                                                             >  2:490948/8‑246 (MQ=255)
gTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTgg                                                                                                                                                                                                                                             <  1:490948/250‑1 (MQ=255)
  aaaTAACGATTATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGgcg                                                                                                                                                                                                                                           <  2:223307/250‑1 (MQ=255)
               gCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGc                                                                                                                                                                                                                              <  1:511104/250‑1 (MQ=255)
                    cACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGttt                                                                                                                                                                                                                         >  1:19052/1‑250 (MQ=255)
                    cACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGAGGGGCGCTGCTGGAGCAACTGAAGCGTCAGAAAAGCTGGCTGAAGAATAAACATAACGACTTACGTGTCTGGGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTCACTGGCCTTACTCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAACGCCGttt                                                                                                                                                                                                                         >  2:288345/1‑250 (MQ=255)
                    cACTGGCGCGCGCGCTACTCATCAAATGCTGTTCAATACCGATCCGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTTGAGCAACTGAAGCGTCAGCCAAGCTGGCTGAAGAATAAACATATCGACTTACGTACCTGCGGTGTTGCCAACTCGAAGGCTATGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGGGCCGttt                                                                                                                                                                                                                         >  2:532947/1‑250 (MQ=255)
                        ggcgtgcgctgtactattctgtttgtGTTCAATACCGTACAGGCTATCGAAGTCTCTGTGAGTGGCGGCGGTGGCTATGGCGGGGCGCTGCTGGAGCAACTGATGCGTCAGCAATGCTGGCTGAAGAATTACCATATCGACTGACGTGTCCGCGTGGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGGCTTAATCTGGAAATCTGGCAGAAAGACCTGGCGCAAGCCAAAGTGCCGTTTAAtct                                                                                                                                                                                                                    <  2:275017/226‑1 (MQ=255)
                               gcgTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAGGAATAAACATATCGACTTACGTGTATGCGGTGTTGCCAACTCGAAGGCCCTGCTCACCAATGTACATGGCCTAAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTaa                                                                                                                                                                                                                       >  1:523634/1‑241 (MQ=255)
                               gcgTTACTCATCAGATGCTGTTCAATACCGATCAGGTGTTCGAAGTGTTTGTGATTGGCGTCGGTGGCGGTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCGGAAGAATAAACATATCGACTTACGTGTCTGCGGTTTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTTGAAAACTGTCAGAAAGAACTGGCGCAAGCCAAAGAGCCTTTTaa                                                                                                                                                                                                                       <  2:523634/241‑1 (MQ=255)
                                    aCTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAATAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCGCCAATGTTCATGGCCTTAATCTGGAAAGCTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTaa                                                                                                                                                                                                         <  1:334601/250‑1 (MQ=255)
                                       catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCATCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                                                      <  1:287593/250‑1 (MQ=255)
                                                cTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCATCGTGa                                                                                                                                                                                             >  2:437363/1‑250 (MQ=255)
                                                  gTTCAATACCGCCTAGGTTATTGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAACCGTCAGCAAAGCTGGCTGAAGAATAAACATAACGCCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTATGCTCACCAATGTACATGGCCTTAATATTGAAAACTGGCAGAAAGAACTGGCGCTAGCCAAAGAGCCGTTTAATCGCGGGCGCTTAATTCGCCTCGTGaaa                                                                                                                                                                                           <  2:299922/250‑1 (MQ=255)
                                                     cAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAat                                                                                                                                                                                       <  1:264212/251‑1 (MQ=255)
                                                            gATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATct                                                                                                                                                                                >  1:37400/1‑251 (MQ=255)
                                                               acccgccgcccgcccgatctgATGGGCGTCGGTGGCGTTGGCGGTGCTCAGCTTGAGTAACTGAAGCGTCAGTAAAGCTGGCTGAAGAATAAACATACCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTTCTCACCAATGTACAATGCCTTAATCGGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATGCGCCTCGTGAAAGAATATCATctgc                                                                                                                                                                              <  2:330665/231‑1 (MQ=255)
                                                                    tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAACGAGCCGTTTAATCTCGGGCGCTTAATTCCCCTCGTGACAGAATATCATCTGCTGAAcc                                                                                                                                                                        >  2:450764/1‑251 (MQ=255)
                                                                     aTCGAGGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACGTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                                                                                                                                        <  2:153870/250‑1 (MQ=255)
                                                                     aTCGAAGCGTTTGAGATTGGCGTCGGAGGCGTTGGCGGTGCGCTGCTGGAGAACCTGAAGCGTCAGCAAGGCTGGCTGAAGACTAAACATATCGACTTACGTGTCTTGGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGGAAACTGGCAGTAAGAACTGGCACAAGCCAAAGAGCCGTTTAATCCCGGGCGCTTAGTTCGCCTCgggtcatataatcctatgctaaaacc                                                                                                                                                                       >  2:286811/1‑226 (MQ=255)
                                                                                  tgATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGAc                                                                                                                                                          >  1:133037/1‑251 (MQ=255)
                                                                                  tgATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTAAATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGACTATCATCTGCTGAACCCGGTCATTGTTGAc                                                                                                                                                          >  2:468547/1‑251 (MQ=255)
                                                                                   gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgc                                                                                                                                                                              >  1:330665/1‑230 (MQ=255)
                                                                                    aTTGTCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGAGTTGCCAACTCGAAGGCTCTGCTCACCCATGTACATGGCCTTAATCTGGAAAACTGGCAGATAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTCATTCGCCTCGTGCAAGAATCTCATCTGCTGAACACGGTCATTGTTGACTg                                                                                                                                                        >  2:100745/1‑251 (MQ=255)
                                                                                         cGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACtt                                                                                                                                                   >  1:203964/1‑251 (MQ=255)
                                                                                             ggTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccag                                                                                                                                               <  1:571384/251‑1 (MQ=255)
                                                                                              gTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagc                                                                                                                                              >  1:211182/1‑251 (MQ=255)
                                                                                                ggCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACGGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagcca                                                                                                                                            <  1:237474/251‑1 (MQ=255)
                                                                                                 gCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTccagccag                                                                                                                                           <  1:321164/251‑1 (MQ=255)
                                                                                                      ggCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGt                                                                                                                                      <  1:258549/251‑1 (MQ=255)
                                                                                                         ggTGCGCTGCTGGAGCAACTTAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGc                                                                                                                                   <  1:428633/251‑1 (MQ=255)
                                                                                                                    ggAGCAACTGAAGCGTCATCAAAGCTGGCTGAAGAATAAACATATCGACTTTCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCCCTTCCAGCCAGGCAGTGGCGGATCAata                                                                                                                          <  1:339390/249‑1 (MQ=255)
                                                                                                                           cTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGCAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACtt                                                                                                                 >  1:69686/1‑251 (MQ=255)
                                                                                                                                 cGTCAGCAAAGCTGGCTGAGGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGAGGCGGagctcgggagggcggcgtgt                                                                                                             >  2:124536/1‑230 (MQ=255)
                                                                                                                                 cGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCgg                                                                                                                                 <  1:124536/229‑1 (MQ=255)
                                                                                                                                  gTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcc                                                                                                                                                 <  2:364523/212‑1 (MQ=255)
                                                                                                                                  gTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcc                                                                                                                                                 >  1:364523/1‑212 (MQ=255)
                                                                                                                                        aaaGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGt                                                                                                   <  1:536373/251‑1 (MQ=255)
                                                                                                                                                    agaaTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCACAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCAAGCCAGGCCGTGGCGGATCAATGTGCCGACTTCCTCCGGCGAGGTTTCCACGGTGTc                                                                                      >  2:76389/1‑251 (MQ=255)
                                                                                                                                                       aTAACCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGTCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAATGTTTCCACGtt                                                                                         <  2:275569/245‑1 (MQ=255)
                                                                                                                                                       aTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGtt                                                                                         >  1:275569/1‑245 (MQ=255)
                                                                                                                                                                  gatttACGTGCCTGCGGTGTTGCCAATTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGTCGCAAGTCAAAGAGCCGTTTAATCTCTGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTTCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa                                                                        <  1:274941/248‑1 (MQ=255)
                                                                                                                                                                  gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt                                                                                                                                                    <  2:536921/177‑1 (MQ=255)
                                                                                                                                                                  gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACt                                                                                                                                                    >  1:536921/1‑177 (MQ=255)
                                                                                                                                                                    cTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAgg                                                                      >  1:95228/1‑251 (MQ=255)
                                                                                                                                                                       aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                                                      <  2:81563/122‑1 (MQ=255)
                                                                                                                                                                       aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                                                      >  1:81563/1‑122 (MQ=255)
                                                                                                                                                                       aCGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCa                                                                   <  1:214223/251‑1 (MQ=255)
                                                                                                                                                                                 gatgtTGCCAACTCGAAGGCTCTGCTCACCAAAGTACATGGCCTTAATCTGGACAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGACGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCCATATGCCGACTTCCTGCGCGGAGGTATACCATTTGTCCCGCCGAACCAAAAAGCCAAAAccccacc                                                         >  2:325010/3‑246 (MQ=255)
                                                                                                                                                                                        ccAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAATGCCAACACCACGTCGATgagt                                                  >  2:180582/1‑248 (MQ=255)
                                                                                                                                                                                               gAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTAcc                                           <  1:96981/251‑1 (MQ=255)
                                                                                                                                                                                                         cTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAACGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGATTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCAATTCCAGCCAGGCAGTGGCGGATCAATATCCCCTCTTCCAGCGCGCAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCaaccggtgcg                                 >  1:318481/1‑242 (MQ=255)
                                                                                                                                                                                                         cTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCCCGTTGTCACGACGCACAAAAAAGCCCACCCCTCGTCGATGGGTTACTACCATCCGTTGCg                                 >  2:145568/1‑251 (MQ=255)
                                                                                                                                                                                                            aCCAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTc                                                                                               >  1:326129/1‑186 (MQ=255)
                                                                                                                                                                                                              cAATGTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATg                            <  2:203964/251‑1 (MQ=255)
                                                                                                                                                                                                                  gTACATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATgcgg                         <  1:468547/250‑1 (MQ=255)
                                                                                                                                                                                                                    aCATGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGcggcgg                      <  1:532947/251‑1 (MQ=255)
                                                                                                                                                                                                                        tGGCCTTAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGaa                    >  1:135607/1‑250 (MQ=255)
                                                                                                                                                                                                                               tAATCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTTCGTTATGCGGCGGaaaaaacac             >  2:448983/1‑247 (MQ=255)
                                                                                                                                                                                                                                 aTCTGGAAAACTGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCgcgg           >  1:493187/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aacgggCAGAAATAACTGGCCCAATCCAAAGTGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTTCTGTATCCGTTCATTTTTGACTGCACTTCCATCCAGGCATTGGCGGATCAATATTCCGACTTCCTGCGAGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAAtt   <  2:497326/247‑1 (MQ=255)
                                                                                                                                                                                                                                          ccgggCAGAAAGAACTGGCGCAAGCCAAAGAGGCCTGTAATCTCGGGCGCATAATTCTACTGGTGAAATAATAGCATCCGCTGAACCCGGTCATCGTTGACTGCACTTCCAGCCCGGAAGTGGCGGATCAATATGCCGCCGTCCTGCGCGAAGGGTTCCACGTTGTCACGCCGATCAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCTGACAAATCGCGGCGTAACTTc  <  2:148809/248‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                                                                                     >  1:453781/1‑166 (MQ=255)
                                                                                                                                                                                                                                            tGGCAGAAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTGCCAGCCAGGCAGTGGTGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                                                                                     <  2:453781/166‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTC  >  NC_000913/1663‑2145

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: