Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,301 G→A S117S (TCC→TCT cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,3010GA97.0% 204.1 / ‑4.3 69S117S (TCC→TCTcyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (37/28);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGC  >  NC_000913/3256057‑3256546
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACAATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGACGCGGATTTCATCTCAAGGATCCTGGCTCTTAACCGAGACt                                                                                                                                                                                                                                                 >  2:452598/1‑251 (MQ=255)
 aaaCGGGACTTCATTGACGAACTTCAGGATCGCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACtt                                                                                                                                                                                                                                                <  1:176497/251‑1 (MQ=255)
     gggACTTCATTGACGAAGTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGACGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGATCTTAACAGAGCCTTTccc                                                                                                                                                                                                                                            >  2:103939/1‑251 (MQ=255)
          ttctttGACGAACTCCAGGAACTCAGCCAGCACCATTCGGGAAAGCCCCGCCAGCGGTGAGTCTTCCCCCCCCTCTGCCACACACGCCTGCTGGGTAAACCCCACACCACCGTGCGTCCCGATATGCACAATGTTGTTATGCCCGCCGACGATGGTGACACGCGCCCCCTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCccgcc                                                                                                                                                                                                                                        <  1:673908/246‑1 (MQ=255)
                 aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCcagc                                                                                                                                                                                                                                  >  1:625655/1‑249 (MQ=255)
                            gatcccagcaggcctctttCTGGAAAGCACCGTCAGCGGAGACTCTTGCTGGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGGGCGTCTCGATATCCAAAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGGGAGAAGAGGATTTCATCGCAAGGTTCCTGTATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCAt                                                                                                                                                                                                                     <  1:236718/234‑1 (MQ=255)
                              tctcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACCCACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGCTCTTAACAGAGACTTTCCCCgccgcccagatcggacgagcgc                                                                                                                                                                                                                     >  2:166869/1‑233 (MQ=255)
                              cctcAGCCAGCGGCGTTCGGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGGTCCTGGATCTTAACAGAGACTTTCCCCGCCGCca                                                                                                                                                                                                                                    <  1:166869/233‑1 (MQ=255)
                                 cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGc                                                                                                                                                                                                                 >  1:129990/1‑250 (MQ=255)
                                 cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGc                                                                                                                                                                                                                 >  1:130218/1‑250 (MQ=255)
                                  agccagcGTCTTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACCCACGCCTCCTGGGTAAACACCCCCCCACCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                                                                                                                                                                                               <  1:103050/251‑1 (MQ=255)
                                     cagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCccgccg                                                                                                                                                                                                                                       >  1:240190/1‑224 (MQ=255)
                                         gtcgtTCTGGAAAGCACCGCCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                                                                                                                                                                                        <  1:290930/251‑1 (MQ=255)
                                               cgggaagcccccggcggCGGAGACTCTTGTCCCCCCTCGCCCACACCCGCCTGCTGGGTAACCCCCCCACCATCGTCCTTCCCGATAGCCAAAATGTGGGTATGCCCGCCGACGAGGGGGACACCCGCCCACTTCGCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTGCACCGCAAGGTTCCTGGATTTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTg                                                                                                                                                                                                   <  2:625655/234‑1 (MQ=255)
                                                 ggAAAGCCCCGTCAGCGGAGCCTCTTGCTCCCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCACCGGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGtt                                                                                                                                                                                                 <  1:309475/250‑1 (MQ=255)
                                                   aaaGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCg                                                                                                                                                                                              >  1:159306/1‑251 (MQ=255)
                                                      gCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                                           <  1:170835/251‑1 (MQ=255)
                                                      gCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACAACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGAATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                                           >  1:425339/1‑251 (MQ=255)
                                                      gCACCGTCAGCGGAGACTCTTGCCCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                                           <  1:314826/251‑1 (MQ=255)
                                                        cccgtcaccggaaccccttcaccccccTCCGCCACCCCCGCCTGCGGGGTAAACCCCCCCCCATCGTGCTCCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACATGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCttt                                                                                                                                                                                         <  2:129990/227‑1 (MQ=255)
                                                        cccgcaagcggagacccttcaccGCCCTCTGCCACCCACGCCTGCTGGTGAAACACCACACCATCGTGCGTATCGATATGCACAATGTGGGTATGCCCGCAGACGATGGTGACACACGCCCACTTATCACGGTGCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCttt                                                                                                                                                                                         <  1:589397/229‑1 (MQ=255)
                                                        aCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCttt                                                                                                                                                                                         >  1:410646/1‑251 (MQ=255)
                                                          cGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                                                                                                                                                                                        >  1:322533/1‑250 (MQ=255)
                                                          cGTCAGCGGAGACTCTTGCTCGCCCTCGGCCACACACGCCTGCTGGGTAAACCCCCCCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACCTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCa                                                                                                                                                                                       <  1:591628/251‑1 (MQ=255)
                                                               gCGGAGCCTCTTGCTCGCCCTCTGCCACCCACGCCTGCTGGTTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACt                                                                                                                                                                                  <  1:560961/251‑1 (MQ=255)
                                                                    gaCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                                                                                                                                                                               >  1:152328/1‑249 (MQ=255)
                                                                     aCTCTTGCACGCCCTCTGCCACACACGCCTGCTGGGGAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCCCTACTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGCTCTTAACCGAGCATTTCCCCGCCGCCAGCACTGCTTTGTCATCGGCAATTTCCTGCGCTGTTGCGTCTTTCAACACTTCCaac                                                                                                                                                                            >  2:220527/1‑249 (MQ=255)
                                                                       tctTGCTCGCCCTCTGCCACACACGCCTGCTAGGTAAACACCACACCATCGTGCGTATCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGAGTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc                                                                                                                                                                          >  1:591364/1‑251 (MQ=255)
                                                                         ttGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCgg                                                                                                                                                                        >  1:464205/1‑251 (MQ=255)
                                                                           gCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCACACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGTATCTTAACAGAGACTTTCCCCGCCGCCCGCAGTGCTTTGGCATCGGCACTTGCCTACGCTGTTGCGTCTTTCCGCACTTCCAGCCCgccc                                                                                                                                                                      >  2:485976/1‑248 (MQ=255)
                                                                            cgcccccTCTCCCACACACGCCTGCTCGGTAAACCCCACCCCATCGTGCGTCTCGATATGCACAATGTTGTTATGCCCGCTGACGATGGTGACACACGCCCTCTTCTCACCGTCCCAGACTTTGGCGCGTGAAAAGAGGATTTCATCGTAAGGTTCCTGATGTTTACCATAGCCTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg                                                                                                                                                                      <  2:454391/246‑1 (MQ=255)
                                                                              cGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTTCGTCTTTCAGCACTTCCACCCCGGCGTTg                                                                                                                                                                   >  2:51271/1‑251 (MQ=255)
                                                                               gCCCTCTGCCACACACGCATGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCCTCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCCGCACTTCCAGCCCGGCGTTgg                                                                                                                                                                  >  2:438926/1‑251 (MQ=255)
                                                                               cccccttgccccccacGCCTGCGGTGAAAACACCACCCCATCGTGGCTCTCGATATTCACAATGTGGGTATCCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCTTCTTTCAGCACTTCCAGCCCGGCGTTgg                                                                                                                                                                  <  2:159306/238‑1 (MQ=255)
                                                                                       ccACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTccc                                                                                                                                                          >  1:439178/1‑250 (MQ=255)
                                                                                        cacacaCGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGt                                                                                                                                                                                                  >  2:440716/1‑210 (MQ=255)
                                                                                        cacacaCGCCTGCTGGGTAAACACCACACCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGt                                                                                                                                                                                                  <  1:440716/210‑1 (MQ=255)
                                                                                        cacacaCGCCTGCTGGGGAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTCATCGCAAGGTTCCTCAATCTTAACAGAGACTTTCCCCGCCGCCAGCAGAGCTTGGGCATCGGCACATGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGATGGGACTTCCAc                                                                                                                                                         >  2:323702/1‑251 (MQ=255)
                                                                                             ccgcccgcggggaaacccccacaccaGCGGGCGTCGGGACATGCACAGGGTTGGTGTGCCCTCCGTGGATGTGGACACCCGCCCACTTCTCACCGTTCCAGTCTGTGGCGCGGTAGAAGAGGATTTCATCGCAGGGTTCCTGGACTTCAACGGAGCCTTCCCCCGCCTCCAGCAGTGCTTTGGTATCGGCAATTGCCTGCGCTGTGGCGTCCTGCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAAc                                                                                                                                                    <  2:439178/233‑1 (MQ=255)
                                                                                              cGCCTGCTGGGTAAACGCCACACCACCGTGCGTTTCGATATGCACAATGTTGGTTTGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACg                                                                                                                                                   <  1:602931/251‑1 (MQ=255)
                                                                                               gCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGc                                                                                                                                                  >  1:587752/1‑251 (MQ=255)
                                                                                                     tgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGccccc                                                                                                                                              >  1:311297/1‑249 (MQ=255)
                                                                                                     tgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACATAACGccccc                                                                                                                                              >  1:352689/1‑249 (MQ=255)
                                                                                                       ggTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATAGCAAGGGTCCTGGATCTTAACAGAGACTTTCCCCTCCGCCAGCAGCGCTTTGGCATAGGCCATTGCCTGCGCTGTTGCGTCTTTCAGCACTGCCAGCCCGGCGTTGGCATTTACACATAACGCCCCCaacg                                                                                                                                          >  2:465312/1‑248 (MQ=255)
                                                                                                                        tgggcggctgtatatGCACATGGTGGGTATGCCCGCCGACGAGGGTGACCCCCCCCCCCTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAg                                                                                                                          <  1:502292/239‑1 (MQ=255)
                                                                                                                        cGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGCGAAGAGGCTTTCATCGCACGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCCTTTCAACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGc                                                                                                                         >  2:340750/1‑251 (MQ=255)
                                                                                                                               ctcGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTTCGCTGTTGCGTCTTTCAGCACTTCCATCCCGGCGTTGGCATTTCCACCCAACGCCCCCATCGCCGCCGCAATCGACAGCCccacc                                                                                                                   >  2:152907/1‑250 (MQ=255)
                                                                                                                                       ccaataTGTTGGTATGCCCGCCGACGATGTTGACCCACGCCCACTTCCCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATTTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTg                                                                                                          <  1:323702/246‑1 (MQ=255)
                                                                                                                                              gTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGa                                                                                                    >  1:227256/1‑250 (MQ=255)
                                                                                                                                                gggTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACg                                                                                                 <  1:67591/250‑1 (MQ=255)
                                                                                                                                                    aTGCCCGCCGACGATGGTGACACCCGCCCACTTCTCACCGTCCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGa                                                                                             <  1:636416/251‑1 (MQ=255)
                                                                                                                                                    aTGCCCGCCGACGATGGCGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCACCGCCGCCAGCAGTGCCTTGGCATTTGTAATTGCTTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCTCCACCATTCCCGTGCCGGGAACGGTg                                                                                              >  1:375259/1‑250 (MQ=255)
                                                                                                                                                             gacgaTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCACGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGAGTTGGCATTTCCACCTACCGCCCCCAGCGCCGCCGCAATCGGCAGCCCCCCCATTCCCGTGCCGGGAACGGTAACCCCCAGAc                                                                                    >  2:618783/1‑251 (MQ=255)
                                                                                                                                                                 aTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCCCCATTCCCGTGCCGGGAAAGGTGACGCCCAgcacgtt                                                                                >  2:122684/1‑245 (MQ=255)
                                                                                                                                                                   ggTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGCCCGtttt                                                                              >  2:242420/1‑249 (MQ=255)
                                                                                                                                                                   ggTGACACACGCCCACTCCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGATTTAACCGCAAGTTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttct                                                                              <  1:536703/251‑1 (MQ=255)
                                                                                                                                                                       acacacGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcat                                                                          >  1:431286/1‑251 (MQ=255)
                                                                                                                                                                        cacccGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                         <  1:438926/247‑1 (MQ=255)
                                                                                                                                                                                    tctcACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTcc                                                                                                              >  2:137295/1‑202 (MQ=255)
                                                                                                                                                                                    tctcACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                                              >  1:512740/1‑250 (MQ=255)
                                                                                                                                                                                    tctcACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGACCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTcc                                                                                                              <  1:137295/202‑1 (MQ=255)
                                                                                                                                                                                           gTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAgg                                                      >  1:607707/1‑251 (MQ=255)
                                                                                                                                                                                            ttCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCACCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGc                                                     >  1:21556/1‑251 (MQ=255)
                                                                                                                                                                                                 aaCTTTGGCACGTGAAACGAGGAATTCGTCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGACGCCAGCAGTGCTTTGGCCTCGGAAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCCCCGACGCACTCGGCCGCCACACCATTCCCGTGCCGGGAACGGGGACGCACAGACAGATCTTCATACGATTTGGCGAAACCCacgcctcg                                                 >  2:660790/2‑243 (MQ=255)
                                                                                                                                                                                                            gggagaAGGGAATTTCAACGCAATGTCCCTGGACCTTAACAGATACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAAcc                                     <  1:70522/249‑1 (MQ=255)
                                                                                                                                                                                                                gaagaGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGCCGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCCGGCTTCTACACGTTCAACCGGAc                                 >  1:352868/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc                   <  1:311422/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGATCTAAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAcc          <  1:441820/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGATGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCCGCCCCACCATTCCCGTGCCCGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACACATGCTTCTAAACGTTAAACCGGACCTCACAGTTCTGCCACCGCAAAcgccgccgc  >  2:454303/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCGCCGCCGC  >  NC_000913/3256057‑3256546

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