Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,059 C→T 5.1% intergenic (‑126/+16) 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,0590CT 5.1% 444.5 / 10.4 136intergenic (‑126/+16)sseB/pepBrhodanase‑like enzyme, sulfur transfer from thiosulfate/aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base C (64/65);  new base T (2/5);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 142 total aligned reads displayed.

TTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAA  >  NC_000913/2654913‑2655205
                                                                                                                                                  |                                                                                                                                                  
ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                                  >  4:821172/1‑149 (MQ=255)
ttCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                                  >  1:826530/1‑149 (MQ=255)
  cGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                <  2:717289/149‑1 (MQ=255)
     tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                             <  3:195804/149‑1 (MQ=37)
     tttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                             <  1:894160/149‑1 (MQ=255)
      ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                            <  3:1291817/149‑1 (MQ=255)
      ttttGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            <  1:607181/149‑2 (MQ=37)
         tGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatat                                                                                                                                         <  2:290092/149‑1 (MQ=255)
           tttCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTa                                                                                                                                       >  2:1127086/1‑149 (MQ=255)
             tCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGt                                                                                                                                     <  4:333263/149‑1 (MQ=255)
              cGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGtt                                                                                                                                    <  1:794194/149‑1 (MQ=255)
                gACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTa                                                                                                                                  >  2:1399983/1‑149 (MQ=255)
                 aCATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTAc                                                                                                                                 >  3:1230783/1‑149 (MQ=255)
                      aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAg                                                                                                                                      <  2:568797/139‑1 (MQ=255)
                      aGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGt                                                                                                                            >  2:154534/1‑149 (MQ=255)
                       gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                             >  1:1325717/1‑131 (MQ=25)
                       gAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                           >  4:25623/1‑149 (MQ=255)
                         aCCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTaa                                                                                                                         >  2:1469257/1‑149 (MQ=255)
                           cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                                       >  1:159175/1‑149 (MQ=255)
                           cacaATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACa                                                                                                                       >  3:65703/1‑149 (MQ=255)
                              aaTTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                                    >  1:1270747/1‑149 (MQ=255)
                               attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                                   >  1:731673/1‑149 (MQ=255)
                               attcattcAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAtt                                                                                                                   >  2:1023596/1‑149 (MQ=255)
                                ttcattcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTa                                                                                                                  <  2:111259/149‑1 (MQ=255)
                                   attcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCt                                                                                                               <  1:64884/149‑1 (MQ=255)
                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGAt                                                                                                                    <  2:734171/143‑1 (MQ=255)
                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                              >  1:1059520/1‑149 (MQ=255)
                                    ttcaACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTa                                                                                                              >  1:244321/1‑149 (MQ=255)
                                       aaCGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCg                                                                                                           >  2:230968/1‑149 (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)
                                                     aaaTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATGTGTTGAATTTGCATAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGt                                                                                             >  2:930587/1‑149 (MQ=255)
                                                       aTGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACACCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTcc                                                                                           <  2:560683/149‑1 (MQ=255)
                                                        tGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTccc                                                                                          <  1:1405227/149‑1 (MQ=255)
                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       <  4:204116/149‑1 (MQ=255)
                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       >  4:727377/1‑149 (MQ=255)
                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       <  1:1023596/149‑1 (MQ=255)
                                                           cTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       <  3:574798/149‑1 (MQ=255)
                                                            tGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAg                                                                                      >  1:536155/1‑149 (MQ=255)
                                                                 cgcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGccc                                                                                 <  2:87715/149‑1 (MQ=255)
                                                                  gcTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCg                                                                                >  3:1125906/1‑149 (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)
                                                                         ttATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                         >  1:563263/1‑149 (MQ=255)
                                                                          tATCAGGCCTACAAGGGATCCGCAATTTGTTGAGTTTGCAGAATTTGTAGGCCGGCTAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAAAAAATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                        >  2:849998/1‑149 (MQ=38)
                                                                               ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                                   <  2:18283/149‑1 (MQ=255)
                                                                               ggCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTc                                                                   >  4:62371/1‑149 (MQ=255)
                                                                                      aaGTGATTCGCAATATCTTGAATTTGCAGAATTTGTAGACCAGATAAGGCGTTTACGCCGCACCCGGCATAGTCGTTTACGTCGTTAACAGATTAGCTATCGTGCGCAGACCAAGTCCCGACGCGCCCGCAGACCACTGTGCAACCGgc                                                            <  1:1469257/149‑1 (MQ=25)
                                                                                      aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       >  1:1167678/1‑122 (MQ=255)
                                                                                      aaGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTa                                                                                       <  2:1167678/122‑1 (MQ=255)
                                                                                        gggATTCGCAATTTGTTGAATTTGCAGGATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCTGTAGCGCCCGCAGACCACTGTTCAACCGgcgc                                                          >  2:678828/1‑149 (MQ=255)
                                                                                           aTTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCttt                                                       >  2:1083175/1‑149 (MQ=255)
                                                                                            ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                      >  3:416055/1‑149 (MQ=255)
                                                                                             tCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTAc                                                     <  1:348031/149‑1 (MQ=255)
                                                                                                     tgttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAATGTTCAACCGGCGCTTTACGGTAAGTc                                             <  4:416055/149‑1 (MQ=255)
                                                                                                       ttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGc                                           >  2:304322/1‑149 (MQ=255)
                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                         <  3:829022/116‑1 (MQ=255)
                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCa                                                                         >  4:829022/1‑116 (MQ=255)
                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                        <  4:322241/149‑1 (MQ=255)
                                                                                                          aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                        <  3:1332788/149‑1 (MQ=255)
                                                                                                           aTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGTCGCATCCGGCATATTAGTTTACGCCGTTAGCAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCCCCCGCAGACCACTGTTCAACCGTCGCTTTACCGTAAGTCGCCgtg                                       >  4:1261054/1‑147 (MQ=39)
                                                                                                                cAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTc                                  >  2:165823/1‑149 (MQ=255)
                                                                                                                 aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGATCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                                 >  1:1209650/1‑149 (MQ=255)
                                                                                                                 aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCg                                 >  1:1297473/1‑149 (MQ=255)
                                                                                                                   aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGat                               >  3:1392410/1‑149 (MQ=255)
                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                           >  1:685809/1‑100 (MQ=255)
                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGAc                                                                           <  2:685809/100‑1 (MQ=255)
                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                          <  2:536155/149‑1 (MQ=255)
                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcag                         >  4:508793/1‑149 (MQ=255)
                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                           <  4:774529/49‑1 (MQ=255)
                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGtt                                                                                                                           >  3:774529/1‑49 (MQ=255)
                                                                                                                            gCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                      >  2:103443/1‑149 (MQ=255)
                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                      <  3:471464/149‑1 (MQ=255)
                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                      <  4:65703/149‑1 (MQ=255)
                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag                     <  4:1309508/149‑1 (MQ=255)
                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTCAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagccag                     >  3:1024026/1‑149 (MQ=255)
                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCagccagcc                   <  2:159175/149‑1 (MQ=255)
                                                                                                                                gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                          <  4:1106316/141‑1 (MQ=255)
                                                                                                                                gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                          >  3:1106316/1‑141 (MQ=255)
                                                                                                                                gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCt                  <  2:331009/149‑1 (MQ=255)
                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCACGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc               <  1:312381/149‑1 (MQ=255)
                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                        <  2:721951/92‑1 (MQ=255)
                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCAc                                                                        >  1:721951/1‑92 (MQ=255)
                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc               <  2:1297473/149‑1 (MQ=255)
                                                                                                                                     ggCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTg             >  2:375536/1‑149 (MQ=255)
                                                                                                                                      gCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGa            <  2:1270747/149‑1 (MQ=255)
                                                                                                                                        gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc               >  3:772098/1‑144 (MQ=255)
                                                                                                                                        gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc               <  4:772098/144‑1 (MQ=255)
                                                                                                                                        gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa          <  3:84831/149‑1 (MQ=255)
                                                                                                                                        gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa          <  3:465876/149‑1 (MQ=255)
                                                                                                                                           tACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGtt       <  1:1127086/149‑1 (MQ=255)
                                                                                                                                             cgccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTtct     <  3:95754/149‑1 (MQ=255)
                                                                                                                                             cgccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTtct     <  4:503588/149‑1 (MQ=255)
                                                                                                                                              gccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTtctc    >  4:1143526/1‑149 (MQ=255)
                                                                                                                                               ccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTca   >  4:1226029/1‑149 (MQ=255)
                                                                                                                                                cgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTcaa  <  4:1306246/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCGTTTTTTGTTTCGGACATAAGAACCACAATTCATTCAACGTTCGGGGCGCAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAGGGATTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAA  >  NC_000913/2654913‑2655205

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: