Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,564 G→A 8.1% intergenic (‑222/+69) yddA ← / ← ydeM putative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,5640GA 8.1% 226.6 / 8.6 74intergenic (‑222/+69)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base A (3/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

ACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT  >  NC_000913/1579420‑1579712
                                                                                                                                                |                                                                                                                                                    
aCAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                                                                  <  2:670205/149‑1 (MQ=255)
  aaTATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                                                                                >  2:373782/1‑149 (MQ=255)
  aaTATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                                                                                <  1:488234/149‑1 (MQ=255)
   atatGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                                                               <  3:219214/149‑1 (MQ=255)
      tGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTc                                                                                                                                            >  4:183818/1‑149 (MQ=255)
        aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          >  2:383715/1‑149 (MQ=255)
        aaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          >  2:453603/1‑149 (MQ=255)
          atatAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                                                                        <  4:401597/149‑1 (MQ=255)
             taAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                                                                     <  4:720792/149‑1 (MQ=255)
                aGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                                                                  >  2:287346/1‑149 (MQ=255)
                aGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                                                                  <  2:97719/149‑1 (MQ=255)
                  gAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcgtc                                                                                                                                <  4:470813/149‑8 (MQ=255)
                  gAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGtt                                                                                                                                >  1:718282/1‑149 (MQ=255)
                      atTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                                                            <  4:169884/149‑1 (MQ=255)
                        ttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                                                                          <  1:80449/149‑1 (MQ=255)
                            ttttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                                                                      <  1:383715/149‑1 (MQ=255)
                            ttttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                                                                      <  1:630908/149‑1 (MQ=255)
                              ttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                                                               >  1:178561/1‑122 (MQ=255)
                              ttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                                                               <  2:178561/122‑1 (MQ=255)
                                           gttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                                                                       >  2:10633/1‑149 (MQ=255)
                                             tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                                                                     <  1:453603/149‑1 (MQ=255)
                                                  tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGtt                                                                                                >  3:498492/1‑149 (MQ=255)
                                                   aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                                                               >  2:387908/1‑149 (MQ=255)
                                                      gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                            >  3:622921/1‑149 (MQ=255)
                                                      gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                            <  4:694370/149‑1 (MQ=255)
                                                        tAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAaga                                                                                          <  1:10633/149‑1 (MQ=255)
                                                           tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                                                                       <  1:591960/149‑1 (MQ=255)
                                                                 atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                                                                                 <  1:119897/149‑1 (MQ=255)
                                                                    ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaa                                                                              <  1:594376/149‑1 (MQ=255)
                                                                    ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaa                                                                              <  4:259521/149‑1 (MQ=255)
                                                                      gCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          <  2:203443/87‑1 (MQ=255)
                                                                      gCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          >  1:203443/1‑87 (MQ=255)
                                                                            gaccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                                                                      <  3:617443/147‑1 (MQ=255)
                                                                            gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGCAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                                                                      <  2:664020/149‑1 (MQ=255)
                                                                            gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                                                                      >  2:550170/1‑149 (MQ=255)
                                                                                     cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                                                             <  3:482196/149‑1 (MQ=255)
                                                                                      aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                            >  2:739752/1‑149 (MQ=255)
                                                                                      aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt                                                            <  4:498492/149‑1 (MQ=255)
                                                                                             gTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                                                               >  4:443470/1‑59 (MQ=255)
                                                                                             gTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                                                               <  3:443470/59‑1 (MQ=255)
                                                                                                  ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATa                                                >  2:543773/1‑149 (MQ=255)
                                                                                                   cAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATaa                                               >  3:339151/1‑149 (MQ=255)
                                                                                                    aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAt                                              >  3:424827/1‑149 (MQ=255)
                                                                                                     ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                                                                                 <  4:677086/49‑1 (MQ=255)
                                                                                                     ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                                                                                 >  3:677086/1‑49 (MQ=255)
                                                                                                       aTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtgtg                                           <  1:287346/149‑1 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                                                                    <  4:411112/54‑1 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                                                                    >  3:411112/1‑54 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                    >  2:145692/1‑118 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                                                                    <  1:145692/118‑1 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAg                                     >  4:157014/1‑149 (MQ=255)
                                                                                                             gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAg                                     >  2:408750/1‑149 (MQ=255)
                                                                                                              cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATAGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGc                                    <  2:203836/149‑1 (MQ=255)
                                                                                                               aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCg                                   >  2:693258/1‑149 (MQ=255)
                                                                                                                 gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGt                                 >  4:176791/1‑149 (MQ=255)
                                                                                                                  ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTa                                <  2:228403/149‑1 (MQ=255)
                                                                                                                      gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt                                                                 <  3:330127/112‑1 (MQ=255)
                                                                                                                      gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt                                                                 >  4:330127/1‑112 (MQ=255)
                                                                                                                       atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTg                           <  4:300043/149‑1 (MQ=255)
                                                                                                                         gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                         <  1:373782/149‑1 (MQ=255)
                                                                                                                         gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                         <  2:671504/149‑1 (MQ=255)
                                                                                                                         gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATAAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                         <  1:505825/149‑1 (MQ=255)
                                                                                                                             atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTCATTATTCTCCATTTGTTTTGCAACGTCCATAATGTGGTAAAGCGGTACTCTGGTCTTa                     <  4:663812/149‑1 (MQ=255)
                                                                                                                             atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTa                     <  1:244270/149‑1 (MQ=255)
                                                                                                                                  gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa                >  1:585375/1‑149 (MQ=255)
                                                                                                                                  gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa                >  1:153440/1‑149 (MQ=255)
                                                                                                                                    ggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                         <  3:155223/136‑1 (MQ=255)
                                                                                                                                    ggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGtt                         >  4:155223/3‑138 (MQ=255)
                                                                                                                                    aaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAAct              >  2:147519/1‑149 (MQ=255)
                                                                                                                                      cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTGCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAActct            <  4:390389/149‑1 (MQ=255)
                                                                                                                                       acCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa           >  4:379054/1‑149 (MQ=255)
                                                                                                                                       acCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa           >  1:646579/1‑149 (MQ=255)
                                                                                                                                       acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa           >  2:698216/1‑149 (MQ=255)
                                                                                                                                         ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTATCTAACTCTAcc         >  2:774210/1‑149 (MQ=255)
                                                                                                                                                gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGt  >  4:71354/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT  >  NC_000913/1579420‑1579712

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: