Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,497,023 A→T 6.3% intergenic (+31/+34) 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,0230AT 6.3% 269.5 / 7.5 79intergenic (+31/+34)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base T (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

AATGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAG  >  NC_000913/2496875‑2497170
                                                                                                                                                    |                                                                                                                                                   
aaTGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCa                                                                                                                                                     >  2:478096/1‑149 (MQ=255)
 aTGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATc                                                                                                                                                   >  1:863126/1‑149 (MQ=255)
   gAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATc                                                                                                                                                   <  3:496187/148‑1 (MQ=255)
      gTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAgg                                                                                                                                                <  4:323935/148‑1 (MQ=255)
      gTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                               >  1:344856/1‑149 (MQ=255)
           tCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCt                                                                                                                                           >  3:1006305/1‑148 (MQ=255)
             aGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGttt                                                                                                                                       <  1:793205/149‑1 (MQ=255)
              ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttt                                                                                                                                       <  2:128095/149‑1 (MQ=255)
                 gTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCa                                                                                                                                                  <  1:397708/135‑1 (MQ=255)
                 gTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCa                                                                                                                                                  >  2:397708/1‑135 (MQ=255)
                   aaCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCAt                                                                                                                                   >  3:212698/1‑148 (MQ=255)
                                 ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtg                                                                                                                     <  4:355828/148‑1 (MQ=255)
                                 ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtg                                                                                                                     <  4:625032/148‑1 (MQ=255)
                                    gagaTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGttt                                                                                                                  <  3:758863/148‑1 (MQ=255)
                                       aTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTc                                                                                                              >  2:654346/1‑149 (MQ=255)
                                            aaGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtt                                                                                                                                        <  1:754882/118‑1 (MQ=255)
                                            aaGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtt                                                                                                                                        >  2:754882/1‑118 (MQ=255)
                                              gTGAAATATTTCTTTATGGGCATTTCTTGTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGttt                                                                                                       >  2:922459/1‑149 (MQ=255)
                                              gTGAAATATTTCTTTACGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGtt                                                                                                        >  4:907901/1‑148 (MQ=255)
                                                  aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCg                                                                                                    <  3:401751/148‑1 (MQ=255)
                                                  aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                   >  2:163218/1‑149 (MQ=255)
                                                  aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                   <  2:857041/149‑1 (MQ=255)
                                                    tatTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTgacga                                                                                                                          >  4:513625/1‑124 (MQ=255)
                                                    tatTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTgacga                                                                                                                          <  3:513625/124‑1 (MQ=255)
                                                      tttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGaa                                                                                                <  4:472427/148‑1 (MQ=255)
                                                       ttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAata                                                                                              <  2:444105/149‑1 (MQ=255)
                                                              tGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgc                                                                                       >  2:542832/1‑149 (MQ=255)
                                                              tGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtg                                                                                        >  3:315745/1‑148 (MQ=255)
                                                               gggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgct                                                                                      >  1:486136/1‑149 (MQ=255)
                                                                ggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgctg                                                                                     <  2:486136/149‑1 (MQ=255)
                                                                    tttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGctgctg                                                                                  >  4:927204/1‑148 (MQ=255)
                                                                     ttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTgg                                                                                 >  4:1000679/1‑148 (MQ=255)
                                                                      tCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGc                                                                                >  4:993057/1‑148 (MQ=255)
                                                                        ttttATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCgg                                                                              <  4:849628/148‑1 (MQ=255)
                                                                         tttATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCggg                                                                             >  3:839099/1‑148 (MQ=255)
                                                                           tATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTTCATTTCTGACGACGGTGTTTTTTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgtag                                                                          >  1:693619/1‑149 (MQ=255)
                                                                            aTGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgtag                                                                          >  3:933599/1‑148 (MQ=255)
                                                                              ggTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTaa                                                                       >  2:51402/1‑149 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTg                                                                                                                                          <  3:692143/79‑1 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                  >  4:233012/1‑119 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                  <  3:233012/119‑1 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                  <  3:977563/119‑1 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                  >  4:977563/1‑119 (MQ=255)
                                                                                 cATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCACCTg                                                                                                                                          >  4:692143/1‑79 (MQ=255)
                                                                                   tCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCg                                                                   <  3:580782/148‑1 (MQ=255)
                                                                                    cGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTc                                                                 <  1:325321/149‑1 (MQ=255)
                                                                                     gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCg                                                                <  2:991753/149‑1 (MQ=255)
                                                                                      tttGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCg                                                                >  4:325429/1‑148 (MQ=255)
                                                                                         gggCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCa                                                             <  3:903457/148‑1 (MQ=255)
                                                                                           gCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACgg                                                          >  1:641315/1‑149 (MQ=255)
                                                                                           gCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACg                                                           <  4:1006546/148‑1 (MQ=255)
                                                                                             cGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGa                                                         >  4:456074/1‑148 (MQ=255)
                                                                                               gTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACa                                                      <  2:619682/149‑1 (MQ=255)
                                                                                                    tttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                  >  4:568407/1‑148 (MQ=255)
                                                                                                        cccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTa                                              <  4:315745/148‑1 (MQ=255)
                                                                                                         ccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                            <  2:447410/149‑1 (MQ=255)
                                                                                                          cTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                            >  3:1007409/1‑148 (MQ=255)
                                                                                                          cTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATc                                           >  2:838548/1‑149 (MQ=255)
                                                                                                           tgatgaTCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                          <  1:275457/149‑1 (MQ=255)
                                                                                                            gatgatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                          >  3:616573/1‑148 (MQ=255)
                                                                                                                  cTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtg                                                                                                                     <  3:642361/67‑1 (MQ=255)
                                                                                                                  cTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtg                                                                                                                     >  4:642361/1‑67 (MQ=255)
                                                                                                                   tAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                  >  1:979219/1‑149 (MQ=255)
                                                                                                                    aaaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                  >  4:136700/1‑148 (MQ=255)
                                                                                                                    aaaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                  >  4:647763/1‑148 (MQ=255)
                                                                                                                     aaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATgg                                 <  3:57223/148‑1 (MQ=255)
                                                                                                                       aGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGcc                               <  4:419616/148‑1 (MQ=255)
                                                                                                                        gCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGa                                    <  3:476887/142‑1 (MQ=255)
                                                                                                                        gCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGa                                    >  4:476887/1‑142 (MQ=255)
                                                                                                                                    aTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                  >  4:466712/1‑148 (MQ=255)
                                                                                                                                    aTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                  <  4:378814/148‑1 (MQ=255)
                                                                                                                                         aaaCAGGCGCCATCAGGCTCTTGTTTGAATTTCTGACGACGGTGTTTATTGCGCGTTTTCGTGATTATGGTTGCTGCTGGCGGGTAGTAAACCGTCTGCACGGAACATCGCTTTAGTCCCTCTGATGGCCTGACGAATACGGTCGCGGt            >  2:959655/1‑149 (MQ=255)
                                                                                                                                         aaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGa                            <  3:347542/133‑1 (MQ=255)
                                                                                                                                         aaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGa                            >  4:347542/1‑133 (MQ=255)
                                                                                                                                         aaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCgg             >  3:32004/1‑148 (MQ=255)
                                                                                                                                          aaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtt           <  1:163218/149‑1 (MQ=255)
                                                                                                                                           aCAGGAGCCATCAGGCTCGTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtt           <  3:851610/148‑1 (MQ=255)
                                                                                                                                            cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt         <  2:344856/149‑1 (MQ=255)
                                                                                                                                            cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt         <  1:478096/149‑1 (MQ=255)
                                                                                                                                               taGCCATCAGGCTCCTGTTTGAATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGGTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCaa      <  2:641315/148‑1 (MQ=255)
                                                                                                                                                    aTCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATAAg  >  4:205538/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAG  >  NC_000913/2496875‑2497170

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