Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,267,779 G→A 6.2% E223K (GAA→AAA)  tyrB → tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,267,7790GA 6.2% 206.4 / 4.0 65E223K (GAA→AAA) tyrBtyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGAAGTAGAAGAGATGCGTACT  >  NC_000913/4267498‑4268063
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                             
cTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGa                                                                                                                                                                                                                                                                                        >  1:292064/1‑289 (MQ=255)
        gAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTAc                                                                                                                                                                                                                                                                                 >  1:312204/1‑288 (MQ=255)
         aaaaCACCGTAGCAATATTGTCCGGGGCTGGATTCGAAGTGAGAACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACg                                                                                                                                                                                                                                                                                <  2:687098/288‑1 (MQ=255)
                                     gggttacgaacggAGTACTTACCCCTGGTATGACGAACCGACTAACGGGGTGCGCTTTAAGGACCTGTTGGCGACGCTGAAAACCTTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCAGAACCCGACGGGGGCCGTTCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCCAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAGGAGGATGCCTACGCTATTAGCGCCATTGCCAGCGCTGGAtt                                                                                                                                                                                                                                                    <  2:515555/276‑1 (MQ=255)
                                       gATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTAc                                                                                                                                                                                                                                                  >  1:164985/1‑288 (MQ=255)
                                         ttCGAAGGAAGTACTTACCCCTGGTATGGCGAAGCAACTAACGGCGTGCGCTGTAATGACCGGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCg                                                                                                                                                                                                                                               <  2:493517/289‑1 (MQ=255)
                                                         aCCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTc                                                                                                                                                                                                                               >  1:716501/1‑289 (MQ=255)
                                                                         agggACTACCGTCGTGCGCTTTAATGACCTGTTGTCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCCCAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtc                                                                                                                                                                                                               <  2:754090/286‑1 (MQ=255)
                                                                                       tGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGAAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGa                                                                                                                                                                                                 >  1:134429/1‑289 (MQ=255)
                                                                                                 tacccTGTTGGACACCGTGCAAACATTACCTACCCGGAGGTATGGGTTATTGCATCCAGGATGCAACAAACCCACGGGGTCCCATCCGACTAGTGATCAGTGGGATGCAGTGATTCAAATCCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGGGCCGGTATCGAAGAGGATGCCTACGCTATTCGCTCCATTGCCAGCGCGGGATTACCCGCTCTGGTGAGCAATTCGTTCTCTAAAATTTTCCGCCTGTACGGCGAGCGCGTcggc                                                                                                                                                                                       <  2:510151/286‑1 (MQ=255)
                                                                                                    ccTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGAAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTcggcgg                                                                                                                                                                                     <  2:134429/288‑1 (MQ=255)
                                                                                                                                         agtgtgTTGCTGCATCCATGTTGACACAACCCAACGGGTGGCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATGTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGcc                                                                                                                                               <  2:684762/287‑1 (MQ=255)
                                                                                                                                            gtgtTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGAc                                                                                                                                                                                   <  2:285315/250‑1 (MQ=255)
                                                                                                                                            gtgtTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGAc                                                                                                                                                                                   >  1:285315/1‑250 (MQ=255)
                                                                                                                                                  ctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAACTTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTgg                                                                                                                                      >  1:289892/1‑289 (MQ=255)
                                                                                                                                                  ctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTCCTgg                                                                                                                                      >  2:158437/1‑289 (MQ=255)
                                                                                                                                                  ctgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTgg                                                                                                                                      >  2:523128/1‑289 (MQ=255)
                                                                                                                                                       tCCATTTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCCGTGGTATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCaa                                                                                                                                 >  1:309125/1‑289 (MQ=255)
                                                                                                                                                       tCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGCCTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCaa                                                                                                                                 >  1:742160/1‑289 (MQ=255)
                                                                                                                                                       tCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCaa                                                                                                                                 >  1:354068/1‑289 (MQ=255)
                                                                                                                                                         cATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtc                                                                                                                                                                                                               >  2:202021/1‑209 (MQ=255)
                                                                                                                                                         cATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtc                                                                                                                                                                                                               <  1:202021/209‑1 (MQ=255)
                                                                                                                                                                    aaCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTGGGTGCCGTTATGGAAGAGGATGCTTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAACTTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATTGTTGAATGTGCCGAAGCCGCTGGCCGCGTCCTGGGGCAATTGCAAGCAAcga                                                                                                                    >  2:747901/1‑287 (MQ=255)
                                                                                                                                                                     aCCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGt                                                                                                                   >  1:66564/1‑289 (MQ=255)
                                                                                                                                                                         aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcgc                                                                                                              <  2:309125/289‑1 (MQ=255)
                                                                                                                                                                         aaCGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcg                                                                                                               >  1:652521/1‑288 (MQ=255)
                                                                                                                                                                               tGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCaa                                                                                                         <  1:158437/288‑1 (MQ=255)
                                                                                                                                                                                          cTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGccc                                                                                             <  2:292064/289‑1 (MQ=255)
                                                                                                                                                                                           tAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCccg                                                                                            <  2:354068/289‑1 (MQ=255)
                                                                                                                                                                                                   aGTGGGATTCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTCTTCCATTCCTCGACACTGCCTATCAAGGTTCCGGTGCCTGTACGTGAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGCGCAATACGTTCTCGAACTCTCTCTCCCTTCACTGCGATCGCGTCGGCGCACGTTCCGTTATGCGTGAAGATGACGCACCCTCTGGCCGCTTCATGGGGCAATTCAAAGCAGAGTTCCGTCGCAACTGACCCACCCCGCTGCAttt                                                                                    >  2:618250/1‑289 (MQ=255)
                                                                                                                                                                                                   aGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAAttt                                                                                    >  1:364162/1‑289 (MQ=255)
                                                                                                                                                                                                     tGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGa                                                                                                                                                                                                                                                      >  1:184227/1‑126 (MQ=255)
                                                                                                                                                                                                     tGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGa                                                                                                                                                                                                                                                      <  2:184227/126‑1 (MQ=255)
                                                                                                                                                                                                     tGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCa                                                                                                                        <  2:351059/252‑1 (MQ=255)
                                                                                                                                                                                                     tGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCa                                                                                                                        >  1:351059/1‑252 (MQ=255)
                                                                                                                                                                                                       ggATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtctc                                                                                                                                                                                                             <  2:202390/165‑1 (MQ=255)
                                                                                                                                                                                                       ggATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtctc                                                                                                                                                                                                             >  1:202390/1‑165 (MQ=255)
                                                                                                                                                                                                          tGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCg                                                                                                                                                           >  2:571838/1‑212 (MQ=255)
                                                                                                                                                                                                          tGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCg                                                                                                                                                           <  1:571838/212‑1 (MQ=255)
                                                                                                                                                                                                          tGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTgc                                                                              >  1:473221/1‑288 (MQ=255)
                                                                                                                                                                                                          tGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTgc                                                                              >  1:392284/1‑288 (MQ=255)
                                                                                                                                                                                                            cGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGcc                                                                                                                                               >  2:204208/1‑222 (MQ=255)
                                                                                                                                                                                                            cGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGcc                                                                                                                                               <  1:204208/222‑1 (MQ=255)
                                                                                                                                                                                                             ggTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGAGGACTTTATGTGATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAg                                                                          <  1:618250/289‑1 (MQ=255)
                                                                                                                                                                                                                gATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAggt                                                                        <  1:523128/288‑1 (MQ=255)
                                                                                                                                                                                                                     taattgtgagatcgcgcgGGGTTATTCCATTCCTGGGGATTGCCTATCAAGGCGTTGGGGCCGGTATGGAAGAGGCTGCATACGCGATTCGAGCCAGTGCCATGGCTGGATTACCCGGTCGGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTGTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGtggtgg                                                                   <  2:638316/274‑1 (MQ=255)
                                                                                                                                                                                                                       agtctcCAAGCCCTCGAGCTTATTCCATTCCTAGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGTATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGCTTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTg                                                                <  2:742160/287‑1 (MQ=255)
                                                                                                                                                                                                                       aTTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTg                                                                <  2:66564/289‑1 (MQ=255)
                                                                                                                                                                                                                         tctcAAAGCCCGCGAGCTTATTCCATTCCTGGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCa                                                              <  2:392284/289‑1 (MQ=255)
                                                                                                                                                                                                                         tctcAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCa                                                              <  2:473221/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTAc                                                                                                                                                                                                                                                  <  1:28960/85‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTAc                                                                                                                                                                                                                                                  >  2:28960/1‑85 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCa                                                              >  2:42792/1‑264 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCa                                                              <  1:42792/264‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGaaa                                      >  1:395083/1‑288 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGaaa                                      >  2:220093/1‑288 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGaaa                                      >  1:620386/1‑288 (MQ=255)
                                                                                                                                                                                                                                                  ttccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAg                                     >  2:305742/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       cGATATTGCCTTTCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCa                                  >  1:663195/1‑287 (MQ=255)
                                                                                                                                                                                                                                                               ccTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAg                                               >  1:326365/1‑266 (MQ=255)
                                                                                                                                                                                                                                                               ccTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAg                                               <  2:326365/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGAAg                   <  1:747901/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGATTTGGTGCCGGTATGAAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtc                                                                                                                                                                                                               >  2:550140/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGATTTGGTGCCGGTATGAAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTtc                                                                                                                                                                                                               <  1:550140/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aTGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGAAGTAGAAGAGATGCGTACt  <  1:725719/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                             
CTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGAAGTAGAAGAGATGCGTACT  >  NC_000913/4267498‑4268063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: