Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,175 T→A 6.9% intergenic (+21/+43) emrB → / ← luxS multidrug efflux pump membrane subunit EmrB/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,1750TA 6.9% 146.1 / 4.2 59intergenic (+21/+43)emrB/luxSmultidrug efflux pump membrane subunit EmrB/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base A (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

TGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTG  >  NC_000913/2814033‑2814316
                                                                                                                                              |                                                                                                                                             
tgcctatcTCCGCCAATGCGTTCTTCTGGCTGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTCGTGTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                               <  2:129833/139‑1 (MQ=255)
tGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                               >  1:839092/1‑143 (MQ=255)
tGAGTATTTCCGCCAACGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGGCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                               <  2:875929/143‑1 (MQ=255)
 tatttcttCCGACAATGGTTTCTTCTGGATGTCAGCCGGGTTCTTCCTCGTTCTGCTGCGGCTGGTGTGGTTTGCTAAACCGCCTTTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                              <  2:89249/137‑1 (MQ=255)
  attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTg                                                                                                                                             >  1:281147/1‑143 (MQ=255)
      tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCGAAACCGCAACTCGACGCAAGTGGCGGCGGAGGCGGTGCGCACTGCGTCCAACTAAGCCAGCTCATTTGAACt                                                                                                                                         >  1:665521/1‑143 (MQ=255)
        tCCGCCAATGAGATCTTCTGGTTGTGAGCCGGGATATTCCTCGTCCTTCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                       <  2:872382/143‑1 (MQ=255)
          cGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCt                                                                                                                                     <  1:538478/143‑1 (MQ=255)
                tGAGATCTTCTGGATGTCAGCCGGGATGCTCCTCGTCCTGCTGGGGCTGGTGTGGCTTGCTAAACCGCCATTTGTCGCAGATGGCGGCGGCGGCGGAGGGCACTACGTACAGCTAGGCCGGGCCATTCGAACTGGCtttttgc                                                                                                                               >  2:493967/1‑141 (MQ=255)
                  ctgtcttctGGATGTCACCCGTGATTGTCCTCGTCCTGCTGGGGCTGTTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaa                                                                                                                             <  2:289192/140‑1 (MQ=255)
                      cgtctGGACGCCTGCCGGGATATTCCCCGTCCTGCTGTGGCTGGTGTAGTTTGCTAATCCGCCATTTGCCGCAGGTGGCGGCGGATGCAGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                         <  1:155487/141‑1 (MQ=255)
                       ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTCCGCCAGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGATGCGCACTAAGTACACCTAAGCCAGCTCATTCGAACTGACTTTTGCCaattact                                                                                                                        >  2:650332/1‑141 (MQ=255)
                         ctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtg                                                                                                                      >  2:573548/1‑143 (MQ=255)
                         ctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTGTTTTCAATCAATtg                                                                                                                      >  1:42077/1‑143 (MQ=255)
                            gATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACACCTAAGCCAGTTCATTTGGACTGGCTTTTTTCAATTACTTGTGa                                                                                                                   >  2:765135/1‑143 (MQ=255)
                             aggTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                                                                  <  2:761173/141‑1 (MQ=255)
                                  aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGAGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGCGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTGAAGTAAATTGTGAAGATAg                                                                                                             >  2:336266/1‑143 (MQ=255)
                                      gggATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                                                         <  2:839092/143‑1 (MQ=255)
                                        gATATTCCTCGGCGTGCTGGGGCTGGGGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGGGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                                                       <  2:174981/143‑1 (MQ=255)
                                           atgccgcGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                                                    <  2:154756/137‑1 (MQ=255)
                                           atTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTATGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                                                    <  1:475129/143‑1 (MQ=255)
                                              ccTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAg                                                                                                 >  1:616441/1‑143 (MQ=255)
                                               cTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                                                                >  2:39617/1‑143 (MQ=255)
                                                  gTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgt                                                                                             <  2:307760/143‑1 (MQ=255)
                                                   tCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgtg                                                                                            >  2:828508/1‑143 (MQ=255)
                                                     ctgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATCAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                          >  2:144744/1‑143 (MQ=255)
                                                      tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGATGGCGGCGGAGGCGGAGCGCACTCAGTACAACTAAGCCAGATCATTTAAACTGGCTTGTTTCAAACAATTGTGGAGATAGTTAACTGACTAGAAGTGAAg                                                                                         >  2:450307/1‑143 (MQ=255)
                                                            ggCTGGTGTGGTGAGGTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTg                                                                                   <  1:291285/143‑1 (MQ=255)
                                                               ttgtgtTGTTTGCTAAACAGTCATTTGGCGCAGGTGGCGGCGGATGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCaa                                                                                <  2:578306/141‑1 (MQ=255)
                                                                  tgtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACtt                                                                             >  1:578934/1‑143 (MQ=255)
                                                                    tgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCACTTAATTGTGAAGATAGTTTACTGGCTAGGTGGGCAGTTCCTGCAACTtct                                                                           >  1:369130/1‑143 (MQ=255)
                                                                    tgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtct                                                                           >  2:912862/1‑143 (MQ=255)
                                                                          gCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGTGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCg                                                                     <  1:129405/143‑1 (MQ=255)
                                                                           cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCgg                                                                    >  1:424596/1‑143 (MQ=255)
                                                                           cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCACGTTCAATTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTATTTCgg                                                                    >  2:363362/1‑143 (MQ=255)
                                                                            tAAACGGCCGTTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                   <  2:71989/143‑1 (MQ=255)
                                                                             gatcctgcATTTGCTGCAGGTGCCTGCGGACGCTGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCa                                                                  <  2:177880/136‑1 (MQ=255)
                                                                               aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                                                                >  1:681594/1‑143 (MQ=255)
                                                                                  gCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                                                             >  1:303596/1‑143 (MQ=255)
                                                                                   ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCATTGCCa                                                            >  1:745780/1‑143 (MQ=255)
                                                                                      tttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCCGTGCCAGtt                                                         >  1:685365/1‑143 (MQ=255)
                                                                                          gcgcAGGTGGCGGCGGAGGCGGTGCGCACTAATTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                                                     >  2:476884/1‑143 (MQ=255)
                                                                                           cgcgcggggcggcGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCg                                                    <  2:500482/136‑1 (MQ=255)
                                                                                           cgcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCg                                                    <  1:644217/143‑1 (MQ=255)
                                                                                             cAGGTGGCTGCGCATGCGTGTCTCACTAAGTACAACTAATCCAGTTCATTTGGACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGtt                                                  <  2:834767/143‑1 (MQ=255)
                                                                                              aGGTGGCGGCGGGGGCGGTGCGCACTAAGTACACCTAAGACAGTTCAGTTGGACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                 <  2:643196/143‑1 (MQ=255)
                                                                                                     ggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGa                                          <  1:496871/143‑1 (MQ=255)
                                                                                                        ggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTACTTGTGCAGATTGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGc                                       <  1:83303/143‑1 (MQ=255)
                                                                                                          aGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGt                                     <  1:765135/143‑1 (MQ=255)
                                                                                                               gTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                <  1:758164/143‑1 (MQ=255)
                                                                                                                    cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGATGCCGCTGATGCGTACGTCACGtt                           >  2:130971/1‑143 (MQ=255)
                                                                                                                     aCTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTc                          <  1:723912/143‑1 (MQ=255)
                                                                                                                     aCTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTc                          <  1:218445/143‑1 (MQ=255)
                                                                                                                          gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGaa                     <  2:281147/143‑1 (MQ=255)
                                                                                                                            aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTGCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATg                   <  1:263369/143‑1 (MQ=255)
                                                                                                                            aCAACTAAGCCAGCTCTTGTGAACTGGCGCTTTTCACTTAATTCTGAAGATATCTTTCTGACTAGATTTGCAGTTCCTGCAACTTCTCTTTCTGCAGTGCCAGTTATTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATg                   <  2:533544/143‑1 (MQ=255)
                                                                                                                               aCTATGCGAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTCTACTGACTAGATGTGCAGTCTCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGCTGCTGTTTATGCGTACGTCACGTTCCAGAATGCTa                <  2:725564/143‑1 (MQ=255)
                                                                                                                                 tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTCATGCGTACGTCACGTTCCACCATGCTAcg              >  2:742837/1‑143 (MQ=255)
                                                                                                                                 aaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcg              >  2:6225/2‑143 (MQ=255)
                                                                                                                                  gaGCCAGTTCATTTGAACTGGCTTGTTTCACTTAAGTGTGAAGATAGTCTACTGACTAGATGTGCAGTTCCTGCATCTTATCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgc             <  2:456338/142‑1 (MQ=255)
                                                                                                                                      cAGTTCATTTGAATTGGCTTTTTTCAATTAATGGAGAAGAGAGGTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCaa         <  1:386987/143‑1 (MQ=255)
                                                                                                                                      cAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAGTGCTACGCGCaa         >  2:38302/1‑143 (MQ=255)
                                                                                                                                            aTTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCt   >  2:300798/1‑143 (MQ=255)
                                                                                                                                             tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAatcctctg  >  1:702847/1‑137 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTG  >  NC_000913/2814033‑2814316

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