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→GTG) thrAfused 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

CGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCT  >  NC_000913/2469‑2954
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
tgAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AACAATGGCGAAAACGCCATGGCCTTTTCTAGCCACTATTATCAGCCGCTGCCGTTGGTGCTgcgggggt                                                                                                                                                                                                                                             >  2:245993/2‑246 (MQ=255)
          gTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCggg                                                                                                                                                                                                                                   <  1:59206/251‑1 (MQ=255)
                  ttttATGGCGAATCTGTCACGGCTCGCCGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGTAGGTGGTGCCTGCCGCGTGAAGATTGCCGACGTGGGTGGTAATGATCCGCTGTTCAAAGGG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGGGCTGCGCgggtaagatcggaaaggcgtcg                                                                                                                                                                                                                           >  2:286157/1‑231 (MQ=255)
                  ttttATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGATAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGata                                                                                                                                                                                                                                            <  1:286157/234‑1 (MQ=255)
                        ggCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTAc                                                                                                                                                                                                                       >  1:110531/1‑249 (MQ=255)
                                      agcTCGACGATCTCGTTACCGCGCGCGTGGCGAAGGCCCGTGATGCAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAGGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATgcagttccatctgcggtggtcttt                                                                                                                                                                                                       >  2:93634/3‑228 (MQ=255)
                                            acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCACTATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                                                                                                                                                 >  1:278142/1‑251 (MQ=255)
                                            acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                                                                                                                                                 <  2:457083/251‑1 (MQ=255)
                                                         cgcgcgcgTGGCGAAGGCCCGGGATGAAGGAAAAGTTTTGCGCTATGTGGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAccc                                                                                                                                                                                    <  2:234362/251‑1 (MQ=255)
                                                              gcgTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCa                                                                                                                                                                               <  2:110531/251‑1 (MQ=255)
                                                                 gggCGAAGGCCGGTGATGAAGGAAAAGTTTTGCGTTATGTGGGCAATATTGATGAAGATGGGGCCTGCCGCTTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAAGGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTAACAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTTCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                                                                                                                                            <  2:26808/250‑1 (MQ=255)
                                                                  ggCGAAGGCCCGTGAGGAAGGAAAAGTGTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                                                                                                                                            <  1:320511/250‑1 (MQ=255)
                                                                       agggccgggatgaAGGAAAAGTTTTGAGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCTGCTGTTCACAGTG‑AAAAATGGCGAAAACGCCCTGGACTTTTATAGCCACTATTATCAGCCGCAGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACCGCTGCCGGTGTCTTTGCCGAGCTGATACGGACCCTCTCATGGCAGtt                                                                                                                                                                       >  2:156832/9‑250 (MQ=255)
                                                                               gatgaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGCAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGGTCTg                                                                                                                                                              >  2:138113/1‑251 (MQ=255)
                                                                                  gaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGCAGATGGTGCCTGCCGCGTGAAGATTGCAGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑CAAACTGGCGAAAAAGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGAGGGCAATGACGGTCCAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGATCTGAc                                                                                                                                                            >  2:456506/1‑250 (MQ=255)
                                                                                    agtgaactttggtcggtATGTTGGAATTATTGTTGAATATGGTGCCGGCCGCGGGAATATGGCTGAAGTGGATGGTAATGATCCGCTGTT‑AAAATGTAAAAAGTACGAAACGGCGCTGTCCTGTTATGGCCACTATTATCAGCAGCTGCCGTTGGTTCTGCTCTTGTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATATGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATg                                                                                                                                                         <  2:12012/235‑1 (MQ=255)
                                                                                          aGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACGTGGTTaa                                                                                                                                                    >  2:85710/1‑250 (MQ=255)
                                                                                           gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                                                                                                                                            >  1:26128/1‑177 (MQ=255)
                                                                                           gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATGGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAAGGGCGAAAACGCCCTGGCCGTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                                                                                                                                            <  2:26128/177‑1 (MQ=255)
                                                                                                      tGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTACCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCg                                                                                                                                       >  1:367696/1‑251 (MQ=255)
                                                                                                          ggCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCtt                                                                                                                                   <  1:408043/251‑1 (MQ=255)
                                                                                                                     tgaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat                                                                                                                        >  1:486221/1‑251 (MQ=255)
                                                                                                                      gaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGGAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat                                                                                                                        >  1:158633/1‑250 (MQ=255)
                                                                                                                       gaGATGGTGCCTGCCGCGTGACGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGGGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGCGTCTGACATGGTTAACGTTTATGCCCCGGCTTCCAGTGCCAata                                                                                                                       >  2:123219/2‑250 (MQ=255)
                                                                                                                               ggcTGCCGCGTGAAGATGGCCGAATGGGATGGTAAGGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAAGGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTTCAGTGCCAATATGAGCGTCg                                                                                                              <  2:163236/249‑1 (MQ=255)
                                                                                                                                           aaGATTGCCGCAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGGCGGGGTTGAtgtgg                                                                                                  >  2:434451/1‑250 (MQ=255)
                                                                                                                                                 gCCGCAGTGGATGGGAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCgggg                                                                                            >  2:302608/1‑251 (MQ=255)
                                                                                                                                                   cGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAACATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACAGCTGACGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATACGAGAGTCGGGTTTGATGTGCTCGGggcg                                                                                          >  2:121437/1‑251 (MQ=255)
                                                                                                                                                             ggTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCAGTTGGTGCTGCGCGGGTATGGTACGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAACCTCTCATGGAAGTTAGGCGTCTGCCATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGAGTCGGGTTTGATGGGATCGGGGCGGCGGTGacac                                                                                >  1:297402/1‑251 (MQ=255)
                                                                                                                                                             cgactcccttccgatctCAAAGTG‑AAAAATGGCGAAAACGCCCTTGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                                                                <  2:225532/235‑1 (MQ=255)
                                                                                                                                                                aatgacccgattgtaaaaGTA‑AAAAATGGCGAAAATGCCCTGGCCTTCTATAGCCACTATGATCAGCCGCTGCCGTTGGTAATGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTGCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCt                                                                              <  2:431870/235‑1 (MQ=255)
                                                                                                                                                                              cAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                                                                >  1:225532/1‑234 (MQ=255)
                                                                                                                                                                               aaaGTG‑AAAAATGGCGAAAACGCCTTGGCCTTTTATATCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGc                                                              <  2:168246/251‑1 (MQ=255)
                                                                                                                                                                               aaaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                                                        >  1:35042/1‑161 (MQ=255)
                                                                                                                                                                               aaaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                                                        <  2:35042/161‑1 (MQ=255)
                                                                                                                                                                                aaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                                                  <  1:354137/166‑1 (MQ=255)
                                                                                                                                                                                aaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                                                  >  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)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
CGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCT  >  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: ‑