Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,760 T→C A808A (GCT→GCC 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,7600TC98.5% 217.8 / NA 67A808A (GCT→GCCthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (35/31);  minor base A (1/0);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTG  >  NC_000913/2520‑3003
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
ctgtGCCGCGCGCGTGGCGAAGGCCGGTGATGAAGGAAAAGGTTTGCGCTATGTTGGCAATATGGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTAAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCTCTATTAGCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACATCTGCCGGTGTCTTTGCCGAtt                                                                                                                                                                                                                                                 <  2:358085/242‑2 (MQ=255)
  ttGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTAc                                                                                                                                                                                                                                           <  1:12261/249‑1 (MQ=255)
             gTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATg                                                                                                                                                                                                                              <  1:506026/251‑1 (MQ=255)
                gagAAGGCCCGTGATGTAGGAAAAGTTGTGCGCTATGGTGGCAATAGTGATGAAGATGGTGCCTGCCGCTTGAAGATTGCCGAAGGGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTGTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCTCGGATATGGTGCGGGCAAGGACGTTACAGCCGCCGTTGTCTTTGCCGATCTTCTACGTTCCCTCTCATGGaa                                                                                                                                                                                                                           <  2:327438/249‑1 (MQ=255)
                     ggCCCGTGATGAAGGAAAAGTTTTGCGATATGTTGGCACTATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTAAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTCTTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAGTGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTAAGTACCCTCTAATGGAAGtt                                                                                                                                                                                                                        >  2:534801/1‑249 (MQ=255)
                        ccGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGCAGTTAgg                                                                                                                                                                                                                     >  2:461513/1‑249 (MQ=255)
                            gatgcAGGAAAAGTTTTGCGCTATGTTGGCACTCTTGCTGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGGGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTCTTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCACTGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGCAGTTAGGAGTCTg                                                                                                                                                                                                               >  2:374618/1‑251 (MQ=255)
                             atgaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTg                                                                                                                                                                                                               >  1:601147/1‑250 (MQ=255)
                               gcaGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGAGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGCCGTTACAGCTGCCGGTGACTTTGCCGCTCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGAc                                                                                                                                                                                                             >  2:407127/3‑250 (MQ=255)
                                       aGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGTTTaaa                                                                                                                                                                                                    <  2:601147/251‑1 (MQ=255)
                                          ttgtctatATGTTGGAAATATTGATGAAGAGGGGGCCTGCCGCGTGAAGATTGCCGAAGGGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGTTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGt                                                                                                                                                                                                  <  2:498391/243‑1 (MQ=255)
                                            tGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATTATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCTGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGttt                                                                                                                                                                                                <  2:370500/250‑1 (MQ=255)
                                                        gCAATATTGATGAAGAGGGTGCCTGCCGAGTGAAGATTGTCGAAGTGGATGGTAATGATCCGCTGTGCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGTTGTCTTTGCCGATCTGCTACGTACCCTCGCATTGAAGGTAGGAGGCTGACATGGTTAAAGTTTATGCCCCGGCTTc                                                                                                                                                                                   <  2:367858/251‑1 (MQ=255)
                                                                 atgaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCa                                                                                                                                                                           >  1:305643/1‑250 (MQ=255)
                                                                      gatggtgccTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGa                                                                                                                                                                     <  2:35346/243‑1 (MQ=255)
                                                                      gatggtgccTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCGTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGa                                                                                                                                                                     <  2:406037/243‑1 (MQ=255)
                                                                              cTGCCGCGTGACGCTTGACGAAGTGGATGGTAATGATACGCTGTTCAAAGTGACACATGGCGAAAAGGACCTGGCCTTTTATAGCCCCTATTATCAGACGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTCCACTCTCATGGAAGTTAGGAGTCTGACTTGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGGCggg                                                                                                                                                             >  2:274040/1‑251 (MQ=255)
                                                                                  cgcgTGAAGATTGCCTAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTGTATAGCCACGATTAGCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTg                                                                                                                                                         <  2:191307/251‑1 (MQ=255)
                                                                                    cgTGAAGATTGCCGGAGTGGCTGGTAATGATTCGCTGTTCAAAGGGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACTTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTGTATGCCCCGGCTTCCATTGCCAATATGAGCGTTGGGTTTGAt                                                                                                                                                       <  2:88181/251‑1 (MQ=255)
                                                                                         aGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAg                                                                                                                                                                                                                      >  1:600672/1‑183 (MQ=255)
                                                                                         aGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAg                                                                                                                                                                                                                      <  2:600672/183‑1 (MQ=255)
                                                                                         aGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCt                                                                                                                                                  >  2:116043/1‑251 (MQ=255)
                                                                                          gagtGCCGAAGTGGCTGGTAATGATCCGCTGTTCAAAGTGAAAAGTGGCGAAAACGGCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTTCCCTCTCATGGAAGTTAGGAGTCTTACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATTTGCGTCGGGTTTGAtgtgttc                                                                                                                                                 <  2:461446/248‑4 (MQ=255)
                                                                                          gATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAAGATGAGCGTCGGGTTTGATGTGCTc                                                                                                                                                 <  2:64017/251‑1 (MQ=255)
                                                                                          aaTTGCCGAAGTTGAGTGTAATGATCCGCTGTACAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCTGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTATGAGTCTGACATGGTTAAAGTGTATGCCCCGGCTTCCAGTGCCAATATGAGCGGCGCGTTTGACGTGCtt                                                                                                                                                 <  2:89774/250‑2 (MQ=255)
                                                                                           aTTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGACAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCAGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCg                                                                                                                                                >  2:294424/1‑251 (MQ=255)
                                                                                              gCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCgggg                                                                                                                                             >  1:560155/1‑251 (MQ=255)
                                                                                                cGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGggcg                                                                                                                                           >  1:471402/1‑251 (MQ=255)
                                                                                                cGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGggcg                                                                                                                                           >  1:267491/1‑251 (MQ=255)
                                                                                                 gAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGggcgg                                                                                                                                          >  1:374446/1‑251 (MQ=255)
                                                                                                       gATGGTAATGATCAGCTGTTCAACGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTg                                                                                                                                     >  2:404839/1‑250 (MQ=255)
                                                                                                        aTGGTAATGATCCGTTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGGCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCTGCGGTGac                                                                                                                                   <  1:136296/251‑1 (MQ=255)
                                                                                                        aTGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCCTGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGac                                                                                                                                   >  1:10388/1‑251 (MQ=255)
                                                                                                             aaTGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTg                                                                                                                              <  1:494236/251‑1 (MQ=255)
                                                                                                              aTGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGt                                                                                                                             <  2:164705/251‑1 (MQ=255)
                                                                                                               tGATCCGCTGTTCAAAGTGAAAATTGGCGAAAACGCCCTGGAATTTTATAGCCCCTATTCTCAGCCGCTGACGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGACGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTCGGCGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGAGGCGGTGACACCTGtt                                                                                                                            >  2:559482/1‑251 (MQ=255)
                                                                                                                     gCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGAGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGt                                                                                                                      >  1:590493/1‑251 (MQ=255)
                                                                                                                        gTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCGTTGCCGATCTGCTACGTACCCTCTCATGGAAGGTATGAGTCTGACATGTTTAAAGGTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCa                                                                                                                   <  1:413118/251‑1 (MQ=255)
                                                                                                                                tGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCgg                                                                                                           >  1:159456/1‑251 (MQ=255)
                                                                                                                                         gCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGCAACCTGTTGATGGTGCATTGCTCGGAGTTGTAGt                                                                                                  >  1:442899/1‑251 (MQ=255)
                                                                                                                                          cGAAAACGCCCTGTCATTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTc                                                                                                 <  2:288070/251‑1 (MQ=255)
                                                                                                                                            aaaaCGCCCTGGGCTGTGATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGAGCGGATATGGTGTGGGCAATGACGCTACAGCTGCCGGTGTCTTTGGCGATCTGCTACGTCCCCTCACATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGGCGCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCACGGGGCGGCGTTGACACCTGTTGATGGTGCATTGCTTGGTGTTGTAGTCAc                                                                                               <  2:520194/251‑1 (MQ=255)
                                                                                                                                                cGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGGGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGtt                                                                                           <  1:243420/251‑1 (MQ=255)
                                                                                                                                                 gCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTg                                                                                          <  2:159456/251‑1 (MQ=255)
                                                                                                                                                     tGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAggc                                                                                      >  1:156756/1‑251 (MQ=255)
                                                                                                                                                           tttataGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCa                                                                                  <  1:461513/247‑1 (MQ=255)
                                                                                                                                                             tataGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCAGTGTCTTTGCCGATCTGCTACTTACCCTCTCATTGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGAGTCGGCGGTGACACCTGTTGCTGGTACATTGCACGTAGATGTAGTCACGGTTGCGGCGTCag                                                                                 >  1:80565/1‑248 (MQ=255)
                                                                                                                                                              ataGCCACTATTCTCAGCCGCTGCAGTTGGTGCTGCGCGGATATGGAGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCCTGGACGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACAACTGTTGATGGGGGATTGCTCGGGGTTGGTGTCACGGTTGAGGCGGCagaggc                                                                             >  2:369610/1‑249 (MQ=255)
                                                                                                                                                                        ttatCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGtctt                                                                   >  1:405373/1‑250 (MQ=255)
                                                                                                                                                                          atCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACCTGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGCGCCTTGCTCGGCGATGTAGTCACGGTTGAGTCGGCAGAGACCTTCAGACTca                                                                  >  2:95315/1‑250 (MQ=255)
                                                                                                                                                                                 gCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGGTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATTTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTc                                                          <  1:231539/251‑1 (MQ=255)
                                                                                                                                                                                    gCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCGCAACAACCTCgg                                                        <  1:192353/250‑1 (MQ=255)
                                                                                                                                                                                    gCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCGTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTGGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCg                                                         <  1:374618/249‑1 (MQ=255)
                                                                                                                                                                                       gatGGTGCTGCGCGGATATGGTACGGGCAATGACGTTACAGCTGCAGGTGTCTTTGCCGATCTGCTCCGTACCCTCTCATGGAAGTTAGGGGTCTGACATGGTTAAAGATTCTGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTAATGGAGCCTTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGGGACATTCAGTCTCAACACCCTCGGCCGc                                                    >  2:499068/3‑251 (MQ=255)
                                                                                                                                                                                       gTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGc                                                    >  2:208023/1‑251 (MQ=255)
                                                                                                                                                                                       gTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGACGGTGTCTTTGCCGATCTGCTACGTACACTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGATTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACCCCTGTTGATGGTGCATTGCTCGGAGATGTAGTAACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACg                                                     >  1:289930/1‑250 (MQ=255)
                                                                                                                                                                                         tGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAAACTCGGACGCtt                                                  >  1:626225/1‑251 (MQ=255)
                                                                                                                                                                                         tGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTAAGTCTCAACAACCTCGGACGCtt                                                  >  1:395792/1‑251 (MQ=255)
                                                                                                                                                                                             gcTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGTGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGGTGATGCGGCAGAGACATTCAGGCTCAACAACCTCGGACTCTTTGcc                                              <  2:541092/250‑1 (MQ=255)
                                                                                                                                                                                              cTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTGAGTCTCAACAACCTCGGGCGCTTTGCCg                                             >  2:170327/1‑251 (MQ=255)
                                                                                                                                                                                                            tGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCACTCATGGAAGTTCGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCAAATATGAGCGTCGGGATTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGCGATGTAGTCACGGATGAGGCGGCTGCGACCTTCAGTCTCAACAAACTCGGGCGCTTTGCCGATAAGCTGccctc                                >  2:285906/1‑247 (MQ=255)
                                                                                                                                                                                                                 gCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGAAGCTTTGCCGATAAGCTGCCGTCAGAAcc                          >  1:373858/1‑251 (MQ=255)
                                                                                                                                                                                                                      gACGTTACAGCTGCCGGTGGCTTTGCCGATCGGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAATGTGTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGAGAAGCTGCCGTCAGAACCACgg                      <  1:499068/250‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTTACAGCTGCCGGTGTCTTTGCAGATCTGCTACGTACCCTCTCATGGAAGTTAGGGGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGATTGATGTGCTCGGGGCGGCGGTGCCACCTGTTGATGGGGCATTGCTCGGAGATGTCGTCACGGTTAAGGCGGCAGAGAAATTCCGTCTTCACAAACTCGGGCGCTTTGCCGATCAGCTGcaggcagcaaaacggg                     >  2:381365/1‑235 (MQ=255)
                                                                                                                                                                                                                                 tGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCCGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGCAAATATCGtt          >  1:310401/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCCGGTGTCGTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATTTGCGCGGGGCGGCGGTGACTCCTGTTGATTGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGTCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGt           <  1:372278/249‑1 (MQ=255)
                                                                                                                                                                                                                                         tCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCCGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGCAAATATCGTTTATCAGTg  >  1:288106/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
CTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTG  >  NC_000913/2520‑3003

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