Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,520 A→T P240P (CCA→CCT thrB → homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,5200AT100.0% 163.7 / NA 54P240P (CCA→CCTthrBhomoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (19/35);  total (19/35)

GACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCT  >  NC_000913/3275‑3750
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
gACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCCCGCCTGCTCTTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTCCTGCCTGGCtt                                                                                                                                                                                                                                   >  2:31102/1‑251 (MQ=255)
gACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCtt                                                                                                                                                                                                                                   <  1:547494/251‑1 (MQ=255)
 aCATCATCAGCCAGCAAGCGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTc                                                                                                                                                                                                                                  >  1:190210/1‑251 (MQ=255)
  catcatCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTCCTGCCTGGCTTcc                                                                                                                                                                                                                                 >  1:608267/1‑251 (MQ=255)
             agcAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggcgcgg                                                                                                                                                                                                                      <  1:382460/251‑1 (MQ=255)
                      cAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCggct                                                                                                                                                                                                             >  1:278036/1‑249 (MQ=255)
                            ttGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGa                                                                                                                                                                                                        <  1:452553/250‑1 (MQ=255)
                               atggggtgCTGTGGGTGCGTGCGTACACGGGGATTAAAGGTTCGACGGCAGAAGCAAGGTGTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCAGTGAGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTAGCGCTGAACCCTTCCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATc                                                                                                                                                                                                    <  2:536996/244‑1 (MQ=255)
                               atgaGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAGCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATc                                                                                                                                                                                                    <  1:141789/251‑1 (MQ=255)
                                 gaGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCg                                                                                                                                                                                                   <  1:633784/250‑1 (MQ=255)
                                    tGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGgcgg                                                                                                                                                                                               >  1:17161/1‑250 (MQ=255)
                                    tGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGGCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCcggcggctgaggacaacggcgc                                                                                                                                                                                               >  1:352715/1‑235 (MQ=255)
                                            ggTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGtttgggcaatcgggacgggcgacc                                                                                                                                                                                       >  2:612735/1‑228 (MQ=255)
                                             ggacggaggaaTCCGGGGATTAAAGTGTCGACGGCTTAAGCCAGGTCTCTGTGCCCGGCTAATTATCGCCGCCAGGATTGCATTGCGCACGTGCGACATCGGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTCAGCTTGCCGCGAAGCTGATGAAAGCTTTTATCGCGGAACCCTACCGTGAACGGTTCCTGCCTGGCTTCCGGCAGGCGCTGCAGGCGTTTGCGAAAATCGGCGAGGTagcgag                                                                                                                                                                                      <  2:28355/241‑1 (MQ=255)
                                                 tGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGAATGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGACGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTGCTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGTGGaaaatgacgcggggagagaaggg                                                                                                                                                                                  >  2:387441/1‑231 (MQ=255)
                                                     gTATCCGGGGATTAAAGTCCCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTAtct                                                                                                                                                                              <  2:17161/251‑1 (MQ=255)
                                                       aTCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCGGGCGGTATctcc                                                                                                                                                                            >  1:400180/1‑251 (MQ=255)
                                                           ggggATTAAAGTCTCGACGGCAGAAGTCAGGGCTATTTTACCGGGGCAGTATCGCCGCCAGGATGGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTAGTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCGAGCGGTATctccggct                                                                                                                                                                        <  2:209199/251‑1 (MQ=255)
                                                             ggatgaacggttcGACGGCAGAAGCCAGGGCTATTTTTCCGGCGCAGTAGCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATctccggctcc                                                                                                                                                                      <  2:400180/240‑1 (MQ=255)
                                                             ggATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCCCGGTAGCGAGCGGGATctccggctcc                                                                                                                                                                      >  1:259656/1‑251 (MQ=255)
                                                                         cGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGGAGCGAGCGGTGTCTCCGGCTCCGGCCCGAccctt                                                                                                                                                          >  2:143563/1‑248 (MQ=255)
                                                                              gCAGAAGCCCGGGCTATTTGACCGGCGCTGGATCGCCGCCAGGAGTGCTTTGCGCATGGGCGACATCTGGCTGGCTTCATTCACGCCGGCTATGCCCGTCAGCCTGAGCTTGCTGCGAAGTTGATGAAAGATGTGAGCGCTGAACCCGACCGTGAACGGTTACTGCCTGGCTTCCGGCATGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCGCGACCGTGTTCGc                                                                                                                                                     <  2:611675/251‑1 (MQ=255)
                                                                                   aGCCAGGGCTATTTTCCCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACAGCTGGCTGGCTTCATTCTCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCtgt                                                                                                                                                <  2:580161/251‑1 (MQ=255)
                                                                                                       cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGcc                                                                                                                            <  1:612735/251‑1 (MQ=255)
                                                                                                              ccgccgcCAGGATTTCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGcgcg                                                                                                                     <  2:239528/250‑1 (MQ=255)
                                                                                                                     cAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGa                                                                                                              <  1:199385/251‑1 (MQ=255)
                                                                                                                     cAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATAGGCGCGGTATCGAGCGGTATCTCCGGCTCCGGCCCGACCTTTTTCGCTCTTTGAGACAAACCGGATACCGCCCAGCGCGTTGCCgc                                                                                                              >  1:24179/1‑250 (MQ=255)
                                                                                                                           tgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcgg                                                                                                                                                                                                                               <  2:488172/132‑1 (MQ=255)
                                                                                                                           tgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcgg                                                                                                                                                                                                                               >  1:488172/1‑132 (MQ=255)
                                                                                                                              attgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                     >  1:222923/1‑251 (MQ=255)
                                                                                                                               ttgcGCACGGGCGACAGTTGGATGGCTTCAGTCAAGCCTGCTATTCCCGACAGCCTGAGTTTGCCGCGAAGCTGATGAAAGATGGTATCGCTGCACCCTGCCGTGAACGGTTACTGCCTGGCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTggg                                                                                                     <  2:65022/250‑1 (MQ=255)
                                                                                                                                 gcgcCCGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTACTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACATTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTaa                                                                                                  <  1:502385/251‑1 (MQ=255)
                                                                                                                                     aCGGGCGAGATCTGGTTGGGTTCATTCACGCCTGCTATCCCCGTCAGCCTGTGCTGTCCGCGAAGCTGATGAAAGATGTTAACGCGGAACCCGAGCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGGAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTTACAAGCCGGATACCGCCCAGCGCGTTGCCTACTGGTTGGGTAAGAAc                                                                                              <  2:292297/251‑1 (MQ=255)
                                                                                                                                                  ggctggCTTCATTCACGCCTGCTATTCCCGTTAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATc                                                                                 <  1:236460/247‑1 (MQ=255)
                                                                                                                                                                  gCCTGCTATTCCCGTCAGCCTGAGCTGGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCa                                                                 <  1:139006/251‑1 (MQ=255)
                                                                                                                                                                     tGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGTAGCGAGCGGTATCTCCCGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTTGGTAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCatat                                                              >  1:365009/1‑251 (MQ=255)
                                                                                                                                                                        ttttCCCGGCAGCCTTGGCTTGCCGCGAATCTGATGAAAGAGTGTATCGCTGAACCCTTCCGTGAACGGTTACTGCCTTGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGTAATCGGCGCTGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTGCCTGCAAAATCAGGAAGGTTTTGTTCATAttt                                                            <  2:409996/248‑1 (MQ=255)
                                                                                                                                                                            cccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCgg                                                       <  2:416242/251‑1 (MQ=255)
                                                                                                                                                                             ccGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGc                                                      <  1:568262/251‑1 (MQ=255)
                                                                                                                                                                                       gAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTACGGCTCCGGCCCGCCCTTGTTCGCTCTGTGTGGAAAGCAGGATACCTCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTAACTGCCAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAcggc                                            >  1:45209/1‑251 (MQ=255)
                                                                                                                                                                                            tGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcgggc                                        <  1:218711/250‑1 (MQ=255)
                                                                                                                                                                                            tGCCGCGAAGCTGATGAAAGATGGTGGCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCTGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCGCCGGCTCCGGCCCGACCTTGTTCGCTCTGGGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcgggc                                        <  2:573873/250‑1 (MQ=255)
                                                                                                                                                                                                  gAAGCTGATGAAAGATGTTATCGCTTAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAg                                 <  1:439399/251‑1 (MQ=255)
                                                                                                                                                                                                  gAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAACCAGCGGGAACTCCGGCTCCGGCCCGGCCTTGTTCGCTCTGTTGGACAAACCGGGTACCGCCCCGCGCGTTTGCGACCGGTTGGGTAAGAAATACCCGCCAAATCAGGGAGGGTTTGGTCATATTTGCCgggcgggtacgggggggggaagg                                  >  2:637205/1‑229 (MQ=255)
                                                                                                                                                                                                   aaGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCGAGt                                <  1:93018/251‑1 (MQ=255)
                                                                                                                                                                                                    aGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATTGGCGCGGTAGCGAGAGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACCAGACGGAAACCGCCCCGCGCGGTGCCGCCTGGATGGGGAAGAACATCCTGCCAACACAGGAAGGGTTTGTTCATATTTTCCGGCTGGGAACggcggggagagcg                                  >  2:315958/1‑241 (MQ=255)
                                                                                                                                                                                                        gatgaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACt                             <  1:303490/249‑1 (MQ=255)
                                                                                                                                                                                                           gaAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaa                        <  1:330807/251‑1 (MQ=255)
                                                                                                                                                                                                              aGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACt                     <  1:186893/251‑1 (MQ=255)
                                                                                                                                                                                                                   ttATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAACCAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGGCTGGGTGGGGAAGAACTACCTGCAAAAACAGGAAGGGTTTTTTCATATTTGCCGGCTGGATACGGGGGGGGCAcagggacatgaaaaatgaaa                 >  1:300655/1‑231 (MQ=255)
                                                                                                                                                                                                                               cccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACaa     <  1:192179/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccgaCCGTGTAAGGTTACTGCCTGGCGTCCGGCAGGCGCGGCCGTCGGTTGCGTAAATCGGCGCGGTAGCTAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTTGTTGTGTAAGACCTACCTGCAAAATCAGGAAGGTTTTGATCATATTTGCCGGCTGGATACGGCGGGCTCTCGAGTACTGGAAAACTAAATGAAACTCTACAAt    <  2:345240/248‑1 (MQ=255)
                                                                                                                                                                                                                                ccTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt    <  1:248952/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aCCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACTGCGGGCGCACGAGTACTGGAAAACTAAATGAAACGCTACAATCt  <  1:553685/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
GACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGATGTTATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCT  >  NC_000913/3275‑3750

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