Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,502 T→C 14.2% A53A (GCA→GCG)  tdcD ← propionate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,5020TC14.2% 185.6 / 20.4 77A53A (GCA→GCG) tdcDpropionate kinase
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (6/5);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

TCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGA  >  NC_000913/3263356‑3263649
                                                                                                                                                   |                                                                                                                                                   
tCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCt                                                                                                                                                  <  2:1470828/150‑1 (MQ=255)
tCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCt                                                                                                                                                  <  1:417076/150‑1 (MQ=255)
   gTAATAATGGCGGACTCGGTAAAACTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                               <  2:624576/150‑1 (MQ=255)
       taatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctcac                                                                                                                                           <  1:1463128/150‑3 (MQ=255)
       taatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctcac                                                                                                                                           <  1:1470194/150‑3 (MQ=255)
       taatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCGTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctcac                                                                                                                                           <  1:515103/150‑3 (MQ=255)
        aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcc                                                                                                                                          >  2:182082/1‑150 (MQ=255)
         atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTccc                                                                                                                                         <  2:471526/150‑1 (MQ=255)
          tGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcccc                                                                                                                                        <  1:281688/150‑1 (MQ=255)
            gCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTAAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGCTCAccacag                                                                                                                                      <  2:297645/150‑5 (MQ=255)
             cGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAt                                                                                                                                     >  2:1373865/1‑150 (MQ=255)
                   cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                               >  2:952543/1‑150 (MQ=255)
                      tAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATa                                                                                                                            <  2:362013/150‑1 (MQ=255)
                        aaatTATTGCCGTCTTGATCGATGCTGTGGCCATTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAg                                                                                                                          <  1:705782/150‑1 (MQ=255)
                                 ccgccgTGAGCGATGCGGTGTCCAATTAAGTCCACAAAGTCATTTAAATTCCGTTTTTCCACTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAttt                                                                                                                 >  2:565150/1‑150 (MQ=255)
                                 ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAttt                                                                                                                 >  2:1384500/1‑150 (MQ=255)
                                       tGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAg                                                                                                           <  2:243438/150‑1 (MQ=255)
                                       tGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAg                                                                                                           <  2:239775/150‑1 (MQ=255)
                                                    ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCgg                                                                                              <  2:355759/150‑1 (MQ=255)
                                                    ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAAGGATAAGAATGCCTTTTACGAGTTAATACCGTCgg                                                                                              >  2:1275808/1‑150 (MQ=255)
                                                     gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                             <  2:594975/150‑1 (MQ=255)
                                                     gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                             >  2:90129/1‑150 (MQ=255)
                                                      ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAAGGATAAGAATGCATTTTACGAGTTAATAACGTCGGCa                                                                                            >  2:372496/1‑150 (MQ=255)
                                                             aGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTg                                                                                     >  2:274104/1‑150 (MQ=255)
                                                              ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                    >  2:291476/1‑150 (MQ=255)
                                                              ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                    >  2:471698/1‑150 (MQ=255)
                                                               gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                   <  1:182082/150‑1 (MQ=255)
                                                                ccACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                   <  1:575453/150‑1 (MQ=255)
                                                                  acacTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACAtt                                                                                <  2:726998/150‑1 (MQ=255)
                                                                     cTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTTGCATTTAt                                                                             >  1:1227432/1‑150 (MQ=255)
                                                                         cATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                         >  1:534317/1‑150 (MQ=255)
                                                                            ttAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                      <  1:1468361/150‑1 (MQ=255)
                                                                                    cGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                              <  2:527476/150‑1 (MQ=255)
                                                                                    cGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                              >  1:530855/1‑150 (MQ=255)
                                                                                       ttttCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCCGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGGTAATACCGTCGGCAATACCTGACATTAATACTTGACAGTCGCTGGcc                                                           >  1:1495843/1‑149 (MQ=255)
                                                                                        tttCCAGTTC‑AAATGCAATTGCCCTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGACAt                                                          <  2:191112/150‑1 (MQ=255)
                                                                                          tgcAGTTC‑AAATGCAATTGCCTTCACTGCACCTTCGGAGCTGTGGTGCGCCCCCGGCGGTGGCTGTCCCGCGTTTGAGGGTAAGAATACATTTTGCGAGGTGATGCCCTAGGCAAAACATGACCTGAATAGTTGCCACTAGCTTGCAtag                                                        >  2:1266663/3‑148 (MQ=255)
                                                                                          tCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                        >  2:1017818/1‑150 (MQ=255)
                                                                                          tCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCg                                                        <  2:49014/150‑1 (MQ=255)
                                                                                            cAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                      >  2:1195063/1‑150 (MQ=255)
                                                                                             aGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGc                                                     <  1:1601906/150‑1 (MQ=255)
                                                                                              gTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCa                                                    <  1:1264699/150‑1 (MQ=255)
                                                                                                 c‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACgg                                                 <  2:404723/150‑1 (MQ=255)
                                                                                                 c‑AAATGCAATTGCCTTCAATGCAACTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACgg                                                 >  2:1265649/1‑150 (MQ=255)
                                                                                                    aaTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTTCAGGCTATACCCCAATTTAGGATAATAATTCATTTTACGAGTTTATAACGTAGGCAAAAACAGAAAATAATAATTGAAAGTAGCTGGCATAGAGCACCGaa                                               >  2:1467718/1‑150 (MQ=255)
                                                                                                             gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                      >  2:1128549/1‑150 (MQ=255)
                                                                                                              ccTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCg                                     <  1:609154/150‑1 (MQ=255)
                                                                                                              ccTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCg                                     <  2:250894/150‑1 (MQ=255)
                                                                                                               cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                                    <  1:1371606/150‑1 (MQ=255)
                                                                                                               cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                                    >  1:65946/1‑150 (MQ=255)
                                                                                                                  cAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTAACAGTAGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                                 >  2:1199405/1‑150 (MQ=255)
                                                                                                                   aaTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGAc                                >  1:499786/1‑150 (MQ=255)
                                                                                                                   aaTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTACCCCAATTAAGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATAACTGACAATAATACTTCACAGTCGCTGGCATCGAGCACGGAAAAAATAATCGAAGAc                                >  2:858291/1‑150 (MQ=255)
                                                                                                                    aTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                               <  2:499786/150‑1 (MQ=255)
                                                                                                                       cACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAc                            >  1:1321903/1‑150 (MQ=255)
                                                                                                                        aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                           <  1:1467718/150‑1 (MQ=255)
                                                                                                                         ccTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACAATAATAATTCACAGTCGCTGGCATAGAGCAAGGAAAACTTAATAGAAGACGAACca                          >  2:1177561/1‑150 (MQ=255)
                                                                                                                              gTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGtt                     <  1:89965/150‑1 (MQ=255)
                                                                                                                               tAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                    <  1:1384500/150‑1 (MQ=255)
                                                                                                                                     tggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc              <  1:90129/150‑1 (MQ=255)
                                                                                                                                     gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc              >  1:994607/5‑150 (MQ=255)
                                                                                                                                     gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc               <  1:321867/145‑1 (MQ=255)
                                                                                                                                     gggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc               >  2:321867/5‑149 (MQ=255)
                                                                                                                                      ggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa             >  2:630573/1‑150 (MQ=255)
                                                                                                                                      ggagAGCCAGCACCGCTGGCTCACCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa             >  1:552445/4‑150 (MQ=255)
                                                                                                                                       gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa            >  2:1225276/3‑150 (MQ=255)
                                                                                                                                       gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa            >  2:302675/3‑150 (MQ=255)
                                                                                                                                        tgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa           <  1:171201/150‑1 (MQ=255)
                                                                                                                                          aGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaac         <  1:565150/150‑1 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca        >  1:1124474/1‑150 (MQ=255)
                                                                                                                                           gCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca        >  1:461236/1‑150 (MQ=255)
                                                                                                                                               gCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg    <  1:291476/150‑1 (MQ=255)
                                                                                                                                                cGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAACTTACATTAATACTTCACAGTCGCTGTCATCGAGCACGGAAAACTTAATCGAAGACTAACCACAGTTAATAACCAAAACAACCgg   <  1:162133/150‑1 (MQ=255)
                                                                                                                                                 ggTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGa  >  2:13290/1‑150 (MQ=255)
                                                                                                                                                 ggTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATAATTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGa  >  1:684323/1‑150 (MQ=255)
                                                                                                                                                  gTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACtt                                          <  1:1334565/109‑1 (MQ=255)
                                                                                                                                                  gTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACtt                                          >  2:1334565/1‑109 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTC‑AAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGA  >  NC_000913/3263356‑3263649

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