Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,052 T→A 5.1% intergenic (‑119/+23) sseB ← / ← pepB rhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,655,0520TA 5.1% 701.2 / 18.0 194intergenic (‑119/+23)sseB/pepBrhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (99/85);  new base A (8/2);  total (107/87)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 202 total aligned reads displayed.

TTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTT  >  NC_000913/2654907‑2655200
                                                                                                                                                 |                                                                                                                                                    
ttCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                  <  2:1606438/150‑1 (MQ=255)
 tCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 <  2:584748/150‑1 (MQ=255)
 tCAAGTTCGTTTTTTGTTTCGGACATAAGAAACACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTTGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                 <  1:2229307/150‑1 (MQ=255)
      ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            >  1:1718229/1‑150 (MQ=255)
      ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            >  2:2481561/1‑150 (MQ=255)
      ttCGTTTTGTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            <  1:2008065/150‑1 (MQ=255)
      gtCGTTTTTAGTGTCGGACAGAAGAACCACAAGTCAGTCACCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                            <  2:1920584/149‑1 (MQ=25)
         gTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                         >  1:1162609/1‑150 (MQ=255)
          ttttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        >  1:3063361/1‑150 (MQ=255)
          ttttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        >  2:1172405/1‑150 (MQ=255)
          ttttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        >  1:295808/1‑150 (MQ=255)
           tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  1:679620/150‑1 (MQ=255)
           tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                       <  1:1587458/150‑1 (MQ=255)
            ttttGTTTCGGACATAAGAACCTCAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCat                                                                                                                                      >  2:1895212/1‑150 (MQ=255)
            ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCat                                                                                                                                      >  1:1270716/1‑150 (MQ=255)
             tttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCata                                                                                                                                     >  2:78029/1‑150 (MQ=255)
              ttGTTTCGGCCAGAAGCACCACAAGGCCGTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGCctac                                                                                                                                    <  2:2490574/150‑2 (MQ=21)
              ttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatat                                                                                                                                    <  2:1822839/150‑1 (MQ=255)
                   tCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGtt                                                                                                                               >  1:1579780/1‑150 (MQ=255)
                      gACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTAc                                                                                                                            >  2:2964549/1‑150 (MQ=255)
                          tAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCg                                                                                                                        >  1:572018/1‑150 (MQ=255)
                          tAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCg                                                                                                                        >  2:1981222/1‑150 (MQ=255)
                          tAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCg                                                                                                                        >  2:1788688/1‑150 (MQ=255)
                           aaGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                       >  1:271153/1‑150 (MQ=255)
                           aaGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                       >  1:270976/1‑150 (MQ=255)
                               aCCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAc                                                                                                                   <  1:2487772/150‑1 (MQ=255)
                                 cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAg                                                                                                                 >  1:2324214/1‑150 (MQ=255)
                                 cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAAATTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAGTAGTTTACGCCGTTAACAg                                                                                                                 >  2:2982468/1‑150 (MQ=255)
                                  acaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGa                                                                                                                >  1:357639/1‑150 (MQ=255)
                                  acaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGa                                                                                                                >  2:2167100/1‑150 (MQ=255)
                                      ttcattcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAg                                                                                                            >  1:943884/1‑150 (MQ=255)
                                         attcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         <  2:1232342/150‑1 (MQ=255)
                                         attcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         <  1:2912158/150‑1 (MQ=255)
                                         attcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                         <  1:2912349/150‑1 (MQ=255)
                                             aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                                     <  2:2739323/150‑1 (MQ=255)
                                                gTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgcg                                                                                                  <  1:2165618/150‑1 (MQ=255)
                                                 ttCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGACTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTGAACAGATTAGCTATCGTgcgc                                                                                                 >  2:384224/1‑150 (MQ=255)
                                                        cgcAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCaa                                                                                          <  1:28066/150‑1 (MQ=255)
                                                        cgcAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCaa                                                                                          <  1:782822/150‑1 (MQ=255)
                                                             aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                                     <  2:2916123/150‑1 (MQ=255)
                                                             aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                                     <  2:2017051/150‑1 (MQ=255)
                                                                ccTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                  <  1:2058700/150‑1 (MQ=255)
                                                                     tGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGcc                                                                             <  1:2462678/150‑1 (MQ=255)
                                                                     tGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGcc                                                                             <  1:1895212/150‑1 (MQ=255)
                                                                        gcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGATTTTGTAGGCCGGTTAAGGCGTTTACGCCGCATCCGGCATATTTGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGc                                                                          <  1:1071179/150‑1 (MQ=255)
                                                                            cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                      <  2:1151413/150‑1 (MQ=255)
                                                                                  tCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGt                                                                >  2:2093187/1‑150 (MQ=255)
                                                                                            caGGGATTCGCAATTTGTTGACTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGgcg                                                      <  2:2621669/149‑1 (MQ=255)
                                                                                              gggATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCt                                                    >  1:3037951/1‑150 (MQ=255)
                                                                                              gggATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCt                                                    >  2:46750/1‑150 (MQ=255)
                                                                                                 aTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                 <  2:3135258/150‑1 (MQ=255)
                                                                                                 aTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                 >  1:1314584/1‑150 (MQ=255)
                                                                                                   tCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                               <  2:2697131/150‑1 (MQ=255)
                                                                                                     gCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGt                                             <  2:3244369/150‑1 (MQ=255)
                                                                                                     gCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGt                                             >  2:1414172/1‑150 (MQ=255)
                                                                                                     gCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGt                                             >  2:884812/1‑150 (MQ=255)
                                                                                                       attttGTTGAATTTGTAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAATCGGCGCTTTACGGTaa                                           <  2:1677421/148‑1 (MQ=255)
                                                                                                       aaTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTaa                                           <  1:2413263/150‑1 (MQ=255)
                                                                                                           tgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCTTTAACAGATTAGCTATCGTGCTCACACCAAGTCCCGTAGCGCCCGCAGACCACTTTTCAACCGGCGCTTTACGGTAAGTCg                                       <  1:3207957/150‑1 (MQ=255)
                                                                                                             ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGcc                                     >  1:72674/1‑150 (MQ=255)
                                                                                                             ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGcc                                     >  1:1547213/1‑150 (MQ=255)
                                                                                                             ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGcc                                     >  1:338049/1‑150 (MQ=255)
                                                                                                               gAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                   <  1:231643/150‑1 (MQ=255)
                                                                                                                aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAg                                  >  1:1385820/1‑150 (MQ=255)
                                                                                                                aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAg                                  >  2:2329191/1‑150 (MQ=255)
                                                                                                                aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAg                                  <  1:1397518/150‑1 (MQ=255)
                                                                                                                aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAg                                  <  1:541101/150‑1 (MQ=255)
                                                                                                                    tGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGGTTTACGGTAAGTCGCCGAGCAGt                              <  2:2459037/150‑1 (MQ=255)
                                                                                                                    tGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                              <  2:2348228/150‑1 (MQ=255)
                                                                                                                    tGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                              <  2:23710/150‑1 (MQ=255)
                                                                                                                    tGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                              <  1:1958018/150‑1 (MQ=255)
                                                                                                                    tGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                              <  1:101425/150‑1 (MQ=255)
                                                                                                                      cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                            <  1:177226/150‑1 (MQ=255)
                                                                                                                       aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGa                           >  2:3020310/1‑150 (MQ=255)
                                                                                                                        gAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGat                          >  2:2183789/1‑150 (MQ=255)
                                                                                                                         aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                         >  1:645229/1‑150 (MQ=255)
                                                                                                                         aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                         >  1:163872/1‑150 (MQ=255)
                                                                                                                          aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGatat                        >  2:1601754/1‑150 (MQ=255)
                                                                                                                            ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGc                      >  1:2894241/1‑150 (MQ=255)
                                                                                                                             tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     >  2:577495/1‑150 (MQ=255)
                                                                                                                             tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     >  2:1133528/1‑150 (MQ=255)
                                                                                                                             tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                     >  2:2670948/1‑150 (MQ=255)
                                                                                                                              gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATGAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcag                    <  1:1841119/150‑1 (MQ=255)
                                                                                                                                 tgCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                 >  2:297004/2‑150 (MQ=255)
                                                                                                                                  gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag                >  1:1452520/1‑150 (MQ=255)
                                                                                                                                   ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccagc               >  2:3197934/1‑150 (MQ=255)
                                                                                                                                   ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccagc               >  2:2543866/1‑150 (MQ=255)
                                                                                                                                    cGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc              >  1:1618796/1‑150 (MQ=255)
                                                                                                                                    cGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGGCGATATGCagccagcc              >  1:2538471/1‑150 (MQ=255)
                                                                                                                                    cGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc              <  1:2748466/150‑1 (MQ=255)
                                                                                                                                    cGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc              <  1:3115337/150‑1 (MQ=255)
                                                                                                                                     ggATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCt             <  2:1523036/150‑1 (MQ=255)
                                                                                                                                       aTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTg           <  2:2593474/150‑1 (MQ=255)
                                                                                                                                       aTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTg           <  2:686782/150‑1 (MQ=255)
                                                                                                                                       aTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTg           <  1:2093187/150‑1 (MQ=255)
                                                                                                                                        tAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCACGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc          <  1:2988460/150‑1 (MQ=255)
                                                                                                                                         aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt         <  2:1273890/150‑1 (MQ=255)
                                                                                                                                         aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt         <  1:2303100/150‑1 (MQ=255)
                                                                                                                                                ttACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGtt  <  1:49834/150‑1 (MQ=255)
                                                                                                                                                ttACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGtt  <  1:49871/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTT  >  NC_000913/2654907‑2655200

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