Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,051 T→A 5.7% intergenic (‑118/+24) 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,0510TA 5.7% 493.5 / 14.6 141intergenic (‑118/+24)sseB/pepBrhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (66/67);  new base A (2/6);  total (68/73)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 143 total aligned reads displayed.

GTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATA  >  NC_000913/2654904‑2655197
                                                                                                                                                   |                                                                                                                                                  
gTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                                   >  1:17883/1‑149 (MQ=255)
gTCTTCAAGTTCGTTTTTTGTATCGGACAAAAGAATCACAATTCATTCAACGTTCGGGGCGCAAATGCATGATGCGCTAAGCTTAACAGTCCTACAAGGGATTCGCAATTTGATGAATTGGCTTAATTTGTAGGCCGGATAAGGCGttt                                                                                                                                                   >  4:905145/1‑149 (MQ=37)
  cTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                 <  4:211006/149‑1 (MQ=255)
   ttCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                >  2:794194/1‑149 (MQ=255)
     cAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                              <  2:405680/149‑1 (MQ=255)
         ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTATGCCGGATGCGGCGTAAAcgccttat                                                                                                                                          <  3:302631/149‑5 (MQ=38)
         ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                          >  1:826530/1‑149 (MQ=255)
           cGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                        <  2:717289/149‑1 (MQ=255)
               ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                    <  3:1291817/149‑1 (MQ=255)
               ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                    <  1:607181/149‑2 (MQ=37)
               ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                    <  2:560497/149‑2 (MQ=37)
                  tGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatat                                                                                                                                 <  2:290092/149‑1 (MQ=255)
                      tCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGt                                                                                                                             <  4:333263/149‑1 (MQ=255)
                       cGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGtt                                                                                                                            <  1:794194/149‑1 (MQ=255)
                         gACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTa                                                                                                                          >  2:1399983/1‑149 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                              >  1:568797/1‑139 (MQ=255)
                               aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                    <  1:865264/149‑1 (MQ=255)
                                gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                     >  1:1325717/1‑131 (MQ=25)
                                  aCCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTaa                                                                                                                 >  2:1469257/1‑149 (MQ=255)
                                   ccACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAc                                                                                                                >  2:312381/1‑149 (MQ=255)
                                    cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                               <  3:1196506/149‑1 (MQ=255)
                                    cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                               >  3:65703/1‑149 (MQ=255)
                                     acaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAg                                                                                                              >  4:1378559/1‑149 (MQ=255)
                                       aaTTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                            <  3:130869/149‑1 (MQ=255)
                                        attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                           >  1:731673/1‑149 (MQ=255)
                                        attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                           >  2:1023596/1‑149 (MQ=255)
                                         ttcattcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTa                                                                                                          <  2:111259/149‑1 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                            >  1:734171/1‑143 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                      >  1:1059520/1‑149 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                      >  2:861345/1‑149 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                      >  4:583635/1‑149 (MQ=255)
                                             ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                      >  1:244321/1‑149 (MQ=255)
                                                aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCg                                                                                                   >  2:230968/1‑149 (MQ=255)
                                                 aCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                                  >  1:463612/1‑149 (MQ=255)
                                                 aCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGt                                                                                                  >  3:503588/1‑149 (MQ=255)
                                                   gTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgc                                                                                                <  2:590217/149‑1 (MQ=255)
                                                            gcAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCaa                                                                                       <  4:1230783/149‑1 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                     >  4:214319/1‑101 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                     <  3:214319/101‑1 (MQ=255)
                                                              aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGt                                                                                     >  4:411740/1‑149 (MQ=255)
                                                                aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACACCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTcc                                                                                   <  2:560683/149‑1 (MQ=255)
                                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                                  <  1:1405227/149‑1 (MQ=255)
                                                                    cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                               <  1:1023596/149‑1 (MQ=255)
                                                                    cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                               >  4:727377/1‑149 (MQ=255)
                                                                     tGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAg                                                                              >  1:536155/1‑149 (MQ=255)
                                                                          cgcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGccc                                                                         <  2:87715/149‑1 (MQ=255)
                                                                           gcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCg                                                                        >  3:1125906/1‑149 (MQ=255)
                                                                            cTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGc                                                                       <  4:733961/149‑1 (MQ=255)
                                                                            cTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGc                                                                       >  1:412566/1‑149 (MQ=255)
                                                                               cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGa                                                                    <  3:1424702/149‑1 (MQ=255)
                                                                               cGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGa                                                                    <  2:850591/149‑1 (MQ=255)
                                                                                gCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                   >  2:1312341/1‑149 (MQ=255)
                                                                                       aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgc                                                                                                >  1:1265120/1‑113 (MQ=255)
                                                                                       aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTgc                                                                                                <  2:1265120/113‑1 (MQ=255)
                                                                                       aGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGtt                                                            >  4:2351/1‑149 (MQ=255)
                                                                                        ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                           >  4:62371/1‑149 (MQ=255)
                                                                                        ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                           <  2:18283/149‑1 (MQ=255)
                                                                                               aaGTGATTCGCAATATCTTGAATTTGCAGAATTTGTAGACCAGATAAGGCGTTTACGCCGCACCCGGCATAGTCGTTTACGTCGTTAACAGATTAGCTATCGTGCGCAGACCAAGTCCCGACGCGCCCGCAGACCACTGTGCAACCGgc                                                    <  1:1469257/149‑1 (MQ=25)
                                                                                               aaGGGCTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGgc                                                    >  4:796089/1‑149 (MQ=255)
                                                                                               aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                               <  2:1167678/122‑1 (MQ=255)
                                                                                                 gggATTCGCAATTTGTTGAATTTGCAGGATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCTGTAGCGCCCGCAGACCACTGTTCAACCGgcgc                                                  >  2:678828/1‑149 (MQ=255)
                                                                                                    aTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCttt                                               >  2:1083175/1‑149 (MQ=255)
                                                                                                           aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                         >  3:1060801/1‑132 (MQ=255)
                                                                                                           aTTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCaa                                                         <  4:1060801/132‑1 (MQ=255)
                                                                                                              tgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAATGTTCAACCGGCGCTTTACGGTAAGTc                                     <  4:416055/149‑1 (MQ=255)
                                                                                                                ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGc                                   >  2:304322/1‑149 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                 >  4:829022/1‑116 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                <  4:322241/149‑1 (MQ=255)
                                                                                                                   aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                <  3:1332788/149‑1 (MQ=255)
                                                                                                                         cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                          >  2:165823/1‑149 (MQ=255)
                                                                                                                         cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                          <  3:727377/149‑1 (MQ=255)
                                                                                                                          aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGATCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                         >  1:1209650/1‑149 (MQ=255)
                                                                                                                          aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                         >  1:1297473/1‑149 (MQ=255)
                                                                                                                             aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                      <  2:112640/149‑1 (MQ=255)
                                                                                                                             aTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGata                      <  1:304322/149‑1 (MQ=255)
                                                                                                                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGatat                     >  2:1482221/1‑149 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                   <  2:685809/100‑1 (MQ=255)
                                                                                                                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                   >  1:685809/1‑100 (MQ=255)
                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcag                 >  4:508793/1‑149 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                   >  3:774529/1‑49 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                   <  4:774529/49‑1 (MQ=255)
                                                                                                                                     gCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca              >  2:103443/1‑149 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca              <  3:471464/149‑1 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag             <  4:1309508/149‑1 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTCAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag             >  3:1024026/1‑149 (MQ=255)
                                                                                                                                        ggatgcggCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTaa                                        <  3:152599/114‑1 (MQ=38)
                                                                                                                                        ggATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc           <  2:159175/149‑1 (MQ=255)
                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                  <  4:1106316/141‑1 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCACGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc       <  1:312381/149‑1 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                <  2:721951/92‑1 (MQ=255)
                                                                                                                                            aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc       <  2:1297473/149‑1 (MQ=255)
                                                                                                                                             aGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt      <  3:821172/149‑1 (MQ=255)
                                                                                                                                              ggCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTg     >  2:375536/1‑149 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGa    <  1:714731/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  <  3:84831/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  <  1:1166487/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  >  4:1245062/1‑149 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  <  3:465876/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  <  1:154534/149‑1 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa  <  2:81939/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTTCAAGTTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATA  >  NC_000913/2654904‑2655197

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