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% 699.8 / 22.1 212intergenic (+23/+24)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (120/78);  new base G (8/4);  total (129/83)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 217 total aligned reads displayed.

ATGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGT  >  NC_000913/176829‑177126
                                                                                                                                                    |                                                                                                                                                     
aTGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATaaa                                                                                                                                                      <  1:2279878/150‑1 (MQ=255)
 tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                     <  1:2088999/150‑1 (MQ=255)
    gCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTAt                                                                                                                                                  >  1:2536353/1‑150 (MQ=255)
      cTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                                >  1:3906955/1‑150 (MQ=255)
      cTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATATGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                                >  1:889718/1‑150 (MQ=255)
        gCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTACTCTGTTGTTCATGTGCCGGATAAACTATCCgg                                                                                                                                              <  2:2562797/150‑1 (MQ=255)
        gCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCgg                                                                                                                                              >  2:2190436/1‑150 (MQ=255)
        gCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCgg                                                                                                                                              <  2:397906/150‑1 (MQ=255)
         cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGCTAAACTATCCGGc                                                                                                                                             >  1:3211072/1‑150 (MQ=255)
         cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                             <  1:1962799/150‑1 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGca                                                                                                                                            <  2:212002/150‑1 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                            <  2:1719356/150‑1 (MQ=255)
          aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                            <  2:3653332/150‑1 (MQ=255)
                    cggcggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctct                                                                                                                                  <  2:1931387/150‑1 (MQ=255)
                    cggcggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctct                                                                                                                                  <  2:522263/150‑1 (MQ=255)
                       cggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATATAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCtt                                                                                                                               >  1:1645800/1‑150 (MQ=255)
                        ggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTc                                                                                                                              >  2:3069252/1‑150 (MQ=255)
                         gTTCCGTTGATTATAACCAAGGTCTGGAAGGTTCTCGTTTCATAGTGACCAACCCCAAAGCCAAAAGCACCAGCGGTTGCGGTTCTTCCTTTAGTATCTAAACTGGTGTTCATGTGCCGGATAAACAATACGGCACAACTTCTCTCTTCa                                                                                                                             >  1:1402679/1‑150 (MQ=255)
                              gTTGATTATACCGAATGTCTGGAATGTTCTCGTATCATCGTGACCAACACGAACGCGAAAAGCACCTGCGGGTGCGGTACTTCCTTTAGTATCTAATCTGTATTTCATGTGCAGGATAAACAATACGGCGCATCTTCTCACATCAACGcc                                                                                                                        >  2:2904226/1‑150 (MQ=255)
                              gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                        >  1:490056/1‑150 (MQ=255)
                              gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                        >  2:2077720/1‑150 (MQ=255)
                              gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                        >  1:2065309/1‑150 (MQ=255)
                               ttGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGccc                                                                                                                       >  1:2496481/1‑150 (MQ=255)
                                  aTTATACCGAAGGTCTGGAAGGTACTCGTTTCATCGTGACAAACCCCAACGCGAAAACCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTACATGTGCCGGATAAACTAACCGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                    <  2:3726547/150‑1 (MQ=255)
                                   ttATACGGAAGGTCTGGAAGGTTCTCGTTTCTTCGTGACCAACCCGATCGCGAAAAGCTCCTGCGGTTGCGGTTCTTCCTTTAGTTTCTAATCTGTTTTTCATGTGCCGGATATACTATCCGGCACATCTTCTCTCTTCACCGCCCATtc                                                                                                                   <  2:136625/150‑1 (MQ=255)
                                   ttATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtc                                                                                                                   >  1:3241054/1‑150 (MQ=255)
                                       aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGCACAt                                                                                                                                         <  1:1650971/124‑1 (MQ=255)
                                         cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                             >  2:2396615/1‑150 (MQ=255)
                                          gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                            >  1:1132030/1‑150 (MQ=255)
                                          gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                            <  2:2752847/150‑1 (MQ=255)
                                           aaGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAg                                                                                                           >  1:2175828/1‑150 (MQ=255)
                                           aaGGTCTGGAAGGGTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAg                                                                                                           <  2:279802/150‑1 (MQ=255)
                                           aaGGTCAGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAg                                                                                                           >  1:954807/1‑150 (MQ=255)
                                            aGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGTTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCAAATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                          <  2:3113755/150‑1 (MQ=255)
                                              gTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgcgc                                                                                                        >  1:3895800/1‑150 (MQ=255)
                                               tCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCa                                                                                                       >  1:2331048/1‑150 (MQ=255)
                                                 tGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCAACCAGCGCaaa                                                                                                     >  1:1243883/1‑150 (MQ=255)
                                                  ggAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                    >  2:748182/1‑150 (MQ=255)
                                                           tcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGc                                                                                           >  2:2542377/1‑150 (MQ=255)
                                                           tcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGc                                                                                           >  2:1754990/1‑150 (MQ=255)
                                                             gTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCtt                                                                                         >  1:1003273/1‑150 (MQ=255)
                                                              tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                        >  1:3080882/1‑150 (MQ=255)
                                                              tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                        >  1:2675266/1‑150 (MQ=255)
                                                              tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGAGGTTGCGGTTCTTCGTTTAGTATGTAATCTGTTGTTCTTGTGGCGGATAAACTATCCGGCACATCTTCTCTCTTCAACTCCCATTCTCATCCAGCGCAATCGTCTGTAGCTTa                                                                                        >  2:3230235/1‑150 (MQ=255)
                                                               ttCATCGTCACCACCCCCAACGCTAAAAGCACCGGGGGTTGGGGTTCTTCCTTTTGTATCTAATCTTTTTTTCTTGTGCGGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTaa                                                                                       <  1:1116551/150‑1 (MQ=255)
                                                                tCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                      >  2:2244624/1‑150 (MQ=255)
                                                                tCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                      <  2:867552/150‑1 (MQ=255)
                                                                tCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                      >  2:928740/1‑150 (MQ=255)
                                                                     gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                 >  2:837534/1‑150 (MQ=255)
                                                                     gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                 >  1:3399766/1‑150 (MQ=255)
                                                                     gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCGTCCTGTAGTAGCTAATCTGGTGTTCATGTGCCGGATAAACTATCCGGCACCTCTTCTCTCGGCCACGCCCAGTCTCATCCAGCGCACCCGGCGGTAGCTGACGATGcc                                                                                 >  2:897238/1‑150 (MQ=255)
                                                                         ccAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACg                                                                             >  1:3654262/1‑150 (MQ=255)
                                                                           aaCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                           <  2:38852/150‑1 (MQ=255)
                                                                            ccccGAACGCTAAAAGCCCCTGCGGTTGCGTTTCTTCCTTTCTTATATAATCTTTTTTTCTTGTGCCGTATAAACTAGCCGGCAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAt                                                                          <  1:367944/149‑1 (MQ=255)
                                                                              ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAg                                                                        >  2:2266697/1‑150 (MQ=255)
                                                                              ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAg                                                                        >  2:253604/1‑150 (MQ=255)
                                                                              ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAg                                                                        >  1:3852687/1‑150 (MQ=255)
                                                                                gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                      <  2:2062678/150‑1 (MQ=255)
                                                                                      gacaaGCACCTGCGGTTGCGGTTCTTCCTTTAGTCTCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGc                                                                <  1:2220291/147‑1 (MQ=255)
                                                                                         aaGCACCTGGGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGa                                                             <  2:710790/150‑1 (MQ=255)
                                                                                           gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAAt                                                           <  1:2261793/150‑1 (MQ=255)
                                                                                           gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAATCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaaa                                                           >  1:3691368/1‑149 (MQ=255)
                                                                                            cACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATc                                                          >  2:3521619/1‑150 (MQ=255)
                                                                                               cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTAACCGGCACATCTTCTCTCTTCAAAGCCCATTATCAACCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAAt                                                       >  1:3440894/1‑150 (MQ=255)
                                                                                                tgcggttgcggtTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATa                                                      >  1:3545576/1‑150 (MQ=255)
                                                                                                  gggttggggttCTTCCTTTAGTATCTATTCTGTTGTTCTTGTGCCGGATAAACTATGCGGCAAATCTTCTCTCTTCAACTCCCTTTCTCTTCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGc                                                    <  1:2364258/149‑1 (MQ=255)
                                                                                                     ttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTc                                                 >  2:3868351/1‑150 (MQ=255)
                                                                                                      tgcgtttCTTCATTTAGTATCTAATCTTTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCa                                                <  2:2675266/150‑1 (MQ=255)
                                                                                                              ttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATg                                        >  1:2898446/1‑150 (MQ=255)
                                                                                                               tcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                       >  2:1109138/1‑150 (MQ=255)
                                                                                                               tcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTACTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                       >  1:1211555/1‑150 (MQ=255)
                                                                                                                  tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                    >  2:102696/1‑150 (MQ=255)
                                                                                                                  tttAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCca                                    >  1:2581461/1‑150 (MQ=255)
                                                                                                                       tATCTATTCTGTTGTTCTTGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATa                               <  1:2973512/150‑1 (MQ=255)
                                                                                                                        aTCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATAc                              <  1:226110/150‑1 (MQ=255)
                                                                                                                         tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACt                             >  1:2696664/1‑150 (MQ=255)
                                                                                                                         tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCAACATACt                             >  2:3590674/1‑150 (MQ=255)
                                                                                                                           tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTgg                           >  1:1046231/1‑150 (MQ=255)
                                                                                                                           tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTgg                           >  2:132119/1‑150 (MQ=255)
                                                                                                                           tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGACCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCCATAGCGTCACGACCCTGCCCATCATACTgg                           >  2:3834904/1‑150 (MQ=255)
                                                                                                                            aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                          >  2:527413/1‑150 (MQ=255)
                                                                                                                            aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                          >  2:1526191/1‑150 (MQ=255)
                                                                                                                             aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCt                         >  2:2961383/1‑150 (MQ=255)
                                                                                                                              tCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                        <  1:1484459/150‑1 (MQ=255)
                                                                                                                                tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                      >  1:2582805/1‑150 (MQ=255)
                                                                                                                                    gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                  <  1:652650/150‑1 (MQ=255)
                                                                                                                                    gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                  >  2:2294066/1‑150 (MQ=255)
                                                                                                                                    gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                  <  1:1556430/150‑1 (MQ=255)
                                                                                                                                     ttCTTGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                 <  2:1536130/150‑1 (MQ=255)
                                                                                                                                       cATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTg               <  2:2207016/150‑1 (MQ=255)
                                                                                                                                          gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa            <  1:3725658/150‑1 (MQ=255)
                                                                                                                                           tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc           >  1:1994828/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg         >  2:3034527/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg         >  1:3264616/1‑150 (MQ=255)
                                                                                                                                               ggATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGACTCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaa       >  2:1923890/1‑150 (MQ=255)
                                                                                                                                                 aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaaa     <  1:3047679/150‑1 (MQ=255)
                                                                                                                                                 aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaaa     >  2:231724/1‑150 (MQ=255)
                                                                                                                                                  tAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt    >  2:571107/1‑150 (MQ=255)
                                                                                                                                                  tAAACTATCCGGAACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAACCGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt    <  1:3635878/150‑1 (MQ=255)
                                                                                                                                                    aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgt  <  1:1754990/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGT  >  NC_000913/176829‑177126

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: