Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 176,979 C→T 5.1% intergenic (+25/+22) erpA → / ← yadS iron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913176,9790CT 5.1% 1005.4 / 27.9 313intergenic (+25/+22)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (154/138);  new base T (5/11);  total (162/151)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 334 total aligned reads displayed.

TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128
                                                                                                                                                     |                                                                                                                                                     
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                       >  2:4103193/1‑150 (MQ=255)
    ccTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATc                                                                                                                                                   >  1:4763844/1‑150 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                 >  2:1662514/1‑150 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                 >  2:2250884/1‑150 (MQ=255)
       gCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCgg                                                                                                                                                <  2:1973456/150‑1 (MQ=255)
       gCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCgg                                                                                                                                                <  1:5799681/150‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                              >  2:4167085/1‑150 (MQ=255)
            atCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGCACAt                                                                                                                                           <  1:5541496/150‑1 (MQ=255)
               tGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGGGGTTCTTCCTTTAGTATCTAATCTGTTGTTCTTGTGCCGGATAAACTATCCGGCACAtctt                                                                                                                                        <  1:488371/150‑1 (MQ=255)
               tGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtctt                                                                                                                                        <  1:392906/150‑1 (MQ=255)
                ggTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATAGTGACCAACCCGAAAGCGAAAAGCACCTGCGGTTGCGGTTCTTGCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtcttc                                                                                                                                       >  1:790198/1‑150 (MQ=255)
                   cggcggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctct                                                                                                                                    <  2:2486786/150‑1 (MQ=255)
                      cggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCtt                                                                                                                                 >  2:1507731/1‑150 (MQ=255)
                                 aTTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                      <  2:3285606/150‑1 (MQ=255)
                                 aTTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                      >  1:4309100/1‑150 (MQ=255)
                                   tataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                    >  2:3950062/1‑150 (MQ=255)
                                    agcccGAAGGTCTGGATGGTTCTCGTTTCATCGTGACCACCCCGACCGCGAAAAGCCCCTGCGGTTGCGGTTCTTCGTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtctc                                                                                                                   <  2:3652118/147‑1 (MQ=255)
                                     taCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCa                                                                                                                  >  2:4468167/1‑150 (MQ=255)
                                     taCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCATCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCa                                                                                                                  <  2:4967263/150‑1 (MQ=255)
                                      aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCAt                                                                                                                 >  1:3355735/1‑150 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                >  1:620657/1‑150 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                >  1:126801/1‑150 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                <  2:1312792/150‑1 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                               >  1:461910/1‑150 (MQ=255)
                                         gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                              <  2:2889336/150‑1 (MQ=255)
                                           aggaatGGAAGTTTCACGTTTCATCGTGACCAACCCGAACGCGAAAACCACCTGCGTTTGCGGTACTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                            <  2:5498510/145‑1 (MQ=255)
                                             gTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgcgc                                                                                                          >  2:4441199/1‑150 (MQ=255)
                                                 ggAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                      >  2:3480135/1‑150 (MQ=255)
                                                  gAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACg                                                                                                     >  2:1014148/1‑150 (MQ=255)
                                                   aaGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGGGGTTCTTCCTTTTGTATCTAATCTGTTGTTCTTGTGCCGTTTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGt                                                                                                    <  1:5364091/150‑1 (MQ=255)
                                                       ttCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGt                                                                                                >  1:2470767/1‑150 (MQ=255)
                                                       ttCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGt                                                                                                >  1:1125497/1‑150 (MQ=255)
                                                        tctgttttAATGGTGACCAACCCGAACGCTAAAAGCACGTGGGGTTGGGTTTCTTCATTTTTTATATAATGTGTTTTTATTGTGCGTTATAAATTACCCGGAACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTa                                                                                               <  1:5899932/145‑1 (MQ=255)
                                                         ctcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAg                                                                                              >  1:2854559/1‑150 (MQ=255)
                                                            gTTTAATCTTGACCATCCCGAAGTCTAAAAGCACCTGGGGTTGCGTTTCTTACTTTAGTATATATTCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCTGTAGCtt                                                                                           <  2:2950297/150‑1 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCtga                                                                                          >  2:3749148/1‑148 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                          >  2:4825901/1‑150 (MQ=255)
                                                              tacATCGTTACCAACCCCAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCGGGATAAACTATCCGCCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTaa                                                                                         <  1:3662272/148‑1 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATg                                                                                     <  2:5978696/150‑1 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATg                                                                                     >  1:2030163/1‑150 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATg                                                                                     >  2:2529799/1‑150 (MQ=255)
                                                                   ggTGACCAACCCTAACGCTAAAAGCACCTGCGGTTGGGTTTCTTCATTTTTTATCTAATCTGTTTTTCATGTGCCGTATAAACTAACCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                    <  1:1423612/149‑1 (MQ=255)
                                                                      gACCACCCCGACCGCGAACCGCACGTGGGGTGGGGGTTCTTCCTTTTGTCTCTAATCGTTTGTTCTTGGGCGGGATAAACTATCGGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCaa                                                                                 <  1:4693744/150‑1 (MQ=255)
                                                                       accaccCCGACCGCGAAAAGCACCTGCGGTTGGGGTTCTTCCTTTATTATCTAATCTGTTGTTCATGTGCCGGATAAACTACCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAAc                                                                                <  1:4363180/150‑1 (MQ=255)
                                                                         cAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                              >  1:933427/1‑150 (MQ=255)
                                                                          accccGACCGCTAAACGCACCTGGGGTTGCGGTTCTTCCTTTTTAATCTAATCTGTTTTTCTTGGGCCGGATAAACTATCCGGCAAATCTTCTCTCTTCAACGCCCATTCTCATCCCGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                             <  1:2088974/148‑1 (MQ=255)
                                                                            cccGAACGCGAAAAGCACCTGCGGTTGCGGTTATTCCTTTAGTATCTAATCTTTTTTTCATGTGCCGTATAAACTATCCGGCAAATCTTCTCTCTTCAACGCCCATTCCCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATa                                                                           <  1:4528206/150‑1 (MQ=255)
                                                                              cGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgc                                                                         >  1:1544774/1‑150 (MQ=255)
                                                                               gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                        >  1:524839/1‑150 (MQ=255)
                                                                               gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                        <  2:3438973/150‑1 (MQ=255)
                                                                               gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                        >  1:1815834/1‑150 (MQ=255)
                                                                                 aCGCTAAATGCACCTTCTTTTGCGGTTCTTCCTTTATTATATAATCTTTTGTTCATGTGCCGTATAAACTATCCGGAAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGAGCAAACGTCGTTAGCTTAAGATGCCAACGAATAgccgc                                                                      <  2:1225745/150‑1 (MQ=255)
                                                                                     gAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATATGTTGTTCATTTGCCGGATAAACTAACCGGCACATCTTCTCGCTTCAACGCCCAATCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGc                                                                  >  1:2616543/1‑150 (MQ=255)
                                                                                      aaaaGCACCTGGGGTTGGGGTTCTTCTTTTATTATCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGcc                                                                 <  1:5326921/150‑1 (MQ=255)
                                                                                        aaGCACCTGCGGTTGCGTTTCTTCCTTTAGTATCTAATGTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGa                                                               <  1:4229193/150‑1 (MQ=255)
                                                                                        aaGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGa                                                               >  1:4940144/1‑150 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAAt                                                         <  1:3170515/150‑1 (MQ=255)
                                                                                                    ttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTc                                                   >  1:5376753/1‑150 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCCACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAc                                              >  1:3991864/1‑150 (MQ=255)
                                                                                                             ttcCTTTTGTATGTAATGTGTTGTTCTTGGGCGGGATAACCTATCGGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATg                                          <  1:5442128/150‑1 (MQ=255)
                                                                                                             ttcCTTTATTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATg                                          <  1:122086/150‑1 (MQ=255)
                                                                                                              tcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                         >  2:4231240/1‑150 (MQ=255)
                                                                                                                cTTTAGTTTCTAATGTGTTGTTCTTGTGCCGGATAAACTATCGGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGccc                                       <  2:2596098/150‑1 (MQ=255)
                                                                                                                 tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                      <  2:708221/150‑1 (MQ=255)
                                                                                                                 tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                      >  1:1027221/1‑150 (MQ=255)
                                                                                                                 tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                      >  1:1836209/1‑150 (MQ=255)
                                                                                                                   tAGTATCTAATCTGTTGTTCTTGTGCCGTATAACCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatc                                    <  2:2130274/150‑1 (MQ=255)
                                                                                                                      tATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATa                                 >  2:3189733/1‑150 (MQ=255)
                                                                                                                         cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                              <  2:595982/150‑1 (MQ=255)
                                                                                                                         cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                              <  1:2523038/150‑1 (MQ=255)
                                                                                                                            aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCt                           <  2:790198/150‑1 (MQ=255)
                                                                                                                             tCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                          >  1:91359/1‑150 (MQ=255)
                                                                                                                              cTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                         <  1:5584775/150‑1 (MQ=255)
                                                                                                                              cTGTTGTTCATGTGCCGTATAAACTAGCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                         <  1:4103193/150‑1 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                        >  1:3843418/1‑150 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                        >  1:4077471/1‑150 (MQ=255)
                                                                                                                                gttgttCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGttt                       <  2:4533913/150‑1 (MQ=255)
                                                                                                                                gttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGttt                       >  2:2453576/1‑150 (MQ=255)
                                                                                                                                 ttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                      >  2:5599569/1‑150 (MQ=255)
                                                                                                                                   gttCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  2:5261224/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  2:5419564/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  1:5938662/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    >  1:3437254/1‑150 (MQ=255)
                                                                                                                                    ttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                   <  1:5835152/150‑1 (MQ=255)
                                                                                                                                      cATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTg                 <  1:4179103/150‑1 (MQ=255)
                                                                                                                                      cATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTg                 <  2:2748668/150‑1 (MQ=255)
                                                                                                                                         gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa              <  1:2358728/150‑1 (MQ=255)
                                                                                                                                          tgCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             >  1:1937422/1‑150 (MQ=255)
                                                                                                                                          tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCAAGCCCATCATACTGGCTGTTTCCAGTGGTAc             >  1:5669643/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            <  2:933427/150‑1 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            >  1:5397363/1‑150 (MQ=255)
                                                                                                                                            ccGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg           >  2:4063354/1‑150 (MQ=255)
                                                                                                                                              ggATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaa         >  1:5084859/1‑150 (MQ=255)
                                                                                                                                               gATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCCTTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaa        >  1:3598155/1‑150 (MQ=255)
                                                                                                                                                aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGccc                                       <  1:258722/118‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:3821121/150‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:3401626/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATttgt  <  2:1563567/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATgtgt  <  1:1724851/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATgtgt  <  2:3645270/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: