Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 176,977 A→G 5.8% intergenic (+23/+24) 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,9770AG 5.8% 339.4 / 8.0 104intergenic (+23/+24)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (58/40);  new base G (3/3);  total (61/43)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 112 total aligned reads displayed.

TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATG  >  NC_000913/176830‑177125
                                                                                                                                                   |                                                                                                                                                    
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                     <  1:1033321/149‑1 (MQ=255)
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                     <  1:903561/149‑1 (MQ=255)
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                     <  1:842045/149‑1 (MQ=255)
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                     <  1:527008/149‑1 (MQ=255)
 gAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                    <  1:853139/149‑1 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                               >  1:152164/1‑149 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                               <  1:823664/149‑1 (MQ=255)
           tatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                          <  4:546105/149‑1 (MQ=255)
           tatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                          <  4:41821/149‑1 (MQ=255)
                        gTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTc                                                                                                                             <  2:877175/149‑1 (MQ=255)
                             gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGc                                                                                                                        >  3:423510/1‑149 (MQ=255)
                              ttGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                       >  3:11802/1‑149 (MQ=255)
                                gATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCa                                                                                                                     <  1:1009494/149‑1 (MQ=255)
                                 aTTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAt                                                                                                                    >  4:112799/1‑149 (MQ=255)
                                  ttATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                   <  4:564748/149‑1 (MQ=255)
                                   tataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtc                                                                                                                  <  4:35221/149‑1 (MQ=255)
                                      aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCa                                                                                                               <  2:550051/149‑1 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                             >  1:832276/1‑149 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                             >  4:1005058/1‑149 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                             >  1:21174/1‑149 (MQ=255)
                                         gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                            >  1:338212/1‑149 (MQ=255)
                                          aaGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                           <  1:28082/149‑1 (MQ=255)
                                                tGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtc                                                                                                                  >  2:951407/1‑136 (MQ=255)
                                                tGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCaa                                                                                                     >  4:52675/1‑149 (MQ=255)
                                                  gAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                   >  4:705754/1‑149 (MQ=255)
                                                  gAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                   >  2:461928/1‑149 (MQ=255)
                                                  gAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                   >  1:837994/1‑149 (MQ=255)
                                                   aaGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACg                                                                                                  >  1:275596/1‑149 (MQ=255)
                                                      gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                               <  2:152164/149‑1 (MQ=255)
                                                      gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                               >  2:79425/1‑149 (MQ=255)
                                                      gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                               >  3:603705/1‑149 (MQ=255)
                                                      gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                               <  4:423510/149‑1 (MQ=255)
                                                       ttCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCgg                                                                                              >  3:59732/1‑149 (MQ=255)
                                                                cATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                     <  3:705754/149‑1 (MQ=255)
                                                                 aTCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                    >  1:108744/1‑149 (MQ=255)
                                                                 aTCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                    >  1:869952/1‑149 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtc                                                                                                                                       >  2:233815/1‑97 (MQ=255)
                                                                  tCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtc                                                                                                                                       <  1:233815/97‑1 (MQ=255)
                                                                   cGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                    <  2:957900/147‑1 (MQ=255)
                                                                   cGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                    >  1:957900/1‑147 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                 <  3:978132/149‑1 (MQ=255)
                                                                      gACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCa                                                                               <  1:951447/149‑1 (MQ=255)
                                                                      gACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCa                                                                               >  4:894907/1‑149 (MQ=255)
                                                                       aCCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCATTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACGTCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCCTAAGATGCCaa                                                                              <  1:870025/149‑1 (MQ=255)
                                                                          aaCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                           <  3:39065/149‑1 (MQ=255)
                                                                           aCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                          <  3:52675/149‑1 (MQ=255)
                                                                           aCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                          <  3:182325/149‑1 (MQ=255)
                                                                           aCCCGAACACGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                          >  1:625324/1‑149 (MQ=255)
                                                                             ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATa                                                                        <  2:338212/149‑1 (MQ=255)
                                                                              cGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                              >  3:795592/1‑78 (MQ=255)
                                                                              cGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                              <  4:795592/78‑1 (MQ=255)
                                                                                aaCGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                     >  4:575436/1‑149 (MQ=255)
                                                                                aaCGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                     >  1:194849/1‑149 (MQ=255)
                                                                                     gAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAg                                                                >  2:245164/1‑149 (MQ=255)
                                                                                      aaaaGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGc                                                               >  1:193381/1‑149 (MQ=255)
                                                                                          gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaa                                                           >  1:325025/1‑149 (MQ=255)
                                                                                             ccTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCa                                                        >  1:405596/1‑149 (MQ=255)
                                                                                             ccTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCa                                                        >  4:39870/1‑149 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCGTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCaa                                                       >  1:813508/1‑149 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                   <  1:644603/105‑1 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                   >  2:644603/1‑105 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCaa                                                       >  3:863188/1‑149 (MQ=255)
                                                                                                 cggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAg                                                    <  2:837994/149‑1 (MQ=255)
                                                                                                 cggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAg                                                    <  1:340165/149‑1 (MQ=255)
                                                                                                      gcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCa                                               >  2:757866/1‑149 (MQ=255)
                                                                                                        ggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                          <  4:325125/56‑1 (MQ=21)
                                                                                                        ggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                          >  3:325125/1‑56 (MQ=21)
                                                                                                        ggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCAAACGAATAGCCGCCAGCCGAATCAATAGCGTCACg                                             >  2:175135/1‑149 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAg                                                                                           <  1:16254/102‑1 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAg                                                                                           >  2:16254/1‑102 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGa                                            <  2:275596/149‑1 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGa                                            <  3:819730/149‑1 (MQ=255)
                                                                                                          ctcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAc                                           <  3:93148/148‑1 (MQ=255)
                                                                                                           tcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAcc                                          <  1:191394/149‑1 (MQ=255)
                                                                                                             ttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCAt                                        >  2:834412/1‑149 (MQ=255)
                                                                                                                cTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGcc                                     >  4:12673/1‑149 (MQ=255)
                                                                                                                 tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                 <  1:1041015/104‑1 (MQ=255)
                                                                                                                 tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                 >  2:1041015/1‑104 (MQ=255)
                                                                                                                 tttAGTATCGAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGccc                                    >  2:272473/1‑149 (MQ=255)
                                                                                                                     gTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatca                                >  3:802191/1‑149 (MQ=255)
                                                                                                                        tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATAc                             >  3:355863/1‑149 (MQ=255)
                                                                                                                          tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                           >  2:155099/1‑149 (MQ=255)
                                                                                                                           aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTgg                          >  1:532292/1‑149 (MQ=255)
                                                                                                                            aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                         >  4:387262/1‑149 (MQ=255)
                                                                                                                            aTCTGTTGTTCATGTGCCGGATAAACTATCCGACACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCTTCATACTGGc                         >  4:898325/1‑149 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                      <  2:559154/133‑1 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                      >  1:559154/1‑133 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                      >  4:1042210/1‑149 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                      >  3:1034616/1‑149 (MQ=255)
                                                                                                                                  agaagaTGTGCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                   >  1:257771/6‑149 (MQ=255)
                                                                                                                                    ttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                 <  3:883163/149‑1 (MQ=255)
                                                                                                                                     tCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                <  4:355863/149‑1 (MQ=255)
                                                                                                                                     tCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                >  1:433436/1‑149 (MQ=255)
                                                                                                                                       aTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTg              <  1:461928/149‑1 (MQ=255)
                                                                                                                                        tgtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTgg             >  3:559840/1‑149 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc          >  2:482876/1‑149 (MQ=255)
                                                                                                                                                  aaaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAt                                                                                   <  3:425217/69‑1 (MQ=255)
                                                                                                                                                  aaaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAt                                                                                   >  4:425217/1‑69 (MQ=255)
                                                                                                                                                  aaaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt   <  2:193381/149‑1 (MQ=255)
                                                                                                                                                   aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtg  >  4:1026758/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATG  >  NC_000913/176830‑177125

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