Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,676 C→T F780F (TTC→TTT 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,6760CT95.0% 171.8 / ‑1.6 60F780F (TTC→TTTthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (32/25);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCG  >  NC_000913/2433‑2920
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
tgaaattgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTTGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGc                                                                                                                                                                                                                                                <  1:249885/250‑1 (MQ=255)
       gaaCCTGTGCTAGCCGCAGAGTTTAACGAGGAGGTGGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCCCGTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGGAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACtatta                                                                                                                                                                                                                                        <  2:28912/251‑1 (MQ=255)
              tgctgcCGGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGACTCTGTCCCAGCACGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGGTTGCCGAAGTGGATGGTAATGGACCGCTGTTCACAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCg                                                                                                                                                                                                                                 >  2:452474/1‑251 (MQ=255)
              tgctgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGGCGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAGAACGACCTGGCCTTTTATAGCCACTATTCTCAGCCg                                                                                                                                                                                                                                 >  1:131538/1‑251 (MQ=255)
                 tgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTg                                                                                                                                                                                                                              >  1:457083/1‑251 (MQ=255)
                    ccGCAGAGTTTAAAGCTGAGGGTGATGTTGCCGCTTTTATGGCGCATCTGTCACAGCTCGACGGTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGGACAGTTTTGCGATATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGACGAAGTGGATGGTAATGATCCGCTGTTCCAAGGGACAAATGGCGCAAACTCCCTGGCCTTTTTGAGCCACTATTAACAGcagactgcg                                                                                                                                                                                                                           >  2:215774/1‑243 (MQ=255)
                        agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTgg                                                                                                                                                                                                                       <  1:173796/251‑1 (MQ=255)
                        agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTgg                                                                                                                                                                                                                       >  1:475592/1‑251 (MQ=255)
                        agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCGGGCCTTTTATAACCACTATTATCAGCCGCGGCCGGTgg                                                                                                                                                                                                                       >  2:100320/1‑251 (MQ=255)
                                    tgAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAACAATGGCGAAAACGCCATGGCCTTTTCTAGCCACTATTATCAGCCGCTGCCGTTGGTGCTgcgggggt                                                                                                                                                                                                           >  2:245993/2‑246 (MQ=255)
                                              gTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCggg                                                                                                                                                                                                 <  1:59206/251‑1 (MQ=255)
                                                      ttttATGGCGAATCTGTCACGGCTCGCCGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGTAGGTGGTGCCTGCCGCGTGAAGATTGCCGACGTGGGTGGTAATGATCCGCTGTTCAAAGGGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGGGCTGCGCgggtaagatcggaaaggcgtcg                                                                                                                                                                                         >  2:286157/1‑231 (MQ=255)
                                                      ttttATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGATAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGata                                                                                                                                                                                                          <  1:286157/234‑1 (MQ=255)
                                                            ggCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTAc                                                                                                                                                                                     >  1:110531/1‑249 (MQ=255)
                                                                          agcTCGACGATCTCGTTACCGCGCGCGTGGCGAAGGCCCGTGATGCAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAGGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATgcagttccatctgcggtggtcttt                                                                                                                                                                     >  2:93634/3‑228 (MQ=255)
                                                                                acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCACTATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                                                                                                               >  1:278142/1‑251 (MQ=255)
                                                                                acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                                                                                                               <  2:457083/251‑1 (MQ=255)
                                                                                             cgcgcgcgTGGCGAAGGCCCGGGATGAAGGAAAAGTTTTGCGCTATGTGGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAccc                                                                                                                                                  <  2:234362/251‑1 (MQ=255)
                                                                                                  gcgTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCa                                                                                                                                             <  2:110531/251‑1 (MQ=255)
                                                                                                     gggCGAAGGCCGGTGATGAAGGAAAAGTTTTGCGTTATGTGGGCAATATTGATGAAGATGGGGCCTGCCGCTTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAAGGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTAACAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTTCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                                                                                                          <  2:26808/250‑1 (MQ=255)
                                                                                                      ggCGAAGGCCCGTGAGGAAGGAAAAGTGTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                                                                                                          <  1:320511/250‑1 (MQ=255)
                                                                                                           agggccgggatgaAGGAAAAGTTTTGAGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCTGCTGTTCACAGTGAAAAATGGCGAAAACGCCCTGGACTTTTATAGCCACTATTATCAGCCGCAGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACCGCTGCCGGTGTCTTTGCCGAGCTGATACGGACCCTCTCATGGCAGtt                                                                                                                                     >  2:156832/9‑250 (MQ=255)
                                                                                                                   gatgaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGCAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGGTCTg                                                                                                                            >  2:138113/1‑251 (MQ=255)
                                                                                                                      gaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGCAGATGGTGCCTGCCGCGTGAAGATTGCAGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGCAAACTGGCGAAAAAGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGAGGGCAATGACGGTCCAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGATCTGAc                                                                                                                          >  2:456506/1‑250 (MQ=255)
                                                                                                                        agtgaactttggtcggtATGTTGGAATTATTGTTGAATATGGTGCCGGCCGCGGGAATATGGCTGAAGTGGATGGTAATGATCCGCTGTTAAAATGTAAAAAGTACGAAACGGCGCTGTCCTGTTATGGCCACTATTATCAGCAGCTGCCGTTGGTTCTGCTCTTGTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATATGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATg                                                                                                                       <  2:12012/235‑1 (MQ=255)
                                                                                                                              aGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACGTGGTTaa                                                                                                                  >  2:85710/1‑250 (MQ=255)
                                                                                                                               gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                                                                                                          >  1:26128/1‑177 (MQ=255)
                                                                                                                               gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATGGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAAGGGCGAAAACGCCCTGGCCGTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                                                                                                          <  2:26128/177‑1 (MQ=255)
                                                                                                                                          tGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTACCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCg                                                                                                     >  1:367696/1‑251 (MQ=255)
                                                                                                                                              ggCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCtt                                                                                                 <  1:408043/251‑1 (MQ=255)
                                                                                                                                                         tgaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat                                                                                      >  1:486221/1‑251 (MQ=255)
                                                                                                                                                          gaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGGAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat                                                                                      >  1:158633/1‑250 (MQ=255)
                                                                                                                                                           gaGATGGTGCCTGCCGCGTGACGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGGGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGCGTCTGACATGGTTAACGTTTATGCCCCGGCTTCCAGTGCCAata                                                                                     >  2:123219/2‑250 (MQ=255)
                                                                                                                                                                   ggcTGCCGCGTGAAGATGGCCGAATGGGATGGTAAGGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAAGGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTTCAGTGCCAATATGAGCGTCg                                                                            <  2:163236/249‑1 (MQ=255)
                                                                                                                                                                               aaGATTGCCGCAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGGCGGGGTTGAtgtgg                                                                >  2:434451/1‑250 (MQ=255)
                                                                                                                                                                                     gCCGCAGTGGATGGGAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCgggg                                                          >  2:302608/1‑251 (MQ=255)
                                                                                                                                                                                       cGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAACATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACAGCTGACGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATACGAGAGTCGGGTTTGATGTGCTCGGggcg                                                        >  2:121437/1‑251 (MQ=255)
                                                                                                                                                                                                 ggTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCAGTTGGTGCTGCGCGGGTATGGTACGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAACCTCTCATGGAAGTTAGGCGTCTGCCATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGAGTCGGGTTTGATGGGATCGGGGCGGCGGTGacac                                              >  1:297402/1‑251 (MQ=255)
                                                                                                                                                                                                 cgactcccttccgatctCAAAGTGAAAAATGGCGAAAACGCCCTTGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                              <  2:225532/235‑1 (MQ=255)
                                                                                                                                                                                                    aatgacccgattgtaaaaGTAAAAAATGGCGAAAATGCCCTGGCCTTCTATAGCCACTATGATCAGCCGCTGCCGTTGGTAATGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTGCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCt                                            <  2:431870/235‑1 (MQ=255)
                                                                                                                                                                                                                  cAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                              >  1:225532/1‑234 (MQ=255)
                                                                                                                                                                                                                   aaaGTGAAAAATGGCGAAAACGCCTTGGCCTTTTATATCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGc                            <  2:168246/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaaGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                      >  1:35042/1‑161 (MQ=255)
                                                                                                                                                                                                                   aaaGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                      <  2:35042/161‑1 (MQ=255)
                                                                                                                                                                                                                    aaGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                >  2:354137/1‑166 (MQ=255)
                                                                                                                                                                                                                    aaGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                <  1:354137/166‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                         >  2:136578/1‑152 (MQ=255)
                                                                                                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                         <  1:136578/152‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTg                                                                                           <  2:255686/182‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTg                                                                                           >  1:255686/1‑182 (MQ=255)
                                                                                                                                                                                                                           aaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGgagg                     >  1:24703/1‑250 (MQ=255)
                                                                                                                                                                                                                            aaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCgg                        >  1:272209/1‑247 (MQ=255)
                                                                                                                                                                                                                               tgtcGAAAACGCCCTGTCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCTCGGATATGGTGCGGTCAATGACGTTACAGCTTCAGGTGTCTTTGCCGATCTGCTACTTACCCTCTCATGTAAGTTAGTACTCTCACATTGTTAAAGTTTAt                                                                                                           <  2:220924/158‑1 (MQ=255)
                                                                                                                                                                                                                               tGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTAt                                                                                                           >  1:220924/1‑161 (MQ=255)
                                                                                                                                                                                                                                ggCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAg                <  1:245993/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGt               <  1:275111/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaaCGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCTGAGATGTAGTCACg           <  1:123219/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTg       >  1:226037/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gggCCTTTTATAGCCACTATTATCAGCCGCTGCCGTGGGTGCTGCGCGGATAGGGTGCGGGCAATGACGTTACAGCTTCCGGTGTCTTTGCCGATCTGTTACGTCCCCTCTCATGGAAGTTAGGAGTCGGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGAGGTGCTCGGGGCGGCGGTGACACCTGTTGAGGGTGCATTGCTCGGAGATTTAGTCACGGTTGAggc   <  2:287142/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAggcg  <  1:1014/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCG  >  NC_000913/2433‑2920

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: