Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 6.6% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide/tripeptide:H(+) symporter DtpC/lysine decarboxylase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT 6.6% 275.5 / 5.1 99intergenic (‑34/+25)dtpC/cadAdipeptide/tripeptide:H(+) symporter DtpC/lysine decarboxylase 1
Reads supporting (aligned to +/- strand):  ref base A (49/43);  new base T (3/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCAT  >  NC_000913/4356314‑4356586
                                                                                                                                    |                                                                                                                                             
agcaTGAGTAAGGCACGCGTGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                                                                    <  1:16703/140‑1 (MQ=255)
agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTGTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                                                                    >  2:219965/1‑143 (MQ=255)
  aaTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTTTGTTCCCAAAAGGGAAGTGGCAAGCCACTTCCCtt                                                                                                                                  >  1:289253/1‑143 (MQ=255)
   aTGGGTAAGGCACGCAGGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                                                                 <  2:463162/143‑1 (MQ=255)
      gggaaGGCACGCATGCCGGAAAAACTGGAGTACTCCCAGAGTTGGATCGCCACGATATAGTATATCGCGCGCGGCGTTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttt                                                                                                                              <  1:520249/140‑5 (MQ=255)
      aGTAAGGCACGAATGCCGTAAAAACGGAAGTGCTCCCAGATTTGGATCGCCACGATATAGGATATCGCGCGCGGCTGTGAGAGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttt                                                                                                                              <  1:225265/143‑5 (MQ=255)
       ggaaGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGGGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttg                                                                                                                             <  2:540753/141‑6 (MQ=255)
          aGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                                                          <  1:227146/143‑1 (MQ=255)
          aGGCACGCATGCCGCAAAAGCCCAAGTACCCCCAGATCGGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAg                                                                                                                           <  1:794326/143‑1 (MQ=255)
            gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCGGTGAGGGTGTTTTCATGTGTTCTCCTTATTTTCAAATATGTCACGTGGCTAGCCACTTCCCTTCTTCGAGcgt                                                                                                                        >  1:270673/1‑141 (MQ=255)
               cGCCTGCCGTAAAAACTGAACTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatt                                                                                                                     <  1:380252/143‑1 (MQ=255)
                gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                                                                    <  1:316646/143‑1 (MQ=255)
                       gTATAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg                                                                                                             <  2:120303/143‑1 (MQ=255)
                       gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg                                                                                                             <  1:691741/143‑1 (MQ=255)
                       gTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg                                                                                                             <  2:392351/143‑1 (MQ=255)
                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGGTATAGTATATCGCGCGCGGCTGTGAGGGTGAGTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCACCTTCCCTTGTACGAGCTAATGAGTTTTGGc                                                                                                            >  2:424146/1‑143 (MQ=255)
                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGATAAAATATGTTCTCCTTATGAGCAAAAAATTGAAGTGGTTATCTACTTCCCTTGTATGAGCTAATTATTTTTTGc                                                                                                            >  1:214471/1‑143 (MQ=255)
                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGAATACATGTATTCTCCTTATGAGCAAAAAAGGGAAGTGGCAATGCACTTCCCTTGTACGCGCTAAGTATTTTTTGc                                                                                                            >  1:764192/1‑143 (MQ=255)
                        tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGGGTTTTCATGTGTTCTCCTTATGAGCAAAGAAGGGAAGTGGCAAGCCACTTCCCGTGTGCGAGCTAAATATTTTTTGc                                                                                                            >  2:303344/1‑143 (MQ=255)
                          aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGGACGAGCTAATTATTTTTTGCtt                                                                                                          >  2:685951/1‑143 (MQ=255)
                                gAAGTACCCCCAGATTTGGATCGCCACGCTATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct                                                                                                    <  1:61879/143‑1 (MQ=255)
                                      cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTCGTCATGGGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCCTGTACCACATAACTATGAGTAGCTCTCTTCGTTCaac                                                                                              >  2:605478/1‑142 (MQ=255)
                                       taccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATa                                                                                             >  2:155819/3‑143 (MQ=255)
                                        cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGTGGGTGTTTTCATGAGTTCTCCTTATGAGCAACCAAGGGAAGTGGCAAGCCACTTCCCTTTTACGAGCTAATTATTTTTTGCTTTCTTCTTTCactcc                                                                                            >  1:204406/1‑139 (MQ=255)
                                         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGCGGGTGTTTACATGTGTACACCGTATGAGCTATGAGGGGAAGTGGCAAGCCACTTCTCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                                                           >  1:517031/1‑143 (MQ=255)
                                         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCCCTTCCCTTGTACGAGCTAATTCCTTTTTGCTTTCCTCTTCCCATAcc                                                                                           >  2:294578/1‑143 (MQ=255)
                                         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGAGTTCTGCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCACTAcc                                                                                           >  1:722855/1‑143 (MQ=255)
                                         ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCGGTGAGGGTGTTTTCATGTGTTCTCCCGATGAGCAAAAAAGGGAAGTGGCAAGCCACGGCCCTGGTACGCGCGGATGCTTTTTTCCTTTCTTAtttagacagc                                                                                           >  2:38534/1‑136 (MQ=255)
                                              tttGGATCGCCACTATACAGTATATCGCGCTCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                      <  2:699563/143‑1 (MQ=255)
                                              tttCGATCGCCACGATATAGTATATCGCGCGTGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                      <  2:71480/143‑1 (MQ=255)
                                              ctcccaTCGCCACGAGATACTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAACAACGGGACGCGGCACGACCGTTCCGTTGTACGAGCTGATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                      <  2:765541/138‑1 (MQ=255)
                                                tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTTCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTGTTCTTTCAACACCTTAACgg                                                                                    >  2:576589/1‑143 (MQ=255)
                                                    tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTGTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTATTTGCTGTCTTCTTTCACTACCTTAACGGtata                                                                                >  2:745015/1‑143 (MQ=255)
                                                     cgtcACGATATAGTGTAACGCGGTAGGCTGTGAGTGTGTTTTCATGTGTTCTCCTTATGAGCAAATAAGGGAGGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAg                                                                               <  2:612178/140‑1 (MQ=255)
                                                     cGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAg                                                                               <  1:470546/143‑1 (MQ=255)
                                                       ctaCGATATAGTATATCGTGCGCGGCTGTGAGGGTGTTATCATGTGTTCTCCATATGAGAAAAAAACAGAAGTGGCAAGCCACTTCCCTTGTACGAGCAAATCATTATATGCTTTCGTCTTTCAATAACTTATCGGTACAgag                                                                             >  2:180254/3‑141 (MQ=255)
                                                       ccACGATATAGTATATCGCGCGCGGCTGTGCGGGTGCTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAGGTGGCAAACCCCTTCCCTTGTACGTGCTAATTATTTTGTGCTGTCTGCTTGCACTGCCATAACAGTATAGCg                                                                             >  2:262518/1‑143 (MQ=255)
                                                            atatAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGACGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCAt                                                                        >  1:255217/1‑143 (MQ=255)
                                                             tataGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGccctg                                                                       >  1:65221/1‑140 (MQ=255)
                                                              ctcggatatCGTGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAGAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCATTACCTTAACGGTATAGCGGCCATCa                                                                      <  2:501399/138‑1 (MQ=255)
                                                              ataGTATATCGAGCGGGGATTTGAGGAGCTTTTCATCTTTTCTCCTTATGAGCAAAAGAGGGAAGTGGCAATCCACTTCCCTTGTACAAGCTAATTATTGTATATTGTCTTGTTTCTATACCTTAACGGTATAGCGGCCATCa                                                                      <  2:60215/143‑1 (MQ=255)
                                                                aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                                                                    <  1:641435/143‑1 (MQ=255)
                                                                  tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                                                                  >  1:825797/1‑143 (MQ=255)
                                                                  tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                                                                  >  1:204805/1‑143 (MQ=255)
                                                                   atatCGCGCGCGGCTGAGAGGGTGTTTACATGTGTTCTCCTTATGAGCAAAAAAGGGAAGGGGCACGCCCCTTCCCTTGTACGAGCTAATCATTTTTTGCTTACATCTTTCAATACCTTCACGGTACCACGGccaccacgctg                                                                 >  2:815852/1‑138 (MQ=255)
                                                                       cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                                                             >  1:717781/1‑143 (MQ=255)
                                                                       cgcgcgcgGCTGTGAGGGTGTGTTCCTGTGTTCTCCTTACGCGCAAACCCGGCAAGTGGCTAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACgg                                                             <  2:491972/143‑1 (MQ=255)
                                                                          gcgcgGCCGGGAGGGTGTTTTCCGGTGGTGTTCTTATGAACAAAGAAGGGGAGTGGCAAGCCACTTCCCTTGTACGGTCTAATTATTTTTTGCTTTCCTCTTTCAATACCTTTGCGGTATAGGGGCCATCAGCCTGACGGTAt                                                          <  2:331356/143‑1 (MQ=255)
                                                                            gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc                                                        >  1:282740/1‑143 (MQ=255)
                                                                            gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTGTTGCTTTCTTCTTTCAATACCTTAACGGTAGAGCGGCCATCGGGGTGACGGTGTGc                                                        >  1:104252/1‑143 (MQ=255)
                                                                             cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACCGTCTAGCGGCCCTCAGCCTGCCGGTACGCa                                                       >  1:241989/1‑143 (MQ=255)
                                                                              ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCCTCAGCCTGACGGTCTGCAc                                                      >  1:291686/1‑143 (MQ=255)
                                                                              ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc                                                      <  1:413871/143‑1 (MQ=255)
                                                                                cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCg                                                    >  1:736117/1‑143 (MQ=255)
                                                                                   tgAGGGTGTTGTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa                                                 <  2:315121/143‑1 (MQ=255)
                                                                                      gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCACTACCTTAACGATGTAACGCCCCTAAGCCTGACCGTATCCACCTTGCata                                              >  1:403227/1‑143 (MQ=255)
                                                                                       ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAatgt                                             >  1:736619/1‑141 (MQ=255)
                                                                                         tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGTTGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg                                           <  2:360622/143‑1 (MQ=255)
                                                                                         tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg                                           >  1:53333/1‑143 (MQ=255)
                                                                                           ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTGTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                         >  2:10743/1‑143 (MQ=255)
                                                                                             ttcgggtgtTCGTCTTATGGGCGAAAAGGTGAAGTGGTAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGttt                                       <  2:642545/138‑1 (MQ=255)
                                                                                             ttaaTGTGTTCTTCGTATGATCAAAAAAGGGGATTGGCATGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGttt                                       <  2:123855/140‑1 (MQ=255)
                                                                                             ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATCCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGttt                                       <  1:330319/143‑1 (MQ=255)
                                                                                              tCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGAAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                      >  1:532861/1‑143 (MQ=255)
                                                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCaaa                                   <  1:719637/143‑1 (MQ=255)
                                                                                                 tgtgTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCaaa                                   <  1:718870/143‑1 (MQ=255)
                                                                                                  gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAg                                  <  1:80580/143‑1 (MQ=255)
                                                                                                   tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGc                                 >  1:205363/1‑143 (MQ=255)
                                                                                                   tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGCTTCAAAGc                                 >  1:205733/1‑143 (MQ=255)
                                                                                                    gtTCTGGTTAGGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTGAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGcc                                <  2:573617/143‑1 (MQ=255)
                                                                                                        gccctATGAGAAAAAAAGTGAAGTGGCAAGCCACCTCCCTTGTACCAGCTAATTATTTTTTGCTTTCTTCTTTCAAGATCTAAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTTAATATCGGTTTCAAAGCCCGGa                            <  2:748893/139‑1 (MQ=255)
                                                                                                          cTTGTGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                          <  2:204406/143‑1 (MQ=255)
                                                                                                            taggAGCAAAAAAGGGACGTGGCAAGCCACTTCTCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGt                        <  2:140911/140‑1 (MQ=255)
                                                                                                            tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTGTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTGTGCACCGTGAATATCGGTTTCAAAGCCCGGAtggg                        >  1:288742/1‑140 (MQ=255)
                                                                                                             aagAGCAACACCGGGAAGTGGCAACCCAACTCCCTTGCACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGCAGGTACGTATCGTCAATATCGGTTTCAAAGCCCGGATAGTg                       <  1:384814/141‑1 (MQ=255)
                                                                                                             aTGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                     <  1:601787/81‑1 (MQ=255)
                                                                                                             aTGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                     >  2:601787/1‑81 (MQ=255)
                                                                                                               gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAg                     >  1:217929/1‑143 (MQ=255)
                                                                                                                aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTgcgc                    >  1:798558/1‑143 (MQ=255)
                                                                                                                aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCAACCAACAGCCTGACAGTATACACCGTGAATATCGGTTTCAAAGCCCGGATAGTGacc                    >  1:702600/1‑141 (MQ=255)
                                                                                                                  cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATAACGGTTTCAACGCCCGGACCGTgcgcgc                  >  2:277926/1‑143 (MQ=255)
                                                                                                                    caaagaGGGAAGTGGCAGTGCACGTCCCTTGTCCGAGCTAATTATTTGTTGCTCTCTTCTTTCACTGCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGcc                 <  2:736581/142‑1 (MQ=255)
                                                                                                                    agaaacGGGACGTGGCATGCGACTTCCCTTGCACGCGCTACTTATTCTTTGCTTGCTTCCTCCAATACCTTAACGGTATAGCGGCCTTCAGCCTGACGGTATGCAGCGTGAATATCTGTTTCAAAGCCCGGATAGTGAGCGcc                 <  2:401995/141‑1 (MQ=255)
                                                                                                                      aaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCCAAGCCCGGATAGTGAGCGCCGa               >  1:373363/1‑143 (MQ=255)
                                                                                                                      aaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTGTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGCATGCACCGTGAATATCGGTTTCAAGGCCCGGGTAGTGAGCGCCGa               >  2:417872/1‑143 (MQ=255)
                                                                                                                       aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGTGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGCATAGTGAGCGCCGAt              >  2:74582/1‑143 (MQ=255)
                                                                                                                       aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGACTAGCGGTGTCAAAGCCCGGATGGTGAGCGCCGag              >  1:290145/1‑142 (MQ=255)
                                                                                                                         aGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGAttt            >  1:93840/1‑143 (MQ=255)
                                                                                                                          gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGGTTCAAAGCCAGGATATGCAGCGCCGATTTc           >  2:147631/1‑143 (MQ=255)
                                                                                                                           ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca          >  1:437165/1‑143 (MQ=255)
                                                                                                                           ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCGTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCGCGGATAGTGAGCGCCGATTTca          <  2:722855/143‑1 (MQ=255)
                                                                                                                           ggAAGGGGCAAGGCACTTCCCTTGTACGAGCTAATTATTTGTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCGGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca          <  1:786090/143‑1 (MQ=255)
                                                                                                                             gaGTGGCAAGCCACTTCCCTTGTACGATCTAATTATTTTTTGCTTTCTTCTTTCATTACCTTAATGGTATAGCGGCCATCAGTCTGACAGTATGCACCGTGAATATCGGTTTCATAGCCCGGATAGTGAGCGCCGATTTcaca        <  2:665800/142‑1 (MQ=255)
                                                                                                                             caGTCGCAAGCCACTTCTCTTGTACGATCTAATTATTTTTTGCTTTCTCCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTTACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcaca        <  2:764192/142‑1 (MQ=255)
                                                                                                                              aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacac       <  2:769815/143‑1 (MQ=255)
                                                                                                                                tGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAg     >  1:690698/1‑143 (MQ=255)
                                                                                                                                 gtcAAGCCACCTCCCTCGTACGAGCTAGTTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGc    <  2:816536/141‑1 (MQ=255)
                                                                                                                                 ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATCGTGAGCGCCGATTTCACACAGc    >  1:23269/1‑143 (MQ=255)
                                                                                                                                   cAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCAt  <  1:88529/143‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                             
AGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCAT  >  NC_000913/4356314‑4356586

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