Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,497,026 A→T 5.1% intergenic (+34/+31) ypdK → / ← alaC inner membrane protein/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,497,0260AT 5.1% 405.0 / 9.9 117intergenic (+34/+31)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (50/61);  new base T (2/4);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 117 total aligned reads displayed.

AAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGT  >  NC_000913/2496879‑2497171
                                                                                                                                                   |                                                                                                                                                 
aaGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAg                                                                                                                                                  >  4:775454/1‑149 (MQ=255)
 aGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAgg                                                                                                                                                 >  4:194202/1‑149 (MQ=255)
  gTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                                >  2:997697/1‑149 (MQ=255)
    cATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTc                                                                                                                                              <  1:597643/149‑1 (MQ=255)
     aTTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTcc                                                                                                                                             <  1:997697/149‑1 (MQ=255)
      ttCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCt                                                                                                                                            >  2:228683/1‑149 (MQ=255)
      ttCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCt                                                                                                                                            >  2:145635/1‑149 (MQ=255)
      ttCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCt                                                                                                                                            <  1:834733/149‑1 (MQ=255)
         aGGCGTAACAGTAAATACGCGGAGAGATGTATAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCAACGTTTGGGCCGGTACCTTTGCCCTGATGATCTAACAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGtt                                                                                                                                         <  3:903056/149‑1 (MQ=255)
          ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGttt                                                                                                                                        <  2:268522/149‑1 (MQ=255)
             gTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttttc                                                                                                                                     <  1:993581/149‑4 (MQ=255)
               aaCAGTAATTACGCGGAGAGATGTAAAGTGAAATAATTCTTTATGGGAATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGTCGCCATCAGGCTCTTGTTTGCAtt                                                                                                                                   >  1:148378/1‑149 (MQ=255)
                aCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCAttt                                                                                                                                  >  1:411964/1‑149 (MQ=255)
                  aGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCt                                                                                                                                <  3:981654/149‑1 (MQ=255)
                    tAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTga                                                                                                                              <  3:194202/149‑1 (MQ=255)
                           gcgGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgt                                                                                                                       >  1:25628/1‑149 (MQ=255)
                           gcgGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgt                                                                                                                       <  2:829753/149‑1 (MQ=255)
                             ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtgt                                                                                                                     >  1:526042/1‑149 (MQ=255)
                             ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTGTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTTCTGTTTGCATTTCTGACGACGgtgt                                                                                                                     <  2:856118/149‑1 (MQ=255)
                                   aTGTAAAGTGAAATATTTCTTTGTGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGGCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTc                                                                                                               >  3:59136/1‑149 (MQ=255)
                                     gTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcg                                                                                                             >  1:1077147/1‑149 (MQ=255)
                                     gTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcg                                                                                                             >  3:25250/1‑149 (MQ=255)
                                      tAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcgc                                                                                                            <  1:932081/149‑1 (MQ=255)
                                          gTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGttt                                                                                                        <  3:392575/149‑1 (MQ=255)
                                           tGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGtttt                                                                                                       >  1:848427/1‑149 (MQ=255)
                                            gAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTc                                                                                                      <  4:153719/149‑1 (MQ=255)
                                              aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                    <  1:761848/149‑1 (MQ=255)
                                              aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                    <  3:775454/149‑1 (MQ=255)
                                              aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                    <  1:613318/149‑1 (MQ=255)
                                                tatTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGa                                                                                                  <  4:347144/149‑1 (MQ=255)
                                                 atTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGaa                                                                                                 <  2:848427/149‑1 (MQ=255)
                                                 atTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGaa                                                                                                 <  1:415081/149‑1 (MQ=255)
                                                      tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGt                                                                                            >  3:172145/1‑149 (MQ=255)
                                                      tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGt                                                                                            >  2:1095907/1‑149 (MQ=255)
                                                        tATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTTCTTTTGCCCTGATGATCTAAAGGCATACTCTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTGGCATTTCTGACGACGGTGTGTATTCCGCGTTTTCGTGAATATGttt                                                                                          >  4:565564/1‑149 (MQ=255)
                                                        tATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGttt                                                                                          <  1:228683/149‑1 (MQ=255)
                                                          tGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgc                                                                                        <  1:145635/149‑1 (MQ=255)
                                                           gggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTcc                                                                                                              >  3:380913/1‑126 (MQ=255)
                                                            ggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgctg                                                                                      >  1:1033746/1‑149 (MQ=255)
                                                                tttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGttt                                                                                          <  3:671264/141‑1 (MQ=255)
                                                                tttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGttt                                                                                          >  4:671264/1‑141 (MQ=255)
                                                                 ttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGc                                                                                 <  2:411964/149‑1 (MQ=255)
                                                                 ttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGc                                                                                 >  4:565519/1‑149 (MQ=255)
                                                                 ttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGc                                                                                 <  2:589629/149‑1 (MQ=255)
                                                                    ttttATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCggg                                                                              >  3:394968/1‑149 (MQ=255)
                                                                      ttATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgta                                                                            <  4:734625/149‑1 (MQ=255)
                                                                           gTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCGATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTaaa                                                                       <  2:671832/149‑1 (MQ=255)
                                                                             cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAAcc                                                                     >  3:507196/1‑149 (MQ=255)
                                                                             cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAAcc                                                                     <  1:101859/149‑1 (MQ=255)
                                                                              aTCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGa                                                                                                  <  2:484404/119‑1 (MQ=255)
                                                                              aTCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGa                                                                                                  >  1:484404/1‑119 (MQ=255)
                                                                              aTCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGttt                                                                                          >  4:1046423/1‑127 (MQ=255)
                                                                              aTCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCg                                                                    >  4:713812/1‑149 (MQ=255)
                                                                               tCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtt                                                                                                                                         >  2:799144/1‑79 (MQ=38)
                                                                               tCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtt                                                                                                                                         <  1:799144/79‑1 (MQ=38)
                                                                                cGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTc                                                                  <  2:802503/149‑1 (MQ=255)
                                                                                     gggCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGGCGACGGTGTTTATTACGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCAc                                                             >  2:909073/1‑149 (MQ=255)
                                                                                     gggCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCAc                                                             <  3:1072783/149‑1 (MQ=255)
                                                                                         cGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGTACGGaa                                                         >  1:539324/1‑149 (MQ=255)
                                                                                               ttttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                   >  4:809348/1‑149 (MQ=255)
                                                                                               ttttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                   >  3:683159/1‑149 (MQ=255)
                                                                                               ttttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                   >  1:435124/1‑149 (MQ=255)
                                                                                               ttttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                   >  4:44577/1‑149 (MQ=255)
                                                                                                tttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCt                                                  <  4:372634/149‑1 (MQ=255)
                                                                                                 ttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTTAATATGTTTGCTGCTGGCGGGTAGTAATCCGACGGCACGGAACATCGCtt                                                 <  4:1076199/149‑1 (MQ=255)
                                                                                                    cccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTaa                                              <  2:1033746/149‑1 (MQ=255)
                                                                                                       tgatgaTCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                           >  1:974816/1‑149 (MQ=255)
                                                                                                       tgatgaTCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                           <  2:264715/149‑1 (MQ=255)
                                                                                                        gatgatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATccc                                          >  1:812261/1‑149 (MQ=255)
                                                                                                        gatgatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATccc                                          <  2:979767/149‑1 (MQ=255)
                                                                                                         atgatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCct                                         <  1:560821/149‑1 (MQ=255)
                                                                                                           gatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCctct                                       >  1:544804/1‑149 (MQ=255)
                                                                                                           gatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCctct                                       <  1:196441/149‑1 (MQ=255)
                                                                                                             tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTgac                                                                                                                             >  4:367384/1‑61 (MQ=255)
                                                                                                             tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTgac                                                                                                                             <  3:367384/61‑1 (MQ=255)
                                                                                                             tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGa                                     <  2:208349/149‑1 (MQ=255)
                                                                                                              cTAAAAGCAACCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGAGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCGCTGAt                                    <  2:812411/149‑1 (MQ=255)
                                                                                                               tAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                   >  2:997374/1‑149 (MQ=255)
                                                                                                               tAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                   >  3:862422/1‑149 (MQ=255)
                                                                                                                 aaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGc                                 >  1:551793/1‑149 (MQ=255)
                                                                                                                  aaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGcc                                >  2:26020/1‑149 (MQ=255)
                                                                                                                   aGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCt                               <  1:16229/149‑1 (MQ=255)
                                                                                                                            tgCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATAc                      <  3:809348/149‑1 (MQ=255)
                                                                                                                            tgCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATAc                      <  3:44577/149‑1 (MQ=255)
                                                                                                                               aaTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgta                                                                            <  3:834871/92‑1 (MQ=255)
                                                                                                                               aaTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgta                                                                            >  4:834871/1‑92 (MQ=255)
                                                                                                                               aaTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCctc                                        <  3:805152/128‑1 (MQ=255)
                                                                                                                               aaTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCctc                                        >  4:805152/1‑128 (MQ=255)
                                                                                                                                aTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTc                  <  1:1051063/149‑1 (MQ=255)
                                                                                                                                  gAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGGCGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTcgc                <  1:909073/149‑1 (MQ=255)
                                                                                                                                     aaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTTTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTc                  >  2:78290/1‑144 (MQ=255)
                                                                                                                                     aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAAt                        >  2:556898/1‑138 (MQ=255)
                                                                                                                                     aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAAt                        <  1:556898/138‑1 (MQ=255)
                                                                                                                                     aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGt             >  4:1021721/1‑149 (MQ=255)
                                                                                                                                       aCAGGAGCCATCAGGCTCCTGTTTGCTTTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGttt           <  2:417542/149‑1 (MQ=255)
                                                                                                                                          ggAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCa        >  3:659971/1‑149 (MQ=255)
                                                                                                                                            aGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAAt      <  2:1043534/149‑1 (MQ=255)
                                                                                                                                             gCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATc     <  4:683159/149‑1 (MQ=255)
                                                                                                                                                aTCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGt  <  3:418669/149‑1 (MQ=255)
                                                                                                                                                aTCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGt  <  1:430667/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGT  >  NC_000913/2496879‑2497171

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