Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 127,442 T→A 7.3% V583E (GTG→GAG)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913127,4420TA 7.3% 216.9 / 3.9 74V583E (GTG→GAG) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base A (4/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAA  >  NC_000913/127302‑127587
                                                                                                                                            |                                                                                                                                                 
tAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATg                                                                                                                                           <  3:13233/149‑1 (MQ=255)
  aGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATgg                                                                                                                                          <  4:84611/148‑1 (MQ=255)
       aCTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTaaa                                                                                                                                      <  8:419040/147‑1 (MQ=255)
           cGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGtt                                                                                                                                 <  5:164538/148‑1 (MQ=255)
            gggCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCg                                                                                                                               >  5:380254/1‑149 (MQ=255)
            gggCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCg                                                                                                                               >  2:273627/1‑149 (MQ=255)
             ggCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGTGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGTAATGGTAAAGAGTTCGt                                                                                                                              <  5:7131/149‑1 (MQ=255)
                gAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGAGTGGAATGGTAAAGAGTTCGTGcc                                                                                                                           >  5:394989/1‑149 (MQ=255)
                     gCAGGGCGGTTGCTTCACCATCTCCATCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGt                                                                                                                       <  4:136676/148‑1 (MQ=255)
                       aGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTg                                                                                                                    <  1:273627/149‑1 (MQ=255)
                        tggCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACACACTTCGAGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGa                                                                                                                   <  4:282681/148‑1 (MQ=255)
                         ggCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGAt                                                                                                                  >  4:388832/1‑149 (MQ=255)
                            ggTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGAtgct                                                                                                               >  7:142372/1‑149 (MQ=255)
                                 cTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGa                                                                                                          <  1:291105/149‑1 (MQ=255)
                                  ttCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGAt                                                                                                         >  1:401846/1‑149 (MQ=255)
                                  ttCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGAt                                                                                                         >  7:59431/1‑149 (MQ=255)
                                    cACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGAt                                                                                                         <  3:388832/147‑1 (MQ=255)
                                     aCCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGAttt                                                                                                       <  8:142372/148‑1 (MQ=255)
                                     aCCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTtc                                                                                                      >  3:171177/1‑149 (MQ=255)
                                      ccATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGAGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATCCTGACGAGTtct                                                                                                     >  8:166654/1‑149 (MQ=255)
                                       cATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTtctc                                                                                                    >  7:96271/1‑149 (MQ=255)
                                       cATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGAGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTtctc                                                                                                    >  7:223552/1‑149 (MQ=255)
                                        aTCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTtctct                                                                                                   <  6:380254/149‑1 (MQ=255)
                                         tctcCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTtctctc                                                                                                  >  5:228284/1‑149 (MQ=255)
                                             cAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTcc                                                                                               <  2:16293/148‑1 (MQ=255)
                                               gCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTc                                                                                            <  4:171177/149‑1 (MQ=255)
                                                     gcggcCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGaccac                                                                                      >  2:53846/1‑149 (MQ=255)
                                                       ggcCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTATCTCTCCTTCGACCACcg                                                                                    >  7:293904/1‑149 (MQ=255)
                                                          cTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGt                                                                                 >  5:403938/1‑149 (MQ=255)
                                                          cTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGAGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGt                                                                                 <  7:312536/149‑1 (MQ=255)
                                                           tGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGCGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGCTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGt                                                                                 >  5:131290/1‑148 (MQ=255)
                                                               tactacCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCAAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGAGTTCTCTCTCCTTCGACCACCGCGTGATCg                                                                            >  3:26332/1‑149 (MQ=255)
                                                                 ctacCCACTTCGCGCCGATTGTGAACGCGCCGTAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGTTAAAGCGTTAGTGCCGAGTCTGATGATGCCGATTTATCTCTCCTTCGACCACCGCGGGATCgcc                                                                          >  6:125181/1‑147 (MQ=255)
                                                                    cccACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACgg                                                                        <  6:228284/148‑1 (MQ=255)
                                                                    cccACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACgg                                                                        <  8:96271/148‑1 (MQ=255)
                                                                    cccACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACg                                                                         <  7:160545/147‑1 (MQ=255)
                                                                    cccACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGAGGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGGCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACg                                                                         <  4:298165/147‑1 (MQ=255)
                                                                      cACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGc                                                                     <  2:319907/149‑1 (MQ=255)
                                                                       aCTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCt                                                                    >  7:122479/1‑149 (MQ=255)
                                                                            gcgcCGATTGTGAACGCGCCGGAAGTGGATATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGATTTCTTGACGCGTCTGATGCTGCCGATTTCTGTCTCCTTCGACCACCGCGTGATCGACGGGGCTGATg                                                                <  6:402374/148‑1 (MQ=255)
                                                                                cGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGcc                                                           >  7:420621/1‑149 (MQ=255)
                                                                                  aTTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGAGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCg                                                         >  4:341221/1‑149 (MQ=255)
                                                                                  aTTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGAGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTATCTCCTTCGACCCCCGCGTTATCGACGGTGCTGATGGTGCACg                                                         >  4:346797/1‑149 (MQ=255)
                                                                                   ttGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGCTTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGt                                                        >  6:283582/1‑149 (MQ=255)
                                                                                      tgAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGttt                                                      <  8:59431/148‑1 (MQ=255)
                                                                                         aCGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCAt                                                   >  3:346315/1‑148 (MQ=255)
                                                                                          cgcgCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTa                                                 >  8:326675/1‑149 (MQ=255)
                                                                                            cgcCGGAAGTGGCTATCCGCGGCGTTTCCAATTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCATCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTAcc                                               <  6:214552/149‑1 (MQ=255)
                                                                                              ccGGATGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACcat                                             <  8:59707/149‑1 (MQ=255)
                                                                                               cGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACcatc                                            <  1:329418/149‑1 (MQ=255)
                                                                                               cGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACcatc                                            <  2:401846/149‑1 (MQ=255)
                                                                                                 gAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACcatcat                                          >  8:256756/1‑149 (MQ=255)
                                                                                                    gTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATg                                                                <  7:1996/124‑1 (MQ=255)
                                                                                                    gTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATg                                                                >  8:1996/1‑124 (MQ=255)
                                                                                                    gTGGCTATCATCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTTATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTaa                                       <  5:125181/149‑1 (MQ=255)
                                                                                                     tGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGCCCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTaa                                       >  2:140134/1‑148 (MQ=255)
                                                                                                     tGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTaac                                      <  3:19122/149‑1 (MQ=255)
                                                                                                       gCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGgtg                                                                                                                                                  >  4:257887/1‑39 (MQ=255)
                                                                                                       gCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGgtg                                                                                                                                                  <  3:257887/39‑1 (MQ=255)
                                                                                                       gCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGGCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTaaca                                     <  1:410000/148‑1 (MQ=255)
                                                                                                       gCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTaacaa                                    >  8:419487/1‑149 (MQ=255)
                                                                                                           tCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCAt                                                   >  7:139334/1‑130 (MQ=255)
                                                                                                           tCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCAt                                                   <  8:139334/130‑1 (MQ=255)
                                                                                                             cTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCt                              >  5:269607/1‑149 (MQ=255)
                                                                                                              tCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTg                             >  5:280934/1‑149 (MQ=255)
                                                                                                               cGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGt                            >  1:162834/1‑149 (MQ=255)
                                                                                                                ggCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGt                            <  7:419487/148‑1 (MQ=255)
                                                                                                                ggCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTATCCTTCGACCTCCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTc                           <  7:166654/149‑1 (MQ=255)
                                                                                                                       ccAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGAGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACAtt                    >  8:62070/1‑149 (MQ=255)
                                                                                                                        cAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTc                   <  6:269607/149‑1 (MQ=255)
                                                                                                                            tCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGGTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCg               <  2:27416/149‑1 (MQ=255)
                                                                                                                                 gATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTgg          >  2:213171/1‑149 (MQ=255)
                                                                                                                                 gATGGAGCCGGAGTGTAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTgg          <  3:341221/149‑1 (MQ=255)
                                                                                                                                 gATGGAGCCGGAGTGGAATTGTAAAGAGTGCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGGGCCCGTTTCATTACCATCATTACCAACACGCGGTCTGACATTCGCCGTCTgg          <  3:346797/149‑1 (MQ=255)
                                                                                                                                  aTGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGt         <  7:341738/149‑1 (MQ=255)
                                                                                                                                         cGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGAGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATgtaa  <  7:62070/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TAAGCTGACTGCGGGCGAAATGCAGGGCGGTTGCTTCACCATCTCCAGCATCGGCGGCCTGGGTACTACCCACTTCGCGCCGATTGTGAACGCGCCGGAAGTGGCTATCCTCGGCGTTTCCAAGTCCGCGATGGAGCCGGTGTGGAATGGTAAAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAA  >  NC_000913/127302‑127587

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