Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 176,978 A→T 5.1% intergenic (+24/+23) 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,9780AT 5.1% 1042.2 / 32.3 304intergenic (+24/+23)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (157/128);  new base T (5/11);  total (162/142)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

Only 100 of 323 total aligned reads displayed.

ATGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTG  >  NC_000913/176829‑177127
                                                                                                                                                     |                                                                                                                                                     
aTGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                       <  2:5142817/150‑1 (MQ=255)
aTGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                       <  2:369933/150‑1 (MQ=255)
 tGAGCCTGCAATATCTTGTCGGCGGTTCCGTTTATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGTATaaaa                                                                                                                                                      <  1:4366777/150‑2 (MQ=255)
  gAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACt                                                                                                                                                     >  1:5987236/1‑150 (MQ=255)
     ccTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATc                                                                                                                                                  >  1:4763844/1‑150 (MQ=255)
     ccTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATc                                                                                                                                                  >  1:526507/1‑150 (MQ=255)
       tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                >  2:1662514/1‑150 (MQ=255)
       tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                >  1:4185848/1‑150 (MQ=255)
       tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                >  2:2250884/1‑150 (MQ=255)
         cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                              >  2:4887499/1‑150 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATAGCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGca                                                                                                                                             <  2:2704492/150‑1 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGca                                                                                                                                             <  1:5388840/150‑1 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                             >  2:4167085/1‑150 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                             >  1:2995892/1‑150 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGGTGCGGTTCGTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGcc                                                                                                                                             >  1:5419564/1‑149 (MQ=255)
                tGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGGGGTTCTTCCTTTAGTATCTAATCTGTTGTTCTTGTGCCGGATAAACTATCCGGCACAtctt                                                                                                                                       <  1:488371/150‑1 (MQ=255)
                tGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtctt                                                                                                                                       <  1:3985138/150‑1 (MQ=255)
                   tCGGCGGTTCCGTTGATTATAACGAAGGTATGGAAGGTTCTCGTTTCAAAGTTAACAAACCCAAAGCGAAAATCAACAGCGGTTGCGGTTCTTCCGTTAGTATATAATATGTTTTTAATTTGCCGGATAAAATAACAGGCCCAACTtctc                                                                                                                                    >  1:1952021/1‑150 (MQ=255)
                              gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                         <  2:1435197/150‑1 (MQ=255)
                                    tataCCGAAGGTCTGGCAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCATTTAGTATCTATTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                   <  2:3044272/150‑1 (MQ=255)
                                    tataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                   >  2:3950062/1‑150 (MQ=255)
                                     ataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCTAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtctc                                                                                                                  <  2:5347634/150‑1 (MQ=255)
                                     ataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtctc                                                                                                                  <  2:3652114/150‑1 (MQ=255)
                                       aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCAt                                                                                                                >  1:3355735/1‑150 (MQ=255)
                                        ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                               >  1:126801/1‑150 (MQ=255)
                                         cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                              >  1:461910/1‑150 (MQ=255)
                                         cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                              >  1:785553/1‑150 (MQ=255)
                                          gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                             >  2:1873495/1‑150 (MQ=255)
                                                   gAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACg                                                                                                    >  2:2862230/1‑150 (MQ=255)
                                                      gtttCTCGTTAAATCGTGACCAACCCGAACGCTAAATGCACCTGGGGTTGGGTTTTTTCCTTTTTTATCTAATGTTTTTTTCTTGTGCCGTATAAACTAGCCGGAAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                                 <  1:2039101/148‑1 (MQ=255)
                                                      ggTTATCGTTTCTTCGTGACCAACCCGAACGCTAAAAGAACGTGGGGTTGGTGTTTTTCTTTTATTTTCTATTGTTTTGTTCTTGTGCCTTATAAACTCTCCGCAACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCg                                                                                                 <  1:1549556/150‑1 (MQ=255)
                                                       gTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCgg                                                                                                >  2:5637694/1‑150 (MQ=255)
                                                        ttCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGt                                                                                               >  1:1125497/1‑150 (MQ=255)
                                                          ctcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAg                                                                                             >  1:2854559/1‑150 (MQ=255)
                                                              tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCtga                                                                                         >  2:3749148/1‑148 (MQ=255)
                                                               ttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTaa                                                                                        <  2:614316/150‑1 (MQ=255)
                                                               tacATCGTTACCAACCCCAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCGGGATAAACTATCCGCCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTaa                                                                                        <  1:3662272/148‑1 (MQ=255)
                                                                 cATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGa                                                                                      <  2:2480968/150‑1 (MQ=255)
                                                                    ggTGACCAACCCTAACGCTAAAAGCACCTGCGGTTGGGTTTCTTCATTTTTTATCTAATCTGTTTTTCATGTGCCGTATAAACTAACCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                   <  1:1423612/149‑1 (MQ=255)
                                                                     gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                  >  1:248232/1‑150 (MQ=255)
                                                                       gACCACCCCGACCGCGAACCGCACGTGGGGTGGGGGTTCTTCCTTTTGTCTCTAATCGTTTGTTCTTGGGCGGGATAAACTATCGGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCaa                                                                                <  1:4693744/150‑1 (MQ=255)
                                                                        accaccCCGACCGCTAAAAGCACCGGGGGTTGGGGTTCTTCCTTTTGTATCAAATCTGTTGTCCTTGTGCCGCATAAACTACCGGCCAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAAc                                                                               <  1:2541698/150‑1 (MQ=255)
                                                                        aCCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAAc                                                                               <  2:1269194/150‑1 (MQ=255)
                                                                          cAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                             >  1:933427/1‑150 (MQ=255)
                                                                          cAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAGTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                             >  1:964484/1‑150 (MQ=255)
                                                                             cccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATa                                                                          <  2:709857/150‑1 (MQ=255)
                                                                                gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTGCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                       >  1:557067/1‑150 (MQ=255)
                                                                                       aaaaGCACCTGCGGTTGCGGTTCTTCCTTTAGTATATAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGcc                                                                <  1:5726531/150‑1 (MQ=255)
                                                                                        aaaGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCg                                                               <  1:5886836/150‑1 (MQ=255)
                                                                                         aaGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGa                                                              >  1:4940144/1‑150 (MQ=255)
                                                                                          aGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaa                                                             >  2:4093977/1‑150 (MQ=255)
                                                                                            cACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATc                                                           <  2:4789374/150‑1 (MQ=255)
                                                                                                  cggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCACCGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGc                                                     <  2:3232463/150‑1 (MQ=255)
                                                                                                     ttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTc                                                  >  1:5376706/1‑150 (MQ=255)
                                                                                                          gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAc                                             >  1:841951/1‑150 (MQ=255)
                                                                                                           ttcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAcc                                            >  2:4420963/1‑150 (MQ=255)
                                                                                                           ttcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAcc                                            >  2:2523038/1‑150 (MQ=255)
                                                                                                            ttttcTTTTTGTATATATTGTTTTTTTCTTGTGCCGTATAAAATATCCGGAAAATCTTCTCTATTCAACGCCCTTTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCa                                           <  1:4071971/148‑1 (MQ=255)
                                                                                                             cttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCAt                                          <  2:3481773/150‑1 (MQ=255)
                                                                                                              ttcCTTTATTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATg                                         <  1:122086/150‑1 (MQ=255)
                                                                                                               tcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                        >  2:4231240/1‑150 (MQ=255)
                                                                                                                 cTTTAGTATCTATTCTGTGGTTCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGccc                                      <  1:5814634/150‑1 (MQ=255)
                                                                                                                  tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                     >  1:1835899/1‑150 (MQ=255)
                                                                                                                  tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                     >  1:1027221/1‑150 (MQ=255)
                                                                                                                    tAGTATCTAATCTGTTGTTCTTGTGCCGTATAACCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatc                                   <  2:2130274/150‑1 (MQ=255)
                                                                                                                      gTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatcat                                 >  1:4325446/1‑150 (MQ=255)
                                                                                                                       tATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATa                                >  1:4082934/1‑150 (MQ=255)
                                                                                                                         tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACt                              >  2:3450487/1‑150 (MQ=255)
                                                                                                                         tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCAACATACt                              >  1:1885907/1‑150 (MQ=255)
                                                                                                                          cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                             <  2:595982/150‑1 (MQ=255)
                                                                                                                            aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                           <  2:5444272/150‑1 (MQ=255)
                                                                                                                            aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                           >  1:3847321/1‑150 (MQ=255)
                                                                                                                             aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCt                          <  2:790198/150‑1 (MQ=255)
                                                                                                                             aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCt                          <  1:4534352/150‑1 (MQ=255)
                                                                                                                               cTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                        <  1:2612334/150‑1 (MQ=255)
                                                                                                                               cTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGt                        <  1:4103323/150‑1 (MQ=255)
                                                                                                                                 gttgttCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGttt                      <  2:4533913/150‑1 (MQ=255)
                                                                                                                                  ttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                     >  2:2385737/1‑150 (MQ=255)
                                                                                                                                   tgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTcc                    <  1:5654405/150‑1 (MQ=255)
                                                                                                                                    gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                   >  1:3437254/1‑150 (MQ=255)
                                                                                                                                    gttCATGGGCCGTATAAACTTTCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCTCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGATTAAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                   <  1:545917/150‑1 (MQ=255)
                                                                                                                                     ttCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                  <  1:5088516/150‑1 (MQ=255)
                                                                                                                                     ttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                  >  1:5919714/1‑150 (MQ=255)
                                                                                                                                          gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa             <  1:2376867/150‑1 (MQ=255)
                                                                                                                                           tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc            >  2:4444334/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg           <  2:933427/150‑1 (MQ=255)
                                                                                                                                            gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg           >  1:4240625/1‑150 (MQ=255)
                                                                                                                                             ccGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg          >  2:4063354/1‑150 (MQ=255)
                                                                                                                                             ccGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg          >  2:1049450/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCa                                                                                                        <  2:5544556/56‑1 (MQ=255)
                                                                                                                                               ggATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaa        >  1:2161991/1‑150 (MQ=255)
                                                                                                                                                gATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCCTTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaa       >  1:3598155/1‑150 (MQ=255)
                                                                                                                                                gATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAACCGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaa       <  2:2935664/150‑1 (MQ=255)
                                                                                                                                                 aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGccc                                      >  2:258722/1‑118 (MQ=255)
                                                                                                                                                    aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgt   >  2:3319797/1‑150 (MQ=255)
                                                                                                                                                    aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgt   <  2:1544774/150‑1 (MQ=255)
                                                                                                                                                     aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg  <  2:3821121/150‑1 (MQ=255)
                                                                                                                                                     aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg  <  1:1780748/150‑1 (MQ=255)
                                                                                                                                                     aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg  <  2:3401626/150‑1 (MQ=255)
                                                                                                                                                     aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg  <  1:5344160/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTG  >  NC_000913/176829‑177127

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: