Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,408,074 A→C I172L (ATC→CTC)  alaA → glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,408,0740AC98.8% 278.9 / ‑3.2 85I172L (ATC→CTC) alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (41/43);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGATCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCTTCCGTCAGG  >  NC_000913/2407841‑2408308
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
ccGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGGGCATTTTCTTTGCGATGAATCCTCTTACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCCCGTACGCGTGGGCCCGtt                                                                                                                                                                                                                                       >  1:25264/1‑239 (MQ=255)
    ggAAGATATTTACATCGGCAATGGTGTTTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGt                                                                                                                                                                                                                                        >  2:123502/1‑234 (MQ=255)
    ggAAGATATTTACATCGGCAATGGTGTTTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGt                                                                                                                                                                                                                                        <  1:123502/234‑1 (MQ=255)
     gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGttatta                                                                                                                                                                                                                                   >  1:701980/1‑238 (MQ=255)
          tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCaa                                                                                                                                                                                                                               <  1:1154957/237‑1 (MQ=255)
             ttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCaa                                                                                                                                                                                                                          <  1:80142/239‑1 (MQ=255)
                cATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCGCGTTATTATCAACCCAAATa                                                                                                                                                                                                                       >  1:20787/1‑239 (MQ=255)
                   cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCCCTGACTGGTTCCCGGACCTCGATGATATTTGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAAACCCAATaaa                                                                                                                                                                                                                     >  1:1004234/1‑237 (MQ=255)
                     gCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCa                                                                                                                                                                                                                  >  1:296611/1‑239 (MQ=255)
                         tGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCg                                                                                                                                                                                                              <  1:1004978/239‑1 (MQ=255)
                             gtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           <  1:311383/238‑1 (MQ=255)
                             gtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           <  1:210256/238‑1 (MQ=255)
                              tATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGAGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGcgcg                                                                                                                                                                                                         >  2:265264/1‑239 (MQ=255)
                               aTCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGcgcg                                                                                                                                                                                                         <  1:555041/238‑1 (MQ=255)
                                      cTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATTTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCCCGTACGCGTGGGCTCGTTATTATCAAACCAAAAAACCCAAACGGCGCGGCATAtt                                                                                                                                                                                                  >  2:854030/1‑238 (MQ=255)
                                              tCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCa                                                                                                                                                                                               >  2:820385/1‑233 (MQ=255)
                                              tCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCa                                                                                                                                                                                               <  1:820385/233‑1 (MQ=255)
                                                          ggCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGGGCATTATCTTTGCGTTGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCTAGCTAAAATTAAGCCTCCTACGCGGCGGCTCGGTTTTTTCTACCCAAAAAAACCAACCGGCGCGGTATATTTCAAAGAGCCTTTAAtggagt                                                                                                                                                                             >  2:1192754/1‑238 (MQ=255)
                                                                     aCAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGGTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGGATATTCCCAAGAGGTTTTAAtggagattgtgggaatt                                                                                                                                                                  >  2:1113305/1‑234 (MQ=255)
                                                                         cGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGGTGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAAc                                                                                                                                                                                                                              >  2:494236/1‑175 (MQ=255)
                                                                         cGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGGTGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAAc                                                                                                                                                                                                                              <  1:494236/175‑1 (MQ=255)
                                                                         cGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCCAATAACCCAACCGGCGCGGTCAATTCCAAAGAGCTTTTAATGGAGATTGTGGAGAGTGCAc                                                                                                                                                              >  2:632387/1‑239 (MQ=255)
                                                                         cGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCAc                                                                                                                                                              >  2:434140/1‑239 (MQ=255)
                                                                         cGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCAc                                                                                                                                                              >  1:94102/1‑239 (MQ=255)
                                                                            ggACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGt                                                                                                                                                            <  1:434140/238‑1 (MQ=255)
                                                                                gAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           <  2:747782/187‑1 (MQ=255)
                                                                                gAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           >  1:747782/1‑187 (MQ=255)
                                                                                  aaTGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATa                                                                                                                                                     <  2:296611/239‑1 (MQ=255)
                                                                                              tGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCCTTATCt                                                                                                                                         >  2:660829/1‑239 (MQ=255)
                                                                                                       ttACCCACTCTGGCCAGCGTCGGTTTCGTTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATg                                                                                                                                <  1:1183891/239‑1 (MQ=255)
                                                                                                             aCTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAAttt                                                                                                                          >  2:863473/1‑239 (MQ=255)
                                                                                                                   gACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACa                                                                                                                    >  2:29795/1‑239 (MQ=255)
                                                                                                                      cgcgGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACaaaa                                                                                                                 <  2:701980/239‑1 (MQ=255)
                                                                                                                        cggcggTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGttatt                                                                                                                                                                                                                                    >  2:185298/1‑122 (MQ=255)
                                                                                                                        cggcggTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGttatt                                                                                                                                                                                                                                    <  1:185298/122‑1 (MQ=255)
                                                                                                                        cggcggTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAAt                                                                                                                >  2:851538/1‑238 (MQ=255)
                                                                                                                         ggcggTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGTTGAATTCCCTGACTGGTTCTCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGCGGGCCCGTTATTATCAACCCAAATAACCCAAACGGCGCGGGATATTCCAAAGAGCTTTTAATGGAGATTTTTGAGATTTCACGCCAGCATAATATCAATAGCTTCGCCGATGAAATTTATGACaaaaat                                                                                                               >  2:124630/1‑236 (MQ=255)
                                                                                                                            ggTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATtctc                                                                                                            <  1:29795/238‑1 (MQ=255)
                                                                                                                                     ttCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTacgacgacg                                                                                                  <  1:1113305/239‑1 (MQ=255)
                                                                                                                                     gtcaaGCGTTAAAGCGTGGCAGTATCATTGCTAGGAATCCGACGACTGGTTCCCGGACCTCTATGATATTCACGATAAAATTACGCCTCGTACGCGTGGGCTGGTTATTATCAACCCAAAGAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAAGGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATATAAGACAAAATTCTCTatgacgacg                                                                                                  <  2:1101024/235‑1 (MQ=255)
                                                                                                                                                     gTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAAt                                                                                  >  1:319684/1‑239 (MQ=255)
                                                                                                                                                      tGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAAtt                                                                                 >  1:650681/1‑239 (MQ=255)
                                                                                                                                                       gCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAAtt                                                                                 >  1:649156/1‑238 (MQ=255)
                                                                                                                                                        cATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           >  2:1128335/1‑115 (MQ=255)
                                                                                                                                                        cATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGgcg                                                                                                                                                                                                           <  1:1128335/115‑1 (MQ=255)
                                                                                                                                                        cATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTgc                                                                               >  1:19544/1‑239 (MQ=255)
                                                                                                                                                              cTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCt                                                                         <  2:948613/239‑1 (MQ=255)
                                                                                                                                                                ttGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGcgcc                                                                            >  2:51376/1‑234 (MQ=255)
                                                                                                                                                                ttGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTgg                                                                       <  1:829057/239‑1 (MQ=255)
                                                                                                                                                                ttGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGAAAAAATTCTCTACGACGACGCTGAGCATCACTCAATTTcgcc                                                                            >  1:548050/1‑234 (MQ=255)
                                                                                                                                                                ttGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGATCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGcgcc                                                                            >  1:575133/1‑234 (MQ=255)
                                                                                                                                                                      tGAATCCTCTGACTGGTTCCCGGGCCTCGATGATATGCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCacct                                                                  <  2:284690/238‑1 (MQ=255)
                                                                                                                                                                                  cTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCa                                                     <  2:650681/239‑1 (MQ=255)
                                                                                                                                                                                     gTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCAtt                                                   >  2:880058/1‑238 (MQ=255)
                                                                                                                                                                                      ttCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGt                                                                                                                                                                                                       <  1:846718/89‑1 (MQ=255)
                                                                                                                                                                                      ttCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGt                                                                                                                                                                                                       >  2:846718/1‑89 (MQ=255)
                                                                                                                                                                                      ttCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCt                                                                         <  2:1064920/215‑1 (MQ=255)
                                                                                                                                                                                      ttCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCt                                                                         >  1:1064920/1‑215 (MQ=255)
                                                                                                                                                                                      ttCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCa                                                     <  1:863473/235‑1 (MQ=255)
                                                                                                                                                                                       tCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTAc                                                 <  2:649156/238‑1 (MQ=255)
                                                                                                                                                                                        cccGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTacgacga                                                                                                    <  2:244823/186‑1 (MQ=255)
                                                                                                                                                                                        cccGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTacgacga                                                                                                    >  1:244823/1‑186 (MQ=255)
                                                                                                                                                                                           ggACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACg                                                                                                                                                             >  2:11029/1‑126 (MQ=255)
                                                                                                                                                                                           ggACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACg                                                                                                                                                             <  1:11029/126‑1 (MQ=255)
                                                                                                                                                                                             aCCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCtt                                                                                                                                                                                       >  1:1101546/1‑98 (MQ=255)
                                                                                                                                                                                             aCCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCtt                                                                                                                                                                                       <  2:1101546/98‑1 (MQ=255)
                                                                                                                                                                                                tCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGa                                       <  1:51376/239‑1 (MQ=255)
                                                                                                                                                                                                tCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGa                                       <  1:124630/239‑1 (MQ=255)
                                                                                                                                                                                                 cGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACgg                                        >  1:369356/1‑237 (MQ=255)
                                                                                                                                                                                                    tgatATTCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGt                                   <  2:795508/239‑1 (MQ=255)
                                                                                                                                                                                                         ttCGCGCTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGaaa                              <  1:265264/239‑1 (MQ=255)
                                                                                                                                                                                                            gcgcTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAAtggagattgtgga                                                                                                                                                                      >  2:1085573/1‑100 (MQ=255)
                                                                                                                                                                                                            gcgcTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAAtggagattgtgga                                                                                                                                                                      <  1:1085573/100‑1 (MQ=255)
                                                                                                                                                                                                             cgcTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTa                                                                                                          <  2:1030444/159‑1 (MQ=255)
                                                                                                                                                                                                             cgcTAAAATTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTa                                                                                                          >  1:1030444/1‑159 (MQ=255)
                                                                                                                                                                                                                 aaaaTTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCCTTACCTTTAACGGACTGTCGAAAGCGTAcc                       >  2:976972/1‑238 (MQ=255)
                                                                                                                                                                                                                 aaaaTTACGCCTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACcg                      <  2:369356/239‑1 (MQ=255)
                                                                                                                                                                                                                      tACGCCTCGTACGCGTTGGCTCGTTATTATCAACCCAAATCACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTTTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTg                 <  2:20787/239‑1 (MQ=255)
                                                                                                                                                                                                                           cTCGTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAgg             <  2:94102/238‑1 (MQ=255)
                                                                                                                                                                                                                              gTACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCtt          >  1:684361/1‑238 (MQ=255)
                                                                                                                                                                                                                               tACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACg                                                                                                                                                             <  1:1196100/90‑1 (MQ=255)
                                                                                                                                                                                                                               tACGCGTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACg                                                                                                                                                             >  2:1196100/1‑90 (MQ=255)
                                                                                                                                                                                                                                 cgcgTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCTTCCGt      <  2:19544/239‑1 (MQ=255)
                                                                                                                                                                                                                                    gTGGGCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCTTCCGTCa    <  2:319684/238‑1 (MQ=255)
                                                                                                                                                                                                                                      gggCTCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCTTCCGTCAgg  <  1:1065205/238‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
CCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTCGCGCTAAAATTACGCCTCGTACGCGTGGGATCGTTATTATCAACCCAAATAACCCAACCGGCGCGGTATATTCCAAAGAGCTTTTAATGGAGATTGTGGAGATTGCACGTCAGCATAATCTCATTATCTTCGCCGATGAAATTTATGACAAAATTCTCTACGACGACGCTGAGCATCACTCAATTGCGCCGCTGGCACCTGACCTGCTGACCATTACCTTTAACGGACTGTCGAAAACGTACCGCGTTGCAGGCTTCCGTCAGG  >  NC_000913/2407841‑2408308

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