Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,357,064 C→A 5.3% D519Y (GAC→TAC)  cadA ← lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,357,0640CA 5.3% 218.1 / 2.7 76D519Y (GAC→TAC) cadAlysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base A (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TTTTTCACACGCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGTGCCAGGTGC  >  NC_000913/4356916‑4357207
                                                                                                                                                    |                                                                                                                                               
tttttCACACGCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTtcgtc                                                                                                                                                 <  1:331018/149‑1 (MQ=255)
   ttCACACGCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAg                                                                                                                                              <  2:527618/149‑1 (MQ=255)
         cGCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCGGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTAttt                                                                                                                                        <  2:638422/149‑1 (MQ=255)
          gCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTc                                                                                                                                       >  1:106032/1‑149 (MQ=255)
          gCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGGCGAggcacttg                                                                                                                                       >  2:560808/1‑143 (MQ=255)
             ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCTCGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAATAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTGTTTCGcc                                                                                                                                    >  1:486640/1‑149 (MQ=255)
              gTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCa                                                                                                                                   <  1:593673/149‑1 (MQ=255)
                tCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCGCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACg                                                                                                                                 <  1:703940/149‑1 (MQ=255)
                 cAGGTCGAACGCACGTTTAAAGTCAGTCAGGGCACGCAGCAGGCTCAGTGCTTTGGTCTTACAAGTACCGTTGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGa                                                                                                                                <  1:360534/149‑1 (MQ=255)
                   ggTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATg                                                                                                                              <  1:653646/149‑1 (MQ=255)
                    gTCGAACGCCCGGTTAAGGTCAGTCAGAGCACGCCGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGGTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGc                                                                                                                             <  2:138459/149‑1 (MQ=255)
                    gTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGc                                                                                                                             >  1:594625/1‑149 (MQ=255)
                     tCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGTCTCAGTGCTTTGGTTTTATCGATACCGATTCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTTGTCGAGGTATTTTGCCTCGATGCt                                                                                                                            >  1:210807/1‑149 (MQ=255)
                     tCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTCTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGGTTTCTCAACAACGATGCCATGTTCGTCGAGGGATTTCCCCCCGATGCt                                                                                                                            >  2:590390/1‑149 (MQ=255)
                        aacgaacgTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCCGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGc                                                                                                                         <  1:27666/149‑1 (MQ=255)
                         aCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGcc                                                                                                                        <  1:71131/149‑1 (MQ=255)
                           gcccGTTTAAAGTCCGCCAGCTCACGCCGCAGGCTTAGTGCTTTGGTCTTATCGACACCGATGCTGAACACGAACAGCAGGTTATACGGCCCGGGTTTCTCAAAAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCgg                                                                                                                      <  1:213270/146‑1 (MQ=255)
                             aCGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCAGGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGaa                                                                                                                    >  2:710555/1‑149 (MQ=255)
                               gTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGACTa                                                                                                                  <  2:131950/149‑1 (MQ=255)
                                tttAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATAc                                                                                                                 >  2:737836/1‑149 (MQ=255)
                                  tAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCa                                                                                                               <  2:127329/149‑1 (MQ=255)
                                    aaGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGTTGCTGAATAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGAAGCCCTGTTTGTGGAGGGATGTCGCCCCGATTATGGCCGGAATAcaaca                                                                                                             >  2:377581/1‑145 (MQ=255)
                                    aaGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCaaa                                                                                                             <  1:523471/149‑1 (MQ=255)
                                          gtcagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAATGGCt                                                                                                       >  2:448618/1‑149 (MQ=255)
                                           tcagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTAGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTc                                                                                                      >  2:83545/1‑149 (MQ=255)
                                           acagAGCACGCAGCAGGCTCAATGATTTGGTCTTATAGATACCGATGCTGAACAGGAACCACAGGCAATACGATATGGTTTTGTTAACAACTAAGACATGTTCGAAGAGGTATTTCGCCACGCTGCTGGCCGGAATACCCAAGACGCTc                                                                                                      >  2:232076/2‑149 (MQ=255)
                                            cagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCa                                                                                                     >  2:732948/1‑149 (MQ=255)
                                              gagCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATg                                                                                                   >  1:715959/1‑149 (MQ=255)
                                                 cACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTg                                                                                                >  1:43171/1‑149 (MQ=255)
                                                   cGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCGCGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGcc                                                                                              <  1:232105/149‑1 (MQ=255)
                                                    gcagcaGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCCGAAGACCAAAGTCACCCAAGGTGCCg                                                                                             >  2:368073/1‑149 (MQ=255)
                                                            cTCATTGCTTTGGTCTTAGCGATACCGATGCTGAACAGCAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTc                                                                                     <  2:22216/149‑1 (MQ=255)
                                                            cTCAGTGCTTCGTTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTc                                                                                     <  2:123356/149‑1 (MQ=255)
                                                            cTCAGGGCTTTGGTCTTAGCGATCCCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTc                                                                                     <  2:682808/149‑1 (MQ=255)
                                                              cAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAAGTATTTCGCCACGATGCTGGCCGGAATACCAAAATCGCTCATGGTGCCGTCTTTTTCCa                                                                                   >  1:334569/1‑149 (MQ=255)
                                                                   cTTTGGTCGTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATcccc                                                                              <  1:293454/149‑1 (MQ=255)
                                                                    tttGGTCTTATCGGTACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATTCCAAAGTTGCTCATGGTGCCGTCTTTTTCCATCCCCg                                                                             <  2:725825/149‑1 (MQ=255)
                                                                     ttGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCgg                                                                            <  2:174529/149‑1 (MQ=255)
                                                                      tGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGa                                                                           <  2:387695/149‑1 (MQ=255)
                                                                       ggTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTAACCGGCAATTTTTTGTCAACAGCGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGGATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAg                                                                          >  2:311533/1‑149 (MQ=255)
                                                                           ttATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGGATACCAAAATCGCTCATGGTGCCGTCTTGTTACATACCCGGAGCcag                                                                      >  2:245139/1‑149 (MQ=255)
                                                                               cGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAgg                                                                  <  1:737836/149‑1 (MQ=255)
                                                                                   cccGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGCGTCAGCAGGGTGa                                                              <  2:251964/148‑1 (MQ=255)
                                                                                    ccGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGAc                                                             >  1:121151/1‑149 (MQ=255)
                                                                                     cGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTTGCCACGATGCTGGCCGGAATACCAATGTCGCTAAAGGTACTCTTGGTTTCCATCCCAGGTCTCAGCAGGGTGACt                                                            <  1:554498/149‑1 (MQ=255)
                                                                                      gATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACtt                                                           <  2:362818/149‑1 (MQ=255)
                                                                                         gCTGAACAGGAACAGCAGGTTATACGGACCTGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGa                                                        <  1:218555/149‑1 (MQ=255)
                                                                                          cTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTAGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGAt                                                       <  1:83545/149‑1 (MQ=255)
                                                                                          cTGAACAGGAACAGCAGGTAATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTAGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGAt                                                       <  1:232076/149‑1 (MQ=255)
                                                                                           tGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATc                                                      >  2:43882/1‑149 (MQ=255)
                                                                                           tGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCACGGAGTCAGCAGGGTGGCTTTGGTc                                                      >  2:24047/1‑149 (MQ=255)
                                                                                             aacgggaacagCAGGTTATACGGACCGGTTTGCTCAACAACGATGCCATGCTCGTCGAGGTATTTCGCCACGATGCCGGCCGGAATACTAAAGTCGCTCATGGTGCAGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCgg                                                    <  2:604891/149‑1 (MQ=255)
                                                                                             aacaggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACTAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCgg                                                    >  2:645070/1‑149 (MQ=255)
                                                                                              acaggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGCCTTTGATCggg                                                   >  1:211678/1‑149 (MQ=255)
                                                                                              acaggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCggg                                                   >  1:228928/1‑149 (MQ=255)
                                                                                              acaggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCggg                                                   <  1:483370/149‑1 (MQ=255)
                                                                                               caggaacggCTGGTTTTACGGCCCGGTTTTCTCAACAACGAGGCCATGATCGTCGCGGTATTTCGCCACGATGCGGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGt                                                  <  2:268351/149‑1 (MQ=255)
                                                                                                    acagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGa                                             <  1:467818/149‑1 (MQ=255)
                                                                                                        cagGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACa                                         <  1:188716/149‑1 (MQ=255)
                                                                                                             tataCGCACCGGTGTTCTCAACAACGCTGCCAGGTTCGTCGAGGGATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGc                                    <  1:241470/149‑1 (MQ=255)
                                                                                                                 cGGACCGGTTTTCCCAACAACGATGCCAGGTTCGTCGAGGTATTTCGCCGAGGTGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGt                                <  1:21639/149‑1 (MQ=255)
                                                                                                                   gACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTa                              <  2:654210/149‑1 (MQ=255)
                                                                                                                        gTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATCCATGTGCTCGTTAtaggt                         >  2:751063/1‑145 (MQ=255)
                                                                                                                        gTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGAt                         <  1:732948/149‑1 (MQ=255)
                                                                                                                        gTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGAt                         >  1:537177/1‑149 (MQ=255)
                                                                                                                         ttttCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCGCGTTCTCGCTg                        >  2:372156/1‑149 (MQ=255)
                                                                                                                             ctcAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCACGATACATGTGCTCGTTATCGATGtttt                    >  1:467413/1‑149 (MQ=255)
                                                                                                                              tcAACACCGATGCCATGCTCGTCGAGGTATTTCGCCATGCTGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGttttt                   <  1:606435/149‑1 (MQ=255)
                                                                                                                               caacaaCGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTg                  >  1:371350/1‑149 (MQ=255)
                                                                                                                                aacaacGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGTTCGGGTCAGGATACATGTGCTCGTTATCGATGTTTTTGa                 >  1:445632/1‑149 (MQ=255)
                                                                                                                                 acaacGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCGGCAGGGTGACTGTAATCGGGTCAAGATACATGTGATCGTTATCGATGTTTTTTaa                >  1:288041/1‑149 (MQ=255)
                                                                                                                                  caacGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAg               <  2:106032/149‑1 (MQ=255)
                                                                                                                                      gATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGt           <  1:674426/149‑1 (MQ=255)
                                                                                                                                      gATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGt           >  1:714223/1‑149 (MQ=255)
                                                                                                                                             tGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCCGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGTGCCAGGt    <  2:416029/149‑1 (MQ=255)
                                                                                                                                              gTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGTGCCAGGTg   >  2:201629/1‑149 (MQ=255)
                                                                                                                                               ttCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGTGCCAGGTGc  >  1:754553/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTTTTCACACGCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATGGTGCCGTCTTTTTCCATCCCCGGAGTCAGCAGGGTGACTTTGATCGGGTCAAGATACATGTGCTCGTTATCGATGTTTTTGAAGCCGTGCCAGGTGC  >  NC_000913/4356916‑4357207

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