Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,231 T→A 7.2% T84S (ACT→TCT)  dcyD ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,2310TA 7.2% 278.8 / 10.6 83T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (33/44);  new base A (4/2);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

AGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAT  >  NC_000913/1999086‑1999378
                                                                                                                                                      |                                                                                                                                                   
aGTGCGTCGCACATTTCAATCTGGGTAATGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTg                                                                                                                                                  >  2:453812/1‑149 (MQ=255)
 gTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGACTga                                                                                                                                                 <  3:496918/149‑2 (MQ=255)
      tCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGACTGCCGcag                                                                                                                                            >  4:701089/1‑148 (MQ=255)
        gCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACAt                                                                                                                                          >  4:819302/1‑149 (MQ=255)
        gCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGACTGCCGcagtc                                                                                                                                          <  3:701089/149‑4 (MQ=255)
         cacaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATg                                                                                                                                         >  1:628584/1‑149 (MQ=255)
            aTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                      <  1:253605/149‑1 (MQ=255)
            aTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                      >  1:801520/1‑149 (MQ=255)
                cAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCTCATGGTTAGAc                                                                                                                                  <  3:681563/149‑1 (MQ=255)
                cAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                                                                  <  4:871961/149‑1 (MQ=255)
                  aTCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTg                                                                                                                                >  2:1045480/1‑149 (MQ=255)
                    cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCACTGCGGCAGTCTGCCGCACATGGTTAGa                                                                                                                                   <  4:250061/149‑1 (MQ=255)
                     tGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                             >  4:917129/1‑149 (MQ=255)
                      gggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTa                                                                                                                                     >  3:200275/1‑140 (MQ=255)
                      gggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTa                                                                                                                                     <  4:200275/140‑1 (MQ=255)
                       ggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCg                                                                                                                           >  4:297420/1‑149 (MQ=255)
                        gTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATg                                                                                                                                         <  1:316128/134‑1 (MQ=255)
                        gTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATg                                                                                                                                         >  2:316128/1‑134 (MQ=255)
                            tGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                                                                      >  2:570677/1‑149 (MQ=255)
                             gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCgg                                                                                                                     >  2:1058587/1‑149 (MQ=255)
                               aCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCa                                                                                                                   <  3:368557/149‑1 (MQ=255)
                               aCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCa                                                                                                                   <  4:982633/149‑1 (MQ=255)
                                   aTCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTaa                                                                                                               >  1:472786/1‑149 (MQ=255)
                                      cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCa                                                                                                            <  2:5071/149‑1 (MQ=255)
                                       aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACt                                                                                                                                 >  2:343804/1‑127 (MQ=255)
                                       aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACt                                                                                                                                 <  1:343804/127‑1 (MQ=255)
                                       aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                           >  2:762607/1‑149 (MQ=255)
                                        gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                          >  4:692980/1‑149 (MQ=255)
                                            caaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTATGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTAt                                                                                                      <  2:904413/149‑1 (MQ=255)
                                               gCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCgg                                                                                                   >  4:314395/1‑149 (MQ=255)
                                                 gATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCa                                                                                                 >  4:128564/1‑149 (MQ=255)
                                                   ttGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                                                                               <  2:973889/149‑1 (MQ=255)
                                                      ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTc                                                                                            <  1:585307/149‑1 (MQ=255)
                                                        gTTGGTTAAATTGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                          >  1:556647/1‑149 (MQ=255)
                                                        gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                          >  3:241846/1‑149 (MQ=255)
                                                             ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGTGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                     <  3:106065/149‑1 (MQ=255)
                                                             ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTACTCAGCGTATCGGCACCTTCACGCaga                                                                                     <  3:652451/149‑1 (MQ=255)
                                                             ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                     >  3:632783/1‑149 (MQ=255)
                                                             ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                     <  2:892275/149‑1 (MQ=255)
                                                                aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                                                                  >  3:954142/1‑149 (MQ=255)
                                                                  tAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                                <  2:829880/149‑1 (MQ=255)
                                                                    gTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                                                                              >  3:311045/1‑149 (MQ=255)
                                                                      tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCATCGTATCGGCACCTTCACGCAGAGCATCTGcc                                                                            <  1:935690/149‑1 (MQ=255)
                                                                          tGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCga                                                                        <  1:275125/149‑1 (MQ=255)
                                                                           gCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGGCCCGGCAGTAATCAGCGTATCGGCATCTTCACGCAGAGCATCTGCCGCgag                                                                       <  1:847359/149‑1 (MQ=255)
                                                                              gTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaaa                                                                    >  3:845956/1‑149 (MQ=255)
                                                                              gTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaaa                                                                    <  4:863508/149‑1 (MQ=255)
                                                                                 gtgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                                                                 <  3:35504/149‑1 (MQ=255)
                                                                                  tgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                                <  3:396553/149‑1 (MQ=255)
                                                                                   gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                               <  3:150092/149‑1 (MQ=255)
                                                                                      aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                          >  2:958330/1‑103 (MQ=255)
                                                                                      aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                          <  1:958330/103‑1 (MQ=255)
                                                                                        tAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTa                                                          <  4:954142/149‑1 (MQ=255)
                                                                                           gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATTGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                          <  4:853287/114‑1 (MQ=255)
                                                                                           gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                          >  3:853287/1‑114 (MQ=255)
                                                                                           gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt                                                       <  4:148369/149‑1 (MQ=255)
                                                                                            aTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTa                                                      >  2:405233/1‑149 (MQ=255)
                                                                                             ttttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTaa                                                     <  4:935800/149‑1 (MQ=255)
                                                                                              tttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAAt                                                    >  1:49999/1‑149 (MQ=255)
                                                                                                   agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                               <  4:311045/149‑1 (MQ=255)
                                                                                                   agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                               >  4:568997/1‑149 (MQ=255)
                                                                                                   agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                               >  3:726585/1‑149 (MQ=255)
                                                                                                   agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                               <  2:472786/149‑1 (MQ=255)
                                                                                                             cGCAGTGCAGACCGAGTTTGACAGCG‑‑‑‑‑ACTGCGGCAGTCTGCAGCACCTGGTTAGTCTGGATCACCCCGGCAGTAATCAACGTATCGGCACTTTCACGCAGAGCATCTGCCCCGTAAAATTTGACCTTACGTAAGTTATTgccgcatatt                                     <  4:864883/149‑6 (MQ=255)
                                                                                                              gCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTg                                    >  1:155157/1‑149 (MQ=255)
                                                                                                                  tGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                <  2:96619/149‑1 (MQ=255)
                                                                                                                   gCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATg                               <  1:1021024/149‑1 (MQ=255)
                                                                                                                     aGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggg                             >  4:816371/1‑149 (MQ=255)
                                                                                                                      gACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                            <  1:1000211/149‑1 (MQ=255)
                                                                                                                         cGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                            <  2:434057/146‑1 (MQ=255)
                                                                                                                         cGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                            >  1:434057/1‑146 (MQ=255)
                                                                                                                           aGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                       <  4:163003/149‑1 (MQ=255)
                                                                                                                           aGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                       <  3:819302/149‑1 (MQ=255)
                                                                                                                           aGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                       <  1:257585/149‑1 (MQ=255)
                                                                                                                                  cAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATcc                >  1:1019126/1‑149 (MQ=255)
                                                                                                                                       aaccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATcc                >  3:205180/8‑149 (MQ=255)
                                                                                                                                           atgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt            >  4:537452/4‑149 (MQ=255)
                                                                                                                                           atgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt            >  3:341755/4‑149 (MQ=255)
                                                                                                                                            aCTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGttt           <  2:1054889/149‑1 (MQ=255)
                                                                                                                                               gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgat        <  4:284268/149‑1 (MQ=255)
                                                                                                                                                   cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaa    <  3:314395/149‑1 (MQ=255)
                                                                                                                                                    aGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAGATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaaa   <  3:505804/149‑1 (MQ=255)
                                                                                                                                                     gTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAt  <  3:816371/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCG‑‑‑‑‑ACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAT  >  NC_000913/1999086‑1999378

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