Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,709 A→G V791V (GTA→GTG 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,7090AG94.5% 154.6 / ‑3.0 56V791V (GTA→GTGthrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (27/25);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

CGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCT  >  NC_000913/2469‑2954
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
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                                                                                                                                                  <  1:354137/166‑1 (MQ=255)
                                                                                                                                                                                aaGTGAAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                                                  >  2:354137/1‑166 (MQ=255)
                                                                                                                                                                                      aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                                                           <  1:136578/152‑1 (MQ=255)
                                                                                                                                                                                      aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                                                           >  2:136578/1‑152 (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)
                                                                                                                                                                                                                             tatCAGCCGTGGCCGTTGGTGCTGCGCGGAGATGGTGCGGGCAATGCCGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTCCCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGTCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTTCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGCCATTCAGcctc                  <  2:70531/250‑1 (MQ=255)
                                                                                                                                                                                                                             tatCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGtctc                  >  1:3633/1‑250 (MQ=255)
                                                                                                                                                                                                                                cAGCCGCTGCCGTTGGTGCTGCGCGGATATGGGGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATACGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTcaaca              <  1:452474/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gCCGTTGGGGCTGCGCGGAGAGGGGGCGGGCAATGACGTTACAGCTGCCTGTGTCTGTGCCGATCTGCTACGGACCCCCTCATGGAAGTTAGGATTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCGCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCg        <  2:291734/249‑1 (MQ=255)
                                                                                                                                                                                                                                             tGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCt  <  1:93634/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
CGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTGAAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCT  >  NC_000913/2469‑2954

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: