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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913176,9790CT 5.5% 711.2 / 20.9 219intergenic (+25/+22)erpA/yadSiron‑sulfur cluster insertion protein/UPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (124/81);  new base T (8/4);  total (132/87)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

Only 100 of 221 total aligned reads displayed.

TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128
                                                                                                                                                     |                                                                                                                                                     
     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)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGc                                                                                                                                               <  1:3679471/150‑1 (MQ=255)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                               >  2:1927846/1‑150 (MQ=255)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                               <  1:1962799/150‑1 (MQ=255)
         aagatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTATTCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                              <  2:2202737/147‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGca                                                                                                                                              <  2:212002/150‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                              <  1:2154029/150‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGca                                                                                                                                              >  2:1750984/1‑150 (MQ=255)
           tatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcaca                                                                                                                                            <  1:1304060/150‑1 (MQ=255)
                      cggTTCCGTTGATTATACCGAAGGTCTGGACGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCtt                                                                                                                                 <  2:3613079/150‑1 (MQ=255)
                      cggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATATAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCtt                                                                                                                                 >  1:1645800/1‑150 (MQ=255)
                       ggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTc                                                                                                                                >  2:3849345/1‑150 (MQ=255)
                        gTTCCGTTGATTATAACCAAGGTCTGGAAGGTTCTCGTTTCATAGTGACCAACCCCAAAGCCAAAAGCACCAGCGGTTGCGGTTCTTCCTTTAGTATCTAAACTGGTGTTCATGTGCCGGATAAACAATACGGCACAACTTCTCTCTTCa                                                                                                                               >  1:1402679/1‑150 (MQ=255)
                           acGTTGATTATACCTAAGGTATGGAAGGTTCTCGTTTCAACTTGACAAACCCGAACGCGAAAAGAACCTGCGGTTGCGGTTCTTCCTTTAGTAACTAATCTGTTCTTCATGTGCCGCATAAACTATCCGGCACATCTTCTCTCTTCAACg                                                                                                                            <  2:358430/149‑1 (MQ=255)
                            cGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGTTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGc                                                                                                                           <  2:3686225/150‑1 (MQ=255)
                             gTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGcc                                                                                                                          >  1:2065309/1‑150 (MQ=255)
                              ttGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGccc                                                                                                                         >  1:2496481/1‑150 (MQ=255)
                                 tttATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTTACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGTATAAACTATCGGGCACATCTTCTCTCTTCAACGCCCAtt                                                                                                                      <  2:6298/149‑1 (MQ=255)
                                  ttATACGGAAGGTCTGGAAGGTTCTCGTTTCTTCGTGACCAACCCGATCGCGAAAAGCTCCTGCGGTTGCGGTTCTTCCTTTAGTTTCTAATCTGTTTTTCATGTGCCGGATATACTATCCGGCACATCTTCTCTCTTCACCGCCCATtc                                                                                                                     <  2:136625/150‑1 (MQ=255)
                                  ttATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtc                                                                                                                     >  1:3241054/1‑150 (MQ=255)
                                     taCCGAAGGTATGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCTAAAAGTACCTGGTGTTGCGTTTCTTCCTTTATTATGTAATCTTTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCa                                                                                                                  <  2:732825/150‑1 (MQ=255)
                                      aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGCACAt                                                                                                                                           >  2:1650971/1‑124 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                >  2:1861858/1‑150 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                               >  2:2396615/1‑150 (MQ=255)
                                        cGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATcc                                                                                                               >  2:2252840/1‑150 (MQ=255)
                                         gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                              >  1:1132030/1‑150 (MQ=255)
                                         gAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCa                                                                                                              >  2:3352677/1‑150 (MQ=255)
                                          aaGGTCAGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAg                                                                                                             >  1:954807/1‑150 (MQ=255)
                                           aGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                            <  2:277954/150‑1 (MQ=255)
                                             gtgtGGATGTTTCTGTTTTCATGGTCACCAACCCGACCGCTAAAAGCACCTGGGGTTGGGTTTCTTACTTTAGTATATAATGTTTTGTTCTTGGGCCGTATAAACTATCGGCCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgcgc                                                                                                          <  1:1750984/147‑1 (MQ=255)
                                              tCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCa                                                                                                         >  1:2331048/1‑150 (MQ=255)
                                                 ggAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                      >  2:748182/1‑150 (MQ=255)
                                                        tctcGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCgg                                                                                                                                                >  2:2742276/1‑101 (MQ=255)
                                                            gTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCtt                                                                                           >  1:1003273/1‑150 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAACGCACGTGCGGTTGCGGTTCTTCCTTTCGTATCTAATCTTTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                          <  1:499225/150‑1 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                          >  1:2675266/1‑150 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTa                                                                                          >  1:3080882/1‑150 (MQ=255)
                                                             tttCATCGTGACCAACCCGAACGCGAAAAGCACCTGAGGTTGCGGTTCTTCGTTTAGTATGTAATCTGTTGTTCTTGTGGCGGATAAACTATCCGGCACATCTTCTCTCTTCAACTCCCATTCTCATCCAGCGCAATCGTCTGTAGCTTa                                                                                          >  2:3230235/1‑150 (MQ=255)
                                                              ttCATCGTCACCACCCCCAACGCTAAAAGCACCGGGGGTTGGGGTTCTTCCTTTTGTATCTAATCTTTTTTTCTTGTGCGGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTaa                                                                                         <  1:1116551/150‑1 (MQ=255)
                                                               tCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                        >  2:928740/1‑150 (MQ=255)
                                                               tCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAg                                                                                        <  2:867552/150‑1 (MQ=255)
                                                                   cGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGc                                                                                    <  2:1319786/150‑1 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                   >  1:3399766/1‑150 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                   >  2:837592/1‑150 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCGTCCTGTAGTAGCTAATCTGGTGTTCATGTGCCGGATAAACTATCCGGCACCTCTTCTCTCGGCCACGCCCAGTCTCATCCAGCGCACCCGGCGGTAGCTGACGATGcc                                                                                   >  2:897238/1‑150 (MQ=255)
                                                                         cAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTATTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                              <  2:2445281/150‑1 (MQ=255)
                                                                         cAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGa                                                                              >  2:569309/1‑150 (MQ=255)
                                                                          aaCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGaa                                                                             <  2:38852/150‑1 (MQ=255)
                                                                           aCCCGAACGCTAAAAGCACCTGGGGTTGGGGTTTTTCCTTTTTTATCTAATTTGTTGTTCATGTGCCGGATAAACTACCCGGCAAATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAt                                                                            <  1:367907/150‑1 (MQ=255)
                                                                               gAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCATTTAGTATCTAATCTGTTGTTCTTGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAgcc                                                                        <  2:2062729/150‑1 (MQ=255)
                                                                                   gcgAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCa                                                                    >  1:2239359/1‑150 (MQ=255)
                                                                                          gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAAt                                                             <  1:2261793/150‑1 (MQ=255)
                                                                                          gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCAATCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGaaa                                                             >  1:3691368/1‑149 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAAt                                                         >  1:3270708/1‑150 (MQ=255)
                                                                                              cTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTAACCGGCACATCTTCTCTCTTCAAAGCCCATTATCAACCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAAt                                                         >  1:3440894/1‑150 (MQ=255)
                                                                                               tgcggttgcggtTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATa                                                        >  1:3545576/1‑150 (MQ=255)
                                                                                                gcggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAg                                                       <  2:3040610/150‑1 (MQ=255)
                                                                                                gcggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAg                                                       <  2:2944514/150‑1 (MQ=255)
                                                                                                 gggttggggttCTTCCTTTAGTATCTATTCTGTTGTTCTTGTGCCGGATAAACTATGCGGCAAATCTTCTCTCTTCAACTCCCTTTCTCTTCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGc                                                      <  1:2364258/149‑1 (MQ=255)
                                                                                                    ttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTc                                                   >  2:3868490/1‑150 (MQ=255)
                                                                                                      gcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCAc                                                 <  2:2243583/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)
                                                                                                                   tAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatc                                    <  1:176388/150‑1 (MQ=255)
                                                                                                                    aGTATCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatca                                   <  1:2182034/150‑1 (MQ=255)
                                                                                                                      tATCTATTCTGTTGTTCTTGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATa                                 <  1:2973512/150‑1 (MQ=255)
                                                                                                                       aTCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATAc                                <  1:226110/150‑1 (MQ=255)
                                                                                                                        tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACt                               >  1:2517456/1‑150 (MQ=255)
                                                                                                                        tCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACt                               >  1:2696664/1‑150 (MQ=255)
                                                                                                                          tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTgg                             >  1:1046231/1‑150 (MQ=255)
                                                                                                                          tAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGACCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCCATAGCGTCACGACCCTGCCCATCATACTgg                             >  2:3834904/1‑150 (MQ=255)
                                                                                                                           aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                            >  2:527408/1‑150 (MQ=255)
                                                                                                                             tCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                          <  1:1484459/150‑1 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                        >  2:2831233/1‑150 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                        >  1:2582805/1‑150 (MQ=255)
                                                                                                                                 ttgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                      <  1:3911453/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    >  2:2294066/1‑150 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  1:1556430/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  1:389478/150‑1 (MQ=255)
                                                                                                                                   gttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCa                    <  1:652650/150‑1 (MQ=255)
                                                                                                                                         gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTa              <  1:3725658/150‑1 (MQ=255)
                                                                                                                                          tgCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             >  2:1979848/1‑150 (MQ=255)
                                                                                                                                          tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             <  2:3000850/150‑1 (MQ=255)
                                                                                                                                          tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             >  1:1994828/1‑150 (MQ=255)
                                                                                                                                          tgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTAc             <  2:2536353/150‑1 (MQ=255)
                                                                                                                                            ccGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACgg           >  2:1533480/1‑150 (MQ=255)
                                                                                                                                                aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaaa       <  1:3047679/150‑1 (MQ=255)
                                                                                                                                                aTAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGaaaa       >  2:231724/1‑150 (MQ=255)
                                                                                                                                                 tAAACTATCCGGAACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAACCGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt      <  1:3635878/150‑1 (MQ=255)
                                                                                                                                                   aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgt    <  1:390457/150‑1 (MQ=255)
                                                                                                                                                   aaCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgt    <  2:3080882/150‑1 (MQ=255)
                                                                                                                                                    aCTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   >  2:2609120/1‑150 (MQ=255)
                                                                                                                                                    aCTATCCGGCAAATCTTCTCTCTTCAACGCCCAGTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAtgtg   <  2:1471899/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128

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