Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,009,005 A→C 5.2% S19A (TCG→GCG)  rhtB ← homoserine, homoserine lactone and S‑methyl‑methionine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,009,0050AC 5.2% 526.3 / 17.1 148S19A (TCG→GCG) rhtBhomoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base A (70/70);  new base C (2/6);  total (72/76)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 148 total aligned reads displayed.

CGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGGTTT  >  NC_000913/4008860‑4009152
                                                                                                                                                 |                                                                                                                                                   
cGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAg                                                                                                                                                  <  2:623449/149‑1 (MQ=255)
 gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAgg                                                                                                                                                 >  3:1134598/1‑149 (MQ=255)
     aaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAg                                                                                                                                                  >  4:890691/1‑144 (MQ=255)
     aaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAg                                                                                                                                                  <  3:890691/144‑1 (MQ=255)
      aaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCGCCAGGCtctg                                                                                                                                            <  3:1020745/149‑4 (MQ=255)
      aaTAGCGGCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGTCGCCAGGCtctg                                                                                                                                            <  4:108172/149‑4 (MQ=255)
        tAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTaaaa                                                                                                                                          >  1:796406/1‑149 (MQ=255)
           cttCCCCATCGCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGTGGCCAACGCTCCCCGATAACCGTGGTTGAGCGAGGTGCTCATAGGGTTGATTGCACCAGAGCCTGGCTCCAGGCTTAAAATGa                                                                                                                                       <  4:966514/147‑1 (MQ=255)
            gTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                      <  2:519913/149‑1 (MQ=255)
            gTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                      >  1:519913/1‑149 (MQ=255)
            gTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGAt                                                                                                                                      <  4:811120/149‑1 (MQ=255)
             tCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATc                                                                                                                                     <  4:365987/149‑1 (MQ=255)
               cccAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGa                                                                                                                                   >  2:818950/1‑149 (MQ=255)
                ccAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGAt                                                                                                                                  >  4:699434/1‑149 (MQ=255)
                ccAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGAt                                                                                                                                  <  4:1095295/149‑1 (MQ=255)
                ccAACCCCACGCCAACCAGCACAATATCAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGAt                                                                                                                                  <  4:1134856/149‑1 (MQ=255)
                     cccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcag                                                                                                                             >  2:319648/1‑149 (MQ=255)
                      ccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATc                                                                                                                                     <  2:497892/140‑1 (MQ=255)
                      ccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATc                                                                                                                                     <  2:439647/140‑1 (MQ=255)
                      ccACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTcagc                                                                                                                            >  1:885609/1‑149 (MQ=255)
                                cAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCaa                                                                                                                  >  2:807984/1‑149 (MQ=255)
                                  gCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAac                                                                                                                >  3:659335/1‑149 (MQ=255)
                                     caATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccac                                                                                                             <  2:702090/149‑1 (MQ=255)
                                      aaTATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAaccacc                                                                                                            >  2:737087/1‑149 (MQ=255)
                                       atatGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAAccacca                                                                                                           <  1:446984/149‑1 (MQ=255)
                                         atGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCAtt                                                                                                         <  3:279202/149‑1 (MQ=255)
                                          tGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTc                                                                                                        <  2:1145488/149‑1 (MQ=255)
                                              tCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAg                                                                                                    >  3:277483/1‑149 (MQ=255)
                                              tCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAg                                                                                                    <  1:656242/149‑1 (MQ=255)
                                                gCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGt                                                                                                  <  4:507699/149‑1 (MQ=255)
                                                 ccAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTc                                                                                                 >  3:524181/1‑149 (MQ=255)
                                                   aGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCat                                                                                               <  4:659335/149‑1 (MQ=255)
                                                   aGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCat                                                                                               <  2:1006049/149‑1 (MQ=255)
                                                    gTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCatg                                                                                              >  2:40610/1‑149 (MQ=255)
                                                     tCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAg                                                                                                                      <  1:527405/124‑1 (MQ=255)
                                                     tCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCatga                                                                                             <  2:248211/149‑1 (MQ=255)
                                                        ggTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCAtgatga                                                                                          >  1:420680/1‑149 (MQ=255)
                                                            tGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTc                                                                                      >  1:612691/1‑149 (MQ=255)
                                                             gAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTcc                                                                                     >  4:398943/1‑149 (MQ=255)
                                                              aaGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTccc                                                                                    >  1:404156/1‑149 (MQ=255)
                                                              aaGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTccc                                                                                    >  2:894540/1‑149 (MQ=255)
                                                                 cccAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgt                                                                                 <  1:1077714/149‑1 (MQ=255)
                                                                  ccAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAac                                                                                                                <  2:619747/117‑1 (MQ=255)
                                                                  ccAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAac                                                                                                                >  1:619747/1‑117 (MQ=255)
                                                                  ccAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtg                                                                                >  3:234072/1‑149 (MQ=255)
                                                                   cAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                               <  3:699434/149‑1 (MQ=255)
                                                                   cAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                               <  3:870711/149‑1 (MQ=255)
                                                                   cAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                               <  1:737087/149‑1 (MQ=255)
                                                                         tAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAg                                                                                                    <  3:381548/122‑1 (MQ=255)
                                                                              gCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACa                                                                                                                                                   >  3:165252/1‑70 (MQ=255)
                                                                              gCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACa                                                                                                                                                   <  4:165252/70‑1 (MQ=255)
                                                                               ccACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtgt                                                                   <  3:1000544/149‑1 (MQ=255)
                                                                               ccACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtgt                                                                   <  1:1021397/149‑1 (MQ=255)
                                                                               ccACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtgt                                                                   <  1:457975/149‑1 (MQ=255)
                                                                               ccACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtgt                                                                   <  2:367992/149‑1 (MQ=255)
                                                                                 aCCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGc                                                                 >  2:774521/1‑149 (MQ=255)
                                                                                 aCCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGc                                                                 <  2:796406/149‑1 (MQ=255)
                                                                                 aCCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGc                                                                 >  3:361879/1‑149 (MQ=255)
                                                                                  ccgcgccgcgATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGcc                                                                >  3:606149/1‑149 (MQ=255)
                                                                                        cgcgATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTAttttt                                                                     >  3:459419/1‑138 (MQ=255)
                                                                                         gcgATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCatga                                                                                             <  4:584892/113‑1 (MQ=255)
                                                                                           gATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGc                                                       <  1:314395/149‑1 (MQ=255)
                                                                                           gATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGGCTATTTTTGTGCCACAATACGc                                                       >  3:327650/1‑149 (MQ=255)
                                                                                             tAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTa                                                     <  4:44962/149‑1 (MQ=255)
                                                                                               aCCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACt                                                   >  2:537380/1‑149 (MQ=255)
                                                                                                     gTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgt                                                                                 <  2:301165/113‑1 (MQ=255)
                                                                                                     gTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgt                                                                                 >  1:301165/1‑113 (MQ=255)
                                                                                                     gTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTa                                                     >  2:1131916/1‑141 (MQ=255)
                                                                                                     gTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTa                                                     <  1:1131916/141‑1 (MQ=255)
                                                                                                     gTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCa                                             <  4:277483/149‑1 (MQ=255)
                                                                                                       tGAGCGAGGTGGTCATACTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGCGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGc                                           <  4:967716/149‑1 (MQ=255)
                                                                                                        gAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCg                                          >  4:216641/1‑149 (MQ=255)
                                                                                                           cGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtg                                       <  4:234072/149‑1 (MQ=255)
                                                                                                            gAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgt                                      <  1:220761/149‑1 (MQ=255)
                                                                                                            gAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgt                                      <  3:398943/149‑1 (MQ=255)
                                                                                                               gtggtCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATa                                                          <  2:1075385/126‑1 (MQ=255)
                                                                                                                  gtCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAgg                                <  1:774521/149‑1 (MQ=255)
                                                                                                                     aTAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCAc                             >  3:987694/1‑149 (MQ=255)
                                                                                                                      tAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACg                            >  2:275994/1‑149 (MQ=255)
                                                                                                                         tgtTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTtgt                                      <  3:916110/136‑1 (MQ=255)
                                                                                                                         tgtTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTa                         >  1:1017033/1‑149 (MQ=255)
                                                                                                                         tgtTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTa                         >  1:272753/1‑149 (MQ=255)
                                                                                                                             gATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAAc                                                                                        >  3:252813/1‑82 (MQ=255)
                                                                                                                             gATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTgtgt                                    >  1:800084/1‑134 (MQ=255)
                                                                                                                              aTTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaaaaa                    <  1:818950/149‑1 (MQ=255)
                                                                                                                                tGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaaaaaaa                  <  2:150970/149‑1 (MQ=255)
                                                                                                                                  cACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACg                >  1:1063100/1‑149 (MQ=255)
                                                                                                                                   aCCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaaaaaaa                  >  2:938070/1‑146 (MQ=255)
                                                                                                                                   aCCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTaaaaaaaa                  <  1:938070/146‑1 (MQ=255)
                                                                                                                                     ttaagCCTGTCGCCAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtg                                                                    >  3:731327/4‑94 (MQ=255)
                                                                                                                                     ttaagCCTGTCGCCAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTtg                                                                    <  4:731327/91‑1 (MQ=255)
                                                                                                                                     cAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtt             >  1:496247/1‑149 (MQ=255)
                                                                                                                                     cAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGAtt             >  1:801961/1‑149 (MQ=255)
                                                                                                                                       gagCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGtt                                        <  1:692353/120‑1 (MQ=255)
                                                                                                                                       gagCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGtt                                        >  2:692353/1‑120 (MQ=255)
                                                                                                                                        agCCTGTCGCCAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                               <  1:585398/80‑1 (MQ=255)
                                                                                                                                        agCCTGTCGCCAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGgtgt                                                                               >  2:585398/1‑80 (MQ=255)
                                                                                                                                            tGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTg      <  4:606149/149‑1 (MQ=255)
                                                                                                                                                gACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGGttt  >  1:496322/1‑149 (MQ=255)
                                                                                                                                                gACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGGttt  <  3:997030/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGGTTGAGCGAGGTGGTCATAGTGTTGATTGCACCAGAGCCTGGCGACAGGCTTAAAATGATCGATGTCAGCAGGTAGGCAAACCACCATTCTAAGGTCATGATGAACTCCCGGTGTGTCTATTTTTGTGCCACAATACGCTACTGTCGCAGCGTTGTGTCAGGCACGCTAAAAAAAACGATTTTACGTGGTTT  >  NC_000913/4008860‑4009152

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